Tfee Solar Gell 
RADIATION 
HANDBOOK 



SOLAR CELL RADIATION- HANDBOOK 


28 June 1973 


Prepared by: 


J. R. Carter, 
H, Y, Tada 


21945-6001- RU- 00 


Contract No. ^953362 

Jet Propulsion Laboratory 
California Institute of Technology 
Pasadena, California 


Approved by: 


Jr. 



TRW 

SYSTEMS GROUP 

ONE SPACE PARK ■ REOONDO BEACH CALIFORNIA 90278 


This work was performed for the Jet Propulsion Laboratory, 
California Institute of Technology, as sponsored by the National 
Aeronautics and Space Administration under NAS 7-100. 


ACKNOWLEDGEMENT 


This document was conceived and supported by the 
Jet Propulsion Laboratory, California Institute of 
Technology. In this connection, the efforts of Dr. 
Bruce E. Anspaugh, project manager, and Mr. John V. 
Goldsmith are gratefully acknowledged. 

The preparation of the document was done within 
the Vulnerability and Hardness Laboratory, TRW Systems. 
In this regard, the labors of many members of this 
organization are appreciated. In addition, Mr. Werner 
Luft, Space Vehicles Division, TRW Systems, provided 
essential unpublished data. 


FOREWORD 


The purpose of this document is to detail a method of predicting 
the degradation of a solar array in a space radiation environment. The 
text contains a discussion of solar cell technology which emphasizes 
the cell parameters which degrade in a radiation environment. The 
experimental techniques used in the evaluation of radiation effects 
are discussed. In Chapter 3, the theoretical aspects of radiation 
damage are discussed, and the experimental data, on which the concept 
of damage equivalent 1-MeV electron fluence is based, are presented. 

In Chapter 4, the methods of developing relative damage coefficients 
from the experimental data are detailed. In this regard, it was found 
necessary to institute two separate equivalent fluences to properly 
describe the changes of solar cell parameters under space proton irra- 
diation. 

Chapter 5 concerns the nature of the space radiation environment 
and contains predicted solar flare proton fluences for the twenty-first 
solar cycle based on a proposed model. In Chapter 6, the method of 
calculating equivalent fluence from electron and proton energy spectrums 
and relative damage coefficients is detailed. In addition, computer- 
calculated equivalent fluence contributions from trapped electrons and 
protons are tabulated for an extensive series of circular earth orbits. 
The estimated annual equivalent fluence contributions due to solar flare 
protons are tabulated for the remainder of the current solar cycle. 

The estimation of degraded solar cell output characteristics from equi- 
valent fluence values is discussed. In Chapter 7, flight data from 
satellites is compared with estimated degradation. 
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CHAPTER 1 


1.0 THEORY OF THE SILICON SOLAR CELL 

In this chapter some elementary concepts of semiconductor theory, 
which are useful in understanding solar cell operation, are described. 

The operation of the silicon solar cell is discussed in terms of an 
equivalent circuit, and the electrical characteristics of the equivalent 
circuit elements are explained in terms of physical quantities. In 
addition, the physical structure of a silicon solar cell is detailed. 

1 . 1 Semiconductor Theory 

1.1.1 Thermal Equilibrium Relationships and Excess Densities . Semicon- 
ductors are a class of materials which have electrical properties and 
physical characteristics intermediate between metals and dielectrics. 

An important characteristic of semiconductor materials is bipolar con- 
duction, where charge transport may occur by conduction band electrons 
or through empty energy states in the valence band which behave electri- 
cally like positively-charged electrons and are referred to as holes. 

The equilibrium concentrations of conduction electrons and holes in 
silicon are determined from thermal considerations by the following 
expression: 

. 1 .21 (eV) 

n Q p Q = 1.5 x 10 33 T 3 e kT (1.1.1) 

= 2.2 x 10 2 ° cm" 6 , for T = 300 K 

_3 

where n Q = the equilibrium concentration of conduction electrons (cm ) 

p Q = the equilibrium concentration of holes (cm ) 

T = temperature (K) 

k = Boltzmann constant (0.8618 x 10’^ eV/K) 

For a highly purified semiconductor, the principal source of charge 
carriers is thermal excitation of electrons from the valence band to the 
conduction band, and the concentration of conduction electrons will equal 
the concentration of holes. This state, in which the electrical properties 
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of a semiconductor are not modified by impurities, may be referred 
to as intrinsic. The electron and hole concentrations in in- 
trinsic silicon, for example, are equal to 1.5 x 10 10 cm at room 
temperature . 

When elements from Column III and V of the periodic table occur in 
substitutional solid solution in silicon, they can be thermally ionized. 

In the case of Column V elements, such as phosphorus or arsenic, the 
ionization results in an electron in the conduction band and a positively- 
charged donor impurity atom in the silicon lattice. Impurities from 
Column III, such as boron, undergo ionization in silicon by accepting a 
thermally-ionized electron from the valence band. This process creates 
a hole in the valence band and a negatively-charged acceptor impurity 
ion. The activation energies for these donor and acceptor atoms in 
silicon are approximately 0.05 eV. For this reason, these equilibrium 
processes go to completion at temperatures near 300 K (kT ^ 0.026 eV), 
and the commonly-used Column III and V impurities in silicon can be 
considered to be completely ionized at room temperature. 

If significant quantities of conduction electrons or holes are pro- 
duced by the addition of impurities, as described above, the semicon- 
ductor may be classed as extrinsic. Extrinsic semiconductors are 
referred to as n-type {i.e., negative type) if the equilibrium concentra- 
tion of conduction electrons exceeds the intrinsic carrier concentration. 
When the equilibrium concentration of holes exceeds the intrinsic car- 
rier concentration of a semiconductor, it is referred to as a p-type 
(ie., positive type). The product of the equil ibrium conduction electron 
and hole concentrations in extrinsic semiconductors remains constant as 
described by equation (1.1.1). Thus, boron-doped, p-type, extrinsic 

silicon with a resistivity of 10 ohm-cm and a hole concentration of 
-3 

1.4 x 10 15 cm must also have a conduction electron concentration 
of 1.6 x 10 5 cm . In this case, the holes are referred to as 
majority carriers and the conduction electrons as minority carriers. 

The concept of Fermi level may also be used to describe several 
aspects of semiconductor theory. The Fermi level of a material is defined 
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as that electron energy state at which the probability of occupancy i:, 
equal to 1/2. The Fermi level Is at the center of the forbidden band 
when silicon is intrinsic. In an n-t.ype semiconductor, the Fermi level 
is above the center of the forbidden band. In a p-type semiconductor 
the Fermi level is below the center of the forbidden band. 


Concentrations of conduction electrons and holes in excess of 
thermal equilibrium values can be Introduced in a semiconductor by 
electrical processes, by the absorption of electromagnetic radiation, 
or in the process of stopping high energy particulate radiation. The 
total instantaneous concentration of carriers during an excitation 
process can be expressed as follows: 

P(t) = p Q + P'(t) (1.1.2) 

n(t) = n Q + n'(t) (1.1.3) 

where p 1 ( t ) and n'(t) are the instantaneous excess hole and electron 
concentrations, which in the general case will be functions of time. 

The absorption in a sample of silicon of electromagnetic radiation, 
referred to as the optical injection of carriers, is fundamental to the 
operation of the solar cell. In the absorption process, an electron- 
hole pair is created for each photon of light absorbed. The densities 
of excess electrons and holes created in this manner obey the following 
equations: 

(1.1.4) 

(1.1.5) 

where g gxt represents the excitation rate per unit volume due to an 
external cause, g^ is the thermal generation rate, and r is the total 
recombination rate. If a net rate of recombination, u, is defined, 

u = r - 9 t h ’ 0.1.6) 

then for the case of holes, for example. 


at** 1- 9 ext* 9th • r ' 
a ? 1 * 1 = w + 9th - r - 
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u. 


(1.1.7) 


dp(t) = 

dt 


9 ext 


It has been found for semiconductors, for the case of small excess 
carrier densities, that is, p(t)« p Q , that a good approximation for u is, 


u ■ (PrfPo’/y 


(1 . 1 . 8 ) 


where p n denotes the case of excess hole density in n-type, or electron- 
conducting material and x p is the lifetime of a hole. The implication 
of this can be seen if the above expression for the time derivative of 
p(t) is integrated, for the case of g gxt = o, with the initial condition 
The result is, 

-t/T„ 


P n (o) - p«> 


p n (t) = p n (o)e 


t > 0, 


(1.1.9) 


and the lifetime is now seen to be the decay time constant governing the 
return of excess holes in n-material if the external source is removed 
at t = 0 . 


An explicit expression for the lifetime, x p , has been developed by 
Hall and Shockley and Read; it is given by the expression, for holes, 

T P * (a P v th N t )_1 ’ (1.1.10) 

where cr p is the cross-section for capture of a hole by what Shockley and 
Read have termed a recombination center, is the thermal velocity of 
an excess carrier and is about 10 7 cm/sec, and N t is the density of the 
recombination centers. These centers, it has been determined, are res- 
ponsible for the recombination of excess carriers, whether injected 
electrically or by electromagnetic, or particulate, radiation. Modifi- 
cation and the creation of additional centers of this type can result 
from the action of high energy radiation in producing lattice displace- 
ments and vacancies, as will be discussed in more detail below. 

1.1.2 Carrier Transport . Current flow or charge transport can occur by 
either of two mechanisms in semiconductors. The drift of charged 
carriers in an electric field is observed in semiconductors as well as 
metals. Such a drift current for the case of holes in a p-type 
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semiconductor can be described as follows: 


Jp = q P Up E. (1.1.11) 

where Jp = hole current density (amperes/cm 2 ) 

q = hole charge (coulomb) 
p = hole concentration (cm ) 

Pp = hole mobility (cm 2 /volt sec) 

E = electric field (volts/cm) 

The coefficients of the electric field (E) in the above expression are 
related to the resistivity of the material in the following manner: 

P (ohm-cm) = — . (1.1.12) 

H P Mg 

Similar expressions can be written to describe conduction electron flow 
and combined expressions can be used if second carrier conduction is 
significant. 


The second mechanism for charge transport in semiconductors is 
carrier diffusion. This process results from the random thermal move- 
ment of particles which exist in a concentration gradient. Such dif- 
fusion is analogous to flow of heat due to thermal gradients and the 
diffusion of atoms and molecules. When a gradient in the concentration 
of holes exists in a semiconductor, a hole flux will flow opposite 
to the gradient. The hole current, for a one-dimensional geometry, is 
shown in the following expression: 


J 


P 




(1.1.13) 


where Jp = hole current density (A/cm 2 ) 

Dp = hole diffusion constant (cm 2 /sec) 

^ = gradient of hole concentration . 
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When both mechanisms contribute to the hole flow, the following equation 
describes the current: 


°p - 1<P % E - D p^> • (1,1.14) 

A similar expression can be written for the conduction electron current 
as follows: 

J „ - 1 ^ E + D n&> • n.i.is) 

In some situations, an expression summing the hole and conduction elec- 
tron current components may be necessary to describe the current. 


The basic equation governing the behavior of charge carriers in 
time and space is the time-dependent continuity equation. This equation 
sums the effects of the divergence of current, carrier recombination, 
and carrier generation. For the non-equilibrium steady state case, the 
total carrier concentrations (n and p) remain constant, and and 
gjjr equal zero. In this case the one-dimensional continuity equations 
for electrons and holes are as follows: 


n - n 


- 1 4- J = 0 

q dx n 


P - P. 


11, - o 

q dx J p °» 


(1.1.17) 


where 9 ext » introduced earlier, is the rate of generation of carriers per 
unit volume. If current flow occurs only by diffusion, equation (1.1.13) 
can be substituted into equation (1.1.17), and a similar substitution 
can be made in equation (1.1.16), leading to the following equations: 


n - n. 


= " 9, 


(1.1.18) 


i He. . p " p q _ 
p dx 2 T 


9 ext 


(1.1.19) 
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1.2 The P-N Junction 


The current-voltage characteristic of a p-n junction is one of the 
factors which determine solar cell response. In this section, the 
general factors which determine diode characteristics will be dis- 
cussed. The carrier concentrations found in a solar cell diode are 
shown graphically in Figure 1.1. The base or p-type region of the device 
has a majority carrier density (p ) of approximately 10 15 cm. 

PO V C 

Because the product of the two carrier concentrations is roughly 10 20 cm , 

_ O 

the minority carrier concentration is 10 s cm" . The surface or diffused 

-3 

layer has a majority carrier concentration approaching 10 20 cm . 
Equilibrium considerations therefore require the minority carrier con- 
centration of this region to be approximately unity. The concentrations 
of holes and conduction electrons differ greatly on different sides of 
the junction. This results in two effects. The minority carriers on 
either side of the junction tend to diffuse across the junction (i.e., 
create a diffusion current). In addition, the very large concentrations 
of conduction electrons existing in the n-type layer form an electro- 
static potential barrier to oppose electron flow from this p-type region. 
As a result, all of the mobile charge carriers (holes and conduction 
electrons) are swept out of the interface region. The ionized impurity 
atoms form a dipole layer in the interface region. This layer is also 
referred to as a space charge layer. In a typical solar cell, the width 
of the n-type diffused layer is roughly 0.4 pm, and the width of the 
space charge region is very roughly 0.5 ym. 

Equation (1.1.17) can be used to determine the behavior of excess 
carriers in the region of a junction. In the case of steady-state 
illumination, 

d 2 n„ - ri 

D £ - P - ■ B P . =-g . x > 0 . 0.2.1) 

" d * 2 

The solution of this equation for a semi -infinite semiconductor with the 
boundary condition, that at x = 0, np = np Q , is 
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The quantity /D t has the dimensions of length and is often referred to 
as the diffusion length (L). The above result indicates that at equilib- 
rium the concentration of conduction electrons in the p-type region will 
approach zero at the junction and will increase exponentially with 
distance from the space charge region. This behavior is shown in Figure 
1.1. At a distance of one diffusion length from the junction, the minority 
carrier concentration is 1/e that of the equilibrium bulk value. Actual 
diffusion lengths found in solar cells can be as large as 200 ym. This 
parameter is of primary importance in the determination of the efficiency 
of a solar cell. 

The equation for the current flow as a function of bias is as follows 
for a p-n junction: 


J = J Q1 (e qV/kT -1) (1.2.3) 

In the case of a large forward bias (V»0), e qV,/k ^ is much larger than 
1 and therefore, 

J = J Q1 e qV/kT (1 .2.4) 

When V< 0, J = -J Q1 . For this reason, J Q1 is also known as the saturation 
current. If the saturation current is assumed to be due to the diffusion 
of minority carriers into the junction, then: 

q D n IT 

Joi ° L n - J °°qn \U 0-2-5) 

n » n 

Based on diffusion limited current, the calculated saturation current 
for an n-p 10 ohm-cm solar cell would be roughly 10 ^ a/cm 2 at room 
temperature. The measured values of saturation currents found in such 
solar cells are considerably higher than the above value. The 
diffusion theory thus does not adequately explain the current voltage 
characteristics of a silicon junction diode. 
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A second theory of the diode current voltage relationship involves 
carrier generation and recombination through defect centers located in 
the space charge region. The diode or rectifier equation predicted by 
this theory is as follows: 1 •"* 

j = J 02 ( e qV/2kT -1) (1.2.6) 


The only difference between equations (1.2.3) and (1.2.6) is the factor 
of 1/2 which appears in the exponent and the form of J 0 2 * The expression 
for J 02 is: 


q W n i 


(1.2.7) 


where W = width of space charge region 


t q = carrier lifetime in space charge region 

n i = intrinsic carrier concentration 1.5 x 10 10 cm ) 


Experimental studies have shown that the generation-recombination model 
and the diffusion model are necessary to describe the diode current flow 
at all voltages. An expression surraning the currents of both models can 
be used to describe the current flow at all voltages. 1,1 

As a result of manufacturing variations, a solar cell junction is 
occasionally shunted by an ohmic resistance, When the value of this 
shunt resistance is less than 10 4 ohms, the shunt current will dominate 
the diode current at forward biases of slightly less than 0.2v. The 
symbol for shunt resistance is R As a result of resistive volume 
elements in current paths to the diode junction, the solar cell 
also has a finite resistance which appears in series with the diode. 

This series resistance (R $ ) is usually less than one ohm and will domi- 
nate the current flow through the diode at large forward biases. A 
model summing both of the above elements is necessary to describe the 
forward voltage-current characteristic of a silicon solar cell in the 
most general case. Such a model is shown in Figure 1.2. In Figure 1.3, 
a generalized current-voltage characteristic is shown for a solar cell 
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Figure 1.3 Typical Dark Solar Cell Current-Voltage 
Characteristic, Forward Bias 
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diode using the above model. Actual solar cells will have considerable 
variation in the shunt and series resistances. 


The junction space charge region of a solar cell has an associated 
capacitance. The capacitance of a conventional solar cell is related to 
the width of the space charge region in the following manner: 

( 1 . 2 . 8 ) 


where C is the capacitance, A is area, e is the dielectric constant 
and W is the width of the space charge region. The acceptor density in 
the p-type region adjacent to the space charge region can be related to 
the capacitance per unit area by: 


N 


a 


2< v a * V c2 

qe A 2 


(1.2.9) 


where is the applied voltage, and is the barrier voltage (0.6 to 
0.8v depending on resistivity of cell base). The above expression assumes 
an abrupt or step junction which is typical of conventional solar cells. 


1.3 Silicon Solar Cell Theory 

When a silicon p-n junction diode is exposed to ionizing radiation 
or light with a photon energy equal to or greater than the band gap of 
silicon, electron-hole pairs are produced in the silicon. Because of 
the gradient of conduction electrons (see Figure 1.1) which exists in 
the p-type region near the space charge region, the conduction electrons 
generated by the radiation diffuse to the junction. When these electrons 
reach the space charge region, they drift in the space charge region 
field to the opposite side of the junction. A similar behavior occurs 
for holes generated in the n-type regions of a solar cell. The dif- 
fusion flux of these generated carriers to the junction is the solar 
cell generated current. Several investigators have developed general 


1-13 



■j 2-1 fl 

expressions for generation current. * ' These expressions are 

solutions of the continuity equations (1.1.18) and (1.1.19) for the case 
of optical carrier generation. The expression for electrons is as 
f ol 1 ows : 


d 2 n 

dx 2 


D n<" 


4i 


+ aN 0 (l - R) e~° d = 0 


(1.3.1) 


where 


1 1 a 

a = absorption coefficient for light of wavelength A, (cm ) 

N q = photon flux density 

R = reflection loss 
x = distance from the junction 
d = distance from front surface 


This equation can be solved to find the minority carrier concentration 
gradient at the edge of the space charge region. The current of car- 
riers into the space charge region can be calculated by evaluating 
current at the edge of the space charge region by use of equation (1.1.15) 
Separate evaluations must be made for the diffused or surface layer and 
electron currents in the bulk response to monochromatic light as 
follows: 1 * 8 

Surface Layer: 


J p U) 


q N n 0 - R) air 


1 - 


(“ l P 4 scj) 5, " h q 4 (’ 4 T 5 ) " 5h e '“ 8 ' (’ 4 V) 

• ■ 


sinh p— + cosh r— 
L P iL P L P 


(1.3.3) 
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Bulk Response: (assuming S = °° at d = b) 


J n (X) 


q N Q (1 - R) aL 


(.< 

n V 


nh 3 - aL cosh n 


T-) 


* °*Tl « 


•ab 


a 2 L 2 
n 


cosh 


b - a 


(1.3.4) 


where a = junction depth (cm) 

b = cel 1 thickness (cm) 

S = surface recombination velocity (cm/sec) . 


Total Response: 


J L (X) = J n (X) + J p (X) (1.3.5) 

The above equations are written for the case of an n-p solar cell. 

It is also assumed that no significant drift fields are present. The 
cell response in A/cm 2 may be normalized to the photon flux density (N Q ). 
In this way, the above equations describe the response of the cell in 
terms of amperes per photon/sec of incident light of a given wave length. 
Solar cell spectral response curves are routinely measured. In these 

experimental measurements, the response is usually normalized to the inci- 

_2 

dent optical power density (watts cm ) rather than photon density rate. 
The calculated response of a typical solar cell in such terms is shown 
in Figure 1.4. 

The previous equations illustrated the role of the minority carrier 
diffusion length in development of the light-generated current of a solar 
cell. These response equations can be folded with the solar spectral 
irradiance and integrated to yield the light-generated solar-cell current 
under solar illumination (see Figure 1.4). 
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Figure 1.4 Calculated Silicon Solar Cell Spectral Response 


The light generated current can be combined with previously dis- 
cussed diode rectifier equations to determine the current-voltage 
characteristic of an illuminated solar cell. The model for an illumi- 
nated solar cell is the same as that shown in Figure 1.2 for a dark 
diode, with the addition of a current source.^ The current source 
(shown dotted in Figure 1.2) represents the light generated current. 

On the basis of the above model, an equation can be written to describe 
the cell current into an external load: 



Several observations can be made regarding the form of the above equation. 
The light generated current is independent of applied voltage and propor- 
tional to the intensity of the incident illumination. The development 
of the light generated current produces a forward bias on the solar cell 
diodes (D1 and D2). The light generated current (1^) will divide between 
the parallel branches containing D1 , D2, R ^ and R $ + R^. The behavior 
of the illuminated solar cell current (I) and voltage (V) as R^ varies 
from zero to infinity is referred to as the I-V characteristic. This 
characteristic is the primary engineering tool used in evaluating solar 
cells. A general expression for the cell current to an external load 
can be obtained by substitution of equations (1.2.4) and (1.2.6) into 
equation (1.3.6). In the case of a good cell under 135 mW/cm 2 solar 
illumination, the current in R $h can be neglected. It has been the prac- 
tice to simplify the two diode currents with the following expression: 



‘ q(V+IR )/nkT 
e s -1 


(1.3.7) 


1-17 



where I is an apparent saturation current, and n is a constant, between 
1 and 2 in value. The resulting expression is often used to describe 
solar cell I-V characteristics : 

I = I L - I Q [e q ( V+IR s)/ nkT -i (1.3.8) 

The development of a solar cell I-V characteristic from the light 
generated current and dark diode characteristic is shown graphically in 
Figure 1.5. An 1^ value of 35 mA/cm 2 is typical of solar cells under 
solar illumination of 135 mW/cm 2 . This 1^ value is shown in Figure 1.5. 

In addition, the dark diode or rectifier characteristics shown in Figure 
1.3 are replotted in Figure 1.5. The diode characteristics are shown 
with and without the series resistance. The illuminated solar cell I-V 
characteristic for a hypothetical cell with R g = 0 is obtained by sub- 
tracting the forward current flowing in D1 , D2, and R ^ from the light 
generated current I L . When R $ is some significant quantity, the dark 
diode characteristic is displaced an amount AV before subtraction from 
I L . The quantity aV is the voltage drop across R g when the solar cell 
diode conducts a forward current equal to +1^. 

It should be understood that this analysis is for a solar cell at 
27°C under solar illumination of 135 mW/cm 2 . The quantity 1^ is propor- 
tional to the light intensity, a function of temperature, and also a 
function of the spectral content of the illumination. The dark diode 
currents 1^ and 1^ are strong functions of temperature. Under the 
assumed conditions of temperature and illumination, 1^ and R $ dominate 
the I-V characteristic of the solar cell. Under other conditions of 
temperature and illumination, the solar cell I-V characteristic may be 
influenced by other factors such as R $ ^ and I^. 

A different set of parameters is used to describe the solar cell 
characteristic for engineering purposes. These are (a) short circuit cur- 
rent I c -(b) open circuit voltage V,. and (c) maximum power P mav . 

The short circuit current is that current produced by the cell when the 
load resistance (R^) approaches zero. In good solar cells, this 
quantity is equal to the light generated current 1^ or A. In cells 
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with high or excessive internal series resistance, or in good cells at 
higher illumination intensities, I will be less than the light 
generated current. The open circuit voltage is the voltage produced 
by the cell when is infinite. In this load condition, all of the 
light generated current is consumed in forward conduction of diodes 
D1 and D2. 





(1.3.9) 


A maximum in the power delivered to the load resistance occurs at some 

point of the solar cell I-V characteristic. The power developed under 

such a load is called the maximum power (P ). A method of determining 

an analytical expression for the I-V characteristic from parameters such 

as I , V . , and V has been described in the literature. 1 -1 1 

sc’ oc max max 

1 .4 Solar Cell Coatings and Contacts 

A silicon solar cell is a composite of several layers of material. 
The layers of n and p-type silicon form the basic cell structure in which 
the current is generated. Additional practical problems are involved in 
maximizing the light which enters the silicon and providing a low 
resistance path for the removal of the generated current from the solar 
cell. When the light passes from one medium to another medium which has 
a different index of fraction, some light is reflected. The amount of 
light reflected can be determined from the following relationship: 



where R = reflectivity (fraction of normal incidence light intensity 
reflected) 

n! = index of refraction, medium 1 
n 2 = index of refraction, medium 2 


The above relationship holds only for normal incidence light. The more 

general case of light incident at an arbitrary angle 0 from the normal 

1 12 

is covered in a basic text on optics. 
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The index of refraction for silicon in the solar cell response region 
varies between 3.5 and 4.0 in the optical region. For this reason, there 
are significant reflectance losses when light is incident on an air-silicon 
interface. The reflectance loss can be reduced by adding an intermediate 
layer. The optimum index of refraction for the intermediate layer is given 
by: 


n op = (n i • n 2 ) 1/2 (1.4.2) 

The optimum index of refraction for an anti reflection coating at an air- 
silicon interface is approximately 1.9 (Silicon monoxide, with an index 
of refraction in the range of 1.8 to 1.9, is used as an anti reflection 
coating on solar cells.) Silicon solar cells used in the space environ- 
ment are always covered by glass to shield against radiation and to 
raise the effective emissivity for better thermal control. Recently, 
efforts have been made to develop an anti reflective coating optimizing 
the transmission from glass adhesive to the silicon. Since the adhesives 
have n values of approximately 1.4, the previous equation indicates that 
an anti reflective coating with n = 2.2 - 2.4 would be optimum for a solar 
cell to be used with cover glass. Titanium oxide and tantalum oxide have 
been proposed and evaluated for this purpose.^' Anti reflecting 

coatings may vary in thickness from 70 to 80 nm. In theory, the coating 
thickness should equal 1/4 the light wavelength for minimum reflection 
loss. 


The contacts of current commercial solar cells are formed by evapora- 
ting titanium and silver metal on the entire back surface and front con- 
tact pattern. The total thickness of this evaporated metallization is 
approximately 5 ym. After the metallization, the cells are usually solder 
dipped. The solder thickness may vary between 0.001 and 0.008 cm 
(0.4 - 3 mils). One of the primary considerations in the selection of the 
contact is the electrical behavior of the metal -semiconductor interface. 
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In general, such interfaces can be described as Schottky barriers. A 
Schottky barrier has a current voltage characteristic of the same form as 
that for a p-n junction. The saturation current for a Schottky barrier 
is as fol lows : ^ 


I 0 = A T 2 exp - U B /kT) (1.4.3) 

-2 -2 

where A = effective Richardson constant (A cm K ) 

= barrier height ( eV ) 

-2 -2 

The quantity A is approximately 100 A cm K and <f>g is approximately 

0.50 (eV) for most metals in contact with p-type silicon. The satura- 

_2 

tion current (I ) at room temperature (T = 300 K) will be between 10 
and 10 A cm . The effect that the Schottky barrier has on the solar 
cell will be related to the forward resistance of the barrier. Since 
the form of the barrier current-voltage characteristic is: 

I = I 0 (e qV/kT -l) (1-4.4) 

The dynamic impedance of the junction is as follows: 

dv . kTZa e -qV/kT (1. 4 .S) 

dI 'o 


It can be seen that impedance of this barrier is inversely proportional 
to the saturation current. Since the saturation current at room tempera- 
ture is very high, the impedance of a Schottky barrier is very low. If 
the barrier potential (<j>g) for a particular metal on silicon is low enough, 
the barrier I-V characteristic will approach low resistance ohmic behavior. 
This is the case for a titanium layer on p-type silicon at room temperature. 
At low temperatures the saturation current of such a Schottky barrier is 
reduced and the diode characteristics become more significant. In this 
case, the Schottky barrier adds a nonlinear voltage drop to the solar cell 
model in series with R $ This problem has received considerable attention 
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in the literature. The problem associated with nonohmic 

contacts can be reduced by producing a heavily doped (p+) layer on the 
silicon interface. In such cases, the space charge region associated with 
the Schottky barrier is generally reduced. Quantum mechanical tunneling 
of the space charge region dominates the behavior of such thin barriers, 
and provides a highly conductive metal semiconductor interface. Since 
the solar cell front contact is applied to a silicon interface which is 
very heavily doped due to phosphorous diffusion, the tunneling mechanism 
assures a low resistance ohmic contact. 
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CHAPTER 2 


2.0 INSTRUMENTATION TECHNIQUES FOR MEASUREMENT OF SOLAR CELL PARAMETERS 

In this section, the various commonly used experimental methods for 
the analysis of radiation effects are discussed. The most commonly used 
measurement in the analysis of radiation effects in silicon solar cells 
is the current-voltage characteristic under illumination. Since solar 
cell response is a strong function of optical wavelength, the light 
source is a major variable in the evaluation of solar cell parameter 
changes. 

2.1 Light Sources and Solar Simulators 

The spectral irradiance of the sun at 1.5X10 11 m (one AU) is 

of primary importance in solar cell analysis for earth orbits. The 

values of solar spectral irradiance proposed by Johnson in 1954 have 

been widely accepted and used for analytical purposes. * Johnson's 

results indicated that the solar constant was 139.5 mW/cm 2 , and also 

that the solar spectrum closely approximates that of a 6000 K 

blackbody. Several high-altitude measurements made in recent years 

2 2 

have been reviewed. ’ The findings indicate a solar constant of 
135.3 ± 2.1 mW/cm 2 . The proposed solar spectral irradiance is tabu- 
lated in Table 2.1. Silicon solar cell response is generally limited 
to the region between 0.3 and 1.2- pm. In this range, the solar power 
density is 104.4 mW/cm 2 . 

Among several solar simulation techniques, the most common method 
is the use of a xenon arc lamp with filters to remove undesired line 
spectra in the near infrared. Unfiltered xenon lamps are also used in 
the pulsed mode, which does not generate the undesired line spectra. 
Unfiltered carbon arcs are also used to simulate solar illumination with 
a reasonable spectral match. A close spectral match to the solar spectrum 
is obtained by the use o,f a xenon-filtered tungsten combination or filtered 
xenon source. These sources match the solar spectrum well enough that 
cell measurements made under these sources can be considered representative 
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TABLE 2.1 

SOLAR SPECTRAL IRRADIANCE - PROPOSED STANDARD CURVE 


- wavelength in micrometers 

E(-) - Solar spectral irradianrc averaged over small bandwidth centered at 1, in W,m \ uni * 
£ jo-'* 1 - Area under the solar spectral irradiance curve in the wavelength range 0 to \ in 

D<o-0 - Percentage of the solar constant associated with wavelengths shorter than X 
Solar Constant - 1 J S 3 W»m 2 
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30 

76. 

,519 
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97.5030 
97.0277 
90.0303 
90.2179 
90. 3724 

90.5047 
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of that in the space environment. Many other types of light sources have 
been used in radiation effects studies. Unfiltered and filtered incan- 
descent tungsten Sources have peak response in the red and near infrared. 
Since this is the wavelength region of the solar cell response which is 
most changed by irradiation, use of these sources will show much more 
severe cell radiation degradations as compared to evaluations with a 
suitable solar simulator. This characteristic severely limits the use 
of tungsten source simulation in the evaluation of radiation effects. 

Filtered xenon arc solar simulators are manufactured by Spectrolab 

and Aerospace Controls. The spectral irradiance for the Spectrolab X-25 

2 3 

simulator is shown in Figure 2.1. * Similar data is shown in Figure 

2.2 for a combination xenon-tungsten source simulator used by Centralab 

Semiconductor Products. In Figure 2.3, the spectral irradiance is 

. 2 5 

shown for a Genarco Model #ME4CWM carbon arc simulator. 

An important recent development in the field of solar simulators 
is the use of pulsed xenon arc lamps for solar cell and solar cell array 
testing. ’ These developments have been prompted by the need for 
a suitable alternative to testing large arrays in natural sunlight on 

the earth's surface. In these systems, an approximately 2 msec pulse 
of light is produced. Solar cell output data can be accumulated during 
about 1 msec of the pulse length. Sophisticated electronic data handling 
systems are necessary to record cell or array outputs and commutate 
external load resistances during the light pulse. Variations in test 
cell current due to light intensity variations are corrected to a normal- 
ized value at the desired illumination. The high intensity peaks in the 
0.8-0. 9 ym region of the xenon arc spectrum are not generated by pulsed 

operation with high current densities. By this means, it is possible 

2 7 PR 

to achieve a reasonably close match to the solar spectrum. * 

Array areas up to eight feet by eight feet can be illuminated with 
excellent temperature control, by such simulators. 
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Figure 2.2 Spectral Output, Centralab Solar Simulator 
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The extent to which a lack of solar spectral match affects a solar 
cell measurement can be estimated if the spectral intensity of the light 
source and the spectral response of the solar cell are known. The light 
generated current of the illuminated cell can be calculated as follows: 


I L (A/cm 2 ) = / R(A).E(A) dA 


( 2 . 1 . 1 ) 


where R(a) = solar cell response as a function of wavelength, A/W 
E(a) = spectral irradiance, W/cm 2 -ym 
dA = an increment of wavelength, ym 


The above equation can be used to determine the light-generated currents 
under solar and simulator illuminations. The generated current under solar 
illumination can be calculated from the generated current under simulator 
illumination if the spectral response of the cell if the spectral irradiance 
of the simulator is known. The relation is as follows: 


/ R < x >' E < X >space dX 

I (space) = I (simulator) (2.1.2) 

/ R(A)-E{X) s1m . dX 

Solar simulator intensities are determined by the short circuit 
current outputs of calibrated primary or secondary standard cells. The 
primary standard cells, commonly in use, were generated by a NASA/JPl 

7 a 7 in 

program of telemetered balloon flights. * Similar programs have 
been conducted by aircraft and high altitude terrestial measurements. 2.11-2.1 4 
When the effects of atmospheric absorption are properly corrected, the 
results are in good agreement with the balloon flight data.^'^ 
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Primary standard cell availability is limited and they are considered 

too valuable for general usage in setting simulator intensities. For this 

reason, secondary standard cells are calibrated for use as working standards. 

Palmer has recently reviewed the methods of generating secondary standard 

cells and concluded that previously proposed methods of calibration may 

? 16 

yield poor results. ' Palmer has proposed the use of alternate methods 

which insure that secondary standard cell calibration accuracy will approach 

2 16 

that of primary standard cells. A solar simulator intensity which pro- 

duces a standard cell response equal to that for free space at one AU 
is referred to as one sun, air mass zero or 135.3 mW/cm 2 (formerly 
139.5 mW/cm 2 ). Solar simulator spectral quality should be monitored 
by use of narrow bandpass or cutoff filters with calibrated spectral 
response detectors. 


2.2 Current-Voltage Characteristics 

The measurement of solar cell current-voltage characteristics is 
the primary means of evaluating the device. The evaluation is made by 
measuring the cell voltage developed and the cell current into load 
resistances varying between zero and infinity. The measurement is 
simple in principle but attention to several practical details is 
necessary to insure accurate results. Solar cell response is a strong 
function of temperature. For this reason, the cell must be in thermal 
equilibrium at a known temperature during the measurement. With ade- 
quate heat sinking and cooling, cells measured under one sun irradiance 
at room temperature can be stabilized at 28°C. To insure that the 
voltages measured are representative of those developed on the cell 
contacts, separate probes are employed to measure cell voltage 
and current. In this way, any voltage drops which occur at the current 
probe-cell interfaces due to contact resistance do not cause errors in 
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the measured cell voltage. The load resistance may be varied manually 

or electronically. The current voltage data is usually plotted with an 

X-Y recorder. The solar cell parameters such as I and V can be read 

sc oc 

directly with digital meters. Multiplier circuits are available which 

produce a voltage proportional to the product of cell voltage and current. 

This output is plotted as a function of cell voltage to directly 

indicate the maximum power and voltage at maximum power. The cell 

series resistance is also determined from current-voltage charac- 

217218 

teri sties at two or more different illumination levels. * ’ * 

2.3 Spectral Response Measurements 

Spectral response measurements are very useful for evaluating 
changes in solar cells due to radiation effects. The spectral response 
(A/W) is a measure of the short-circuit current density generated by 
the cell under various monochromatic illuminations of a known power 
density. The spectral response is often reported in terms of relative 
units when absolute values of light intensities are not determined. 

Various schemes have been used to measure the spectral response of 
solar cells. High resolution monochromators are used when extreme 
accuracy is desired. When less accuracy is needed, narrow bandpass 
filters can be used as sources of monochromatic light. When a mono- 
chrome tor is used, there are two methods to normalize the solar cell 
output to the light intensity. Tungsten light sources are usually used 
in monochromators, and the entrance slit width can be varied to control 
the optical power density illuminating the cell under test. In some 
systems, the entrance slit width can be automatically controlled to 
maintain a constant optical power density on the solar cell. An alter- 
nate approach is to maintain a constant slit width and allow the optical 
power density on the cell to vary with wavelength, and attenuate the 
cell response at each wavelength measured by use of calibrated voltage 
dividers. 

One disadvantage of these methods of measurement is that the solar 
cell response is determined at very low minority carrier injection levels. 
Solar cells irradiated with neutrons and protons have response 
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characteristics which are dependent upon the concentration of injected 

minority carriers. In such cases the cell must be illuminated with a 

light source similar in intensity and spectral content to the intended 

space environment during the spectral response evaluation. A scheme 

for measuring spectral response under approximate solar illumination 

2 19 

has been suggested by Stofel . ' 

2.4 Irradiation Methods 

The evaluation of solar cell radiation effects requires a wide 
range of specialized equipment and instrumentation. Charged particle 
accelerators are the primary sources for space radiation simulation. 

The range of electron energies of interest is 0.3 to 10 MeV. Electron 
energies of 0.3 to 3 MeV are usually obtained with Van de Graaff elec- 
trostatic or Dynamatron accelerators. Higher electron energies can be 
reached with linear electron accelerators. Proton energies from 0.1 to 
3 MeV are obtained with Van de Graaff accelerators. Proton energies 
greater than 10 MeV can be obtained from cyclotrons. For lower energy 
protons, it is necessary to transport the proton beam and perform the 
irradiation in vacuum to avoid excessive energy losses. A survey of 

2 20 

accelerator facilities has been published but is currently out of date. 

The Space Radiation Effects Laboratory at Newport News, Virginia is oper- 

2 21 

ated for NASA specifically for space radiation simulation. * The facili- 
ties include accelerators for all electron energies of interest and a 
600 MeV synchrocyclotron. Accelerators invariably produce irradiation 
rates which are many orders of magnitude greater than those of space 

environments. Real time irradiations of solar cells have been done 

222 223 

using beta emitting sources. ’ ’ ’ These sources generate a spec- 

trum of electron energies and fluxes similar to that of some space 
environments. 

A successful experiment must include accurate knowledge of the 
particle energy, measurement of cross sectional beam intensity at the 
irradiation area, as well as the intensity during the irradiation. 

Although there are several methods of accomplishing the above measure- 
ments, all can be done with a Faraday cup. A design of a Faraday cup 
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p p A 

suitable for accelerators in the 1 MeV range is shown in the literature. 

The desirable characteristics of a Faraday cup are as follows: 

a. Shielding thickness must exceed particle energy range. 

b. A high cup length-to-diameter ratio is desirable. 

c. Use low Z (atomic number) materials to reduce secondary 
electron emission and bremsstrahlung production. 

d. Cup must be in vacuum or potted. 

e. Cup should be a reentrant cavity. 

f. Cup should be screened to suppress secondary electron 
emission if necessary. 

g. Cup should have remote X-Y translation to facilitate beam 
mapping. 

A Faraday cup requires a current measuring instrument which operates 
-10 “6 

in the range of 10 to 10 amperes and integrates charge. Instruments 
of this nature are produced by Keithley Instruments, Cleveland, Ohio 
and El cor Products, Silver Spring, Maryland. 

The particle energy can be determined by means of a range-energy 
measurement. In this measurement, increasing thickness of absorbers 
are introduced into a constant flux beam, while the flux of particles 
exiting the absorber is monitored with a Faraday cup or a radiation- 
degraded solar cell. If the beam is monoenergetic, a plot of cup current 
(or cell photocurrent) versus absorber thickness is extrapolated to zero 
current at an absorber thickness which is equal to the projected range 
of the mean particle energy of the beam. This technique is satisfactory 
for electrons and high energy protons. Since the beam current must 
remain constant as absorber thickness is increased, a second independent 
means of monitoring beam current must be available. Van de Graaff 
generators are equipped with a generating voltmeter which produces a 
dc voltage proportional to the potential difference on the accelerator 
tube. A check calibration at one operating energy is sufficient to 
insure accurate calibration. Corrections must be made for energy loss 
in the accelerator exit window and in the atmosphere between the exit 
port and the target. 
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2.5 Diffusion Length Measurement 

The importance of minority carrier diffusion length (or lifetime) 
in the study of solar cells was discussed in Section 1.0. The changes 
in this parameter are of primary importance in the solar cell degrada- 
tion in radiation environments. A method using gamma ray irradiation for 
experimental measurement of minority carrier diffusion length was suggested 

p pC 

by Gremmelmaier. This technique requires the uniform generation of 
electron-hole pairs throughout the active volume of the p-n junction de- 
vice. Under these conditions the generated current density (J^) is 
expressed as follows: 


J L = e q o (L p + W + L n } (2.5.1) 

where e = electronic charge 

q Q = generation rate of electron-hole pairs 
Ip = hole diffusion length in an n-type layer 
W - width of space charge layer 
l_ n = electron diffusion length in a p-type layer 


Since l and W are usually very small compared to L n , they may be 
neglected, and the measured short circuit current becomes proportional 
to the diffus-ion length in the p-type base region (L n ). The generation 
rate, determined for this uniform radiation> thus allows accurate deter- 
mination of diffusion length from the measured short circuit current. 


There are several experimental methods of uniformly injecting electron- 
hole pairs. In addition to the use of gamma radiation, high energy elec- 

p pc p 27 p 28 

trons, high energy protons and infrared light ’ * have been used 

to achieve uniform injection. To achieve uniform injection in a solar 
cell irradiated with 1 Mev electrons, it is necessary to introduce a 
0.030 cm (0.012 in.) aluminum shield immediately in front of the cell during 
a normal incidence front irradiation. The details of this procedure and 
the experimental evaluation of the generation rate have been covered by 

p pc 

Rosenzweig. 
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The experimental measurement of diffusion length by the above 
methods has several inherent limitations. Since the diffusion length 
is that distance from which 1/e of injected minority carriers will 
diffuse to the junction during their lifetime, the diffusion length 
concept involves both minority carrier lifetime and diffusion. Minor- 
ity carrier lifetime, in the most general case, could vary throughout 
the active region of a solar cell. In practice this situation arises 
when solar cells are irradiated with low energy protons which do not 
penetrate the entire active volume of the cell. Diffusion lengths of 
solar cells with a nonuniform minority carrier lifetime in the active 
base region cannot be measured by the above methods. Surface recombina- 
tion at the solar cell back contact can also cause errors in measured 
diffusion lengths. These errors are negligible for cells in which the 
thickness exceeds two or three times the diffusion length. 


The measurement of diffusion length by the above methods also assumes 
the external cell-current generated is collected entirely by diffusion of 
excess minority carriers to the junction. Some recent designs of solar 
cell structure utilize "drift fields." In such cases, excess minority 
carrier collection is aided by the presence of an electric field in the 
base region; and the short circuit current under conditions of uniform 
pair production cannot be related to diffusion length by the above equa- 
tion. 


An additional limitation arises if 1 MeV electron or other high 
energy radiations are used in the diffusion length measurement. The 
radiation flux must be kept low to minimize damage to the cell during 
the measurement. The 1 MeV electron beam current during such a measure- 
ment is approximately TO -9 A/cm 2 . The generation rate of excess minority 
carriers produced by this electron flux is considerably lower than that 
produced by solar illumination at 135 mW/cm 2 . In most cases the diffu- 
sion length or lifetime is not dependent upon the concentration of excess 
minority carriers. In such cases, the diffusion length measured with low 
levels of injected minority carriers is the same as that for a cell under 
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one sun illumination. Silicon solar cells irradiated with protons and 

neutrons exhibit injection level dependence of the diffusion length, 

and must be illuminated with simulated solar illumination to allow accu- 

2 29 

rate measurement of the diffusion length. The schemes used by Denney ' 

2 19 

for proton irradiated cells and Stofel , et al. ’ for neutron irradiated 
cells have shown that the diffusion lengths of such cells measured under 
approximate solar illumination are roughly two times greater than that 
measured under low injection level conditions. 

2.6 Statistical and Error Analysis 

The analysis of radiation effects in solar cells involves the col- 
lection and evaluation of large amounts of experimental data. Such data 
are often presented in graphical or tabular form without the use of 
statistical analysis. The basic statistical tools, which have application 
in the analysis of radiation effects data, along with the nature and 
causes of errors in the data, will be discussed. 

A common situation in the analysis of irradiated solar cells involves 
measurements of engineering or physical parameters on a group of solar 
cells exposed to a particular environment or radiation fluence. The 
characteristics of such a data set are a central tendency and a distri- 
bution about this central value. The most suitable measure of the cen- 
tral tendency is the arithmetic mean (x ) . The arithmetic mean is defined 
as follows: 

n 

x = l ^ (2.6.1) 

i=l 

where x. = measured value on cell i 
n = number of cells measured 

The degree to which measured values are dispersed or scattered can 
be described in several ways. One such measure is the range or differ- 
ence between the largest and smallest measured values in a group. A 
second measure of the dispersion in a group of measurements is the 
standard deviation. The standard deviation (a) is defined as follows: 
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( 2 . 6 . 2 ) 



(x i -3T) 2 I 


1 

2 


n - 1 


An additional measure of dispersion is the variance, which is defined as 
the square of the standard deviation (a 2 ). 

Solar cell evaluations usually involve sample sizes of less than 30 
and more often from 3 to 10. When such small sample sizes are used, the 
mean of experimental results may not be representative of a similar experi- 
ment involving a large number of samples. The maximum difference which 
may occur between the mean value of a small sample and the mean value of 
a large sample having a normal or Gaussian distribution can be calculated 
for any desired probability. This difference is referred to as the con- 
fidence limits of a mean value. The confidence limits expressed relative 
to the mean value are as follows: 

x±t -^r- (2.6.3) 

V" 

where t c is the critical percentile value for Students' t distribution. 
Values of t c may be obtained from statistical manuals or handbooks. For 
example, if 95% confidence levels are desired for a sample size of 5, t c 
is equal to 2.78. 

When solar cell experimental data are collected for increasing radia- 
tion fluences, it is desirable to display the data graphically in a manner 
which reveals fundamental empirical relationships. As an example, it 
has been empirically observed that many parameters of irradiated solar 

cells such as I , V , and P can be described by equations similar to 
sc oc max 

the following relation: 

P n»x (t) = P maxo ' C log 0 + * /p x> (2 ' 6 ' 4) 
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where P ($) = maximum power to load at fluence $> 
max r 

P = maxtmum power to load before irradiation 
maxo r 

$ = radiation fluence 

$ - a constant dependent on radiation energy and 

x cell type 

C <= a constant dependent on radiation type 

Computer programs can be employed to fit equations to experimental 
data. In most cases, however, curves can be fitted to experimental data 
by manual means as well. 

The coefficient of correlation is a measure of the degree to which 
an analytical expression can describe variations in experimental data. 
The coefficient of correlation is defined as follows: 


r 


_ + 




n 

I (X, 

i=l 




1 

2 


(2.6.5) 


where x e ^ n = value of x from fitted analytical expression 
x = mean value of experimental data 
x,j = experimentally observed value 

If the fitted analytical expression or curve has a correlation coefficient 
of 1 , it explains all of the experimental ly observed variation, If an 
unexplained variation exists, the correlation coefficient will be between 
-1 and +1. It is usually possible to describe experimental solar cell 
data by expressions which have correlation coefficients very close to one. 
A few examples of the application of statistical methods to solar cell 
data have appeared in the literature.^"^’ ^>31» 2.32 
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In addition to the random errors that affect data, there are syste- 
matic errors which are inherent to the operation of instruments. These 
errors involve limitations of the accuracy with which electron fluence, 
electron energy, solar simulator intensity, etc., can be determined with 
particular instruments. When an experiment is designed, it is desirable 
to insure that the random and systematic instrumentation errors are neg- 
ligible or minor with respect to the accidental errors. In such cases, 
the results can be analyzed assuming that the dispersion or distribution 
in a group of measurements is due only to the accidental error and 
characteristic deviation in the sample group due to manufacturing dif- 
ferences . 
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CHAPTER 3 


3.0 RADIATION EFFECTS 

The behavior of solar cells in a radiation environment can be 
described in terms of the changes in the engineering output parameters 
of the devices. This approach limits the understanding of the physical 
changes which occur in the device. Since other environmental factors 
may need consideration, an understanding of a physical model provides 
a basis for estimates of the behavior in a complex environment. In 
addition, solar arrays of the future will become more complex and may 
utilize materials which are affected by different aspects of radiation 
damage. For these reasons, the engineer should be aware of all aspects 
of the interaction of radiation with matter, and understand the physical 
models which describe the processes. 

3. 1 The Theory of Radiation Damage 

The radiation usually of interest in the study of degradation of 
materials and devices consists of energetic or fast massive particles 
(i.e., electrons, protons, neutrons, or ions). The origin of these par- 
ticles may be particle accelerators, the natural space radiation 
environment, nuclear reactions, or secondary mechanisms such as Compton 
electrons produced by gamma rays. Because they have mass, energy and 
possibly charge, these particles or other particles generated by them 
can interact in several ways with materials. The dominant interactions 
are: 

a. Inelastic Collisions with Atomic Electrons . Inelastic 
collisions with bound atomic electrons are usually the 
predominant mechanism by which an energetic charged par- 
ticle loses kinetic energy in an absorber. In such a 
collisions, electrons experience a transition to an 
execited state (excitation) or to an unbound state 
(ionization) . 

b. Elastic Collisions with Atomic Nuclei . Energetic charged 
particles may have coulombic interactions with the positive 
charge of the atomic nucleus through Rutherford scattering. 

In some cases the amount of energy transferred to the atom 
will displace it from its position in a crystalline lattice. 
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This energetic displaced atom may in turn undergo similar 
collisions with other atoms of the material. Energetic 
particles may also interact directly by a hard sphere colli- 
sion with the nucleus. The probability of this type of 
event is usually less than that for Rutherford scattering, 
except at higher energies. If sufficient energy is trans- 
ferred to displace an atom from its lattice site, that atom 
will probably be energetic enough to displace many other 
atoms . 

c. Inelastic Collisions with Atomic Nuclei. This general 
category of interaction includes several processes which 
are important in radiation damage studies. Highly ener- 
getic protons undergo inelastic collisions with the atomic 
nucleus. In this process the energetic proton interacts 
with the nucleus and leaves the nucleus in an excited or 
activated state. The excited nucleus emits energetic 
nucleons and the recoiling nucleus is displaced from its 
lattice site. This recoiling nucleus in turn causes more 
displacements. This process is also referred to as spal- 
lation. Collisions between neutrons of thermal energy 
and nuclei can also be included in this group. However, 
these interactions are of little importance in solar array 
degradation. 

The major types of radiation damage phenomena in solids which are 
of interest to the solar array designer are ionization and atomic dis- 
placement. It is important to classify an effect into one of these two 
categories, if possible, because the general behavior of each phenomenon 
has been characterized to a large extent. 

Ionization 


Ionization is the process of removing orbital electrons from an atom. 
It may occur in atoms or molecules in gases, liquids or solids. The 
measure of the intensity of an ionizing radiation is the Roentgen. This 
unit is defined by a charge generation of 2.58 x 10' 4 coulomb/kilogram 
of air. The measure of the absorbed dose in any material of interest 
is usually defined in terms of the absorbed energy per unit mass. The 
accepted unit of absorbed dose is the rad (100 ergs/gm or 0.01 joules/kg). 
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Through the use of the concept of absorbed dose, various radiation 
exposures can be reduced to absorbed dose units which will reflect the 
degree of ionization damage in the material of interest. This concept 
can be applied to electron, gamma and X-ray radiations of all energies. 

For electrons, the particle fluence is multiplied by the electronic stop- 
ping power of the electron energy of interest to determine the absorbed 
dose. In this manner, the effects of an exposure to fluxes of trapped 
electrons of various energies in space can be reduced to an absorbed dose. 

In general, this practice is also applicable to proton irradiations; how- 
ever, some caution must be exercised. In some types of materials, the 
effects of the ionization caused by heavy particles is confined to the 
vicinity of the particle track. If homogeneous ionization is produced 
by protons in the absorber material of interest, one can convert proton 
fluences to absorbed doses and sum them with doses from other radiations. 

The variations of stopping power and range for electrons and pro- 
tons of various energies can be seen in Figures 3.1 and 3.2. The 
data presented are for silicon and have been normalized for density in 
the usual manner. The stopping power and range of a fast particle are 
not strong functions of the atomic number of the absorber material. 

For this reason, the data in Figures 3.1 and 3.2 can be used for 
materials with a similar atomic number with a negligible error. 

Radiation may affect solar cell array materials by several ionization- 
related effects. The reduction of transmittance in solar cell cover slides 
is an important effect of ionizing radiation. The darkening is caused 
by the formation of color centers in glass or oxide materials. The color 
centers form when ionizing radiation excites an orbital electron to the 
conduction band. These electrons become trapped by impurity atoms in the 
oxide to form charged defect complexes which can be relatively stable at 
room temperature. Radiation produces similar darkening effects which may 
be observed in cover slide coatings and adhesives. 

Radiation produces many ionization-related effects in organic materials. 
These changes all result from the production of ions, free electrons, and 
free radicals. As a result of these actions, transparent polymers are 
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Figure 3.2 Stopping Power and Range Curves for Protons in Silicon. 
Reference, Janni , AFWL-TR-65-1 50, 1966. 
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darkened and crosslinking between main-chain members may drastically 
alter the mechanical properties. The contemplated use of polymeric 
materials in flexible extended solar arrays will require the array 
designer to have knowledge of the ionization-related radiation effects 
in those materials. 

The use of silicon dioxide as a surface passivation coating and 
dielectric material in silicon devices is the reason for a wide range 
of ionization-related radiation effects in electronic hardware. The 
development of trapped charges in these oxides may cause increased 
leakage currents, decreased gain, and surface channel development in 
bipolar transistors and increased threshold voltages in MOS field 
effect transistors. The presence of ionizing radiation in silicon 
excites the electrons of the valence band to the conduction band, 
creating electron-hole pairs in much the same way that carrier pairs 
are generated by visible light. Although an optical photon of energy 
equal or greater than 1.1 eV will create an electron-hole pair, roughly 
three times this amount of energy must be absorbed from a high energy 
particle to produce the same carriers. In silicon devices, the electron- 
hole pairs which are generated by ionizing radiation cause photocurrents 
in the same manner as solar illumination. 

Atomic Displacement 

The loss of energy by fast electrons and protons caused by collision 
processes with the electrons of an absorber or target material account 
for a large fraction of the dissipated energy. For electrons and protons 
in tne range of 0.1 to 10 MeV, these electron collisions determine the 
particle range in an absorber. Despite this fact, a different type of 
collision process is the basis for the damage which permanently degrades 
silicon solar cells in the space environment. The basis for this damage 
is the displacement of silicon atoms from their lattice sites by fast 
particles in the crystalline absorber. These displaced atoms undergo 
other reactions and finally form stable defects which produce significant 
changes in the equilibrium carrier concentrations and the minority 
carrier lifetime. 
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The displacement of an atom from a lattice site requires a certain 
miminum energy similar to that of other atomic movements. The energy 
of sublimation for a silicon atom is 4.9 eV. The energy for the forma- 
tion of a vacancy in the silicon lattice is 2.3 eV. The displacement 
of an atom involves the formation of a vacancy, the formation of an 
interstitial atom and other electronic and phonon losses. It is reason- 
able to expect that the energy of displacement is several times larger 
than the energy of formation for a vacancy. Seitz has estimated that 

3 1 

the displacement energy is roughly four times the sublimation energy. * 

Electron threshold energies of 145 keV and 125 keV have been reported by 

3 2 3 3 3 4 

various investigators. * * ’ The following equation relates the 

3 1 

electron threshold energy to the displacement energy. ' 

m Q E. 

E d = 2 — — — (E t + 2 m e c2) (3.1.1) 

M m g c 2 ' c 

where E d = displacement energy (MeV) 

E t = threshold energy (MeV) 
m g = electron mass (1/1836) 

M « atomic weight, Si (28) 

m e c 2 = electronic mass-energy equivalence, 0.511 MeV 

The reported threshold energies indicate displacement energies of 12,9 eV 
or 11.0 eV, respectively . 


Although proton threshold energies have not been determined, they 
can be calculated from the classical form of the above equation: 


E 


d 


4 • M • M 

2 e 

(M p + M) 2 4 


(3,1.2) 


where M p = proton mass, 1. The above values of displacement energies 
indicate proton or neutron thresholds of 97.5 or 82.5 eV in silicon. 
Since particles below the threshold energies cannot produce displacement 
damage, the space environment energy spectrums are effectively cut off 
below these values. 
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For particles above the threshold energy, the probability of an 
atomic displacement can be described in terms of a displacement cross 
section. Using this concept, the number of displacements can be esti- 
mated from: 

N d = n a a v $ (3.1.3) 

where = number of displacements per unit volume 

n = number of atoms per unit volume of absorber (5 x 10 22 
silicon atoms/cm 3 ) 

a = displacement cross section (cm 2 ) 

O = average displacements per primary displacement 
$ = radiation fluence (particles/cm 2 ) 

The displacement cross sections for fast electrons of various ener- 
gies can be calculated from the relativistic generalization of the Ruther- 

3 1 

ford scattering cross section equation, ’ For silicon, the calculated 
displacement cross section for 1 MeV electrons is about 68 x 10 cm“ 
and increases only 10% for electron energies of 5 MeV and greater. The 
electron displaced silicon atom may receive enough energy to in turn 
displace other silicon atoms. The mechanism for these secondary displace- 
ments is Rutherford interactions for silicon atoms of energies greater 
than 10 3 eV and hard sphere collisions for lower energy atoms. Although 
different theories of the production of secondary displacements have been 
presented, their results are very similar. Using the model of Kinchin 
and Pease, * the average number of displacements in silicon is 1.53 for 
a 1 MeV electron. The electron energy variation of the various parameters 
is shown below. 
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TABLE 3.1 

SILICON DISPLACEMENT PARAMETERS, VARIOUS ELECTRON ENERGIES 


Electron 

Energy 

(Mev) 

a 

(10~ 24 cm 2 ) 

V 

av 

(10“ 24 cm 2 ) 

n av 

a 

(cm" 1 ) 

AN d /A$ 

5.2 

1 

68 

1.53 

104 

5.2 

1.0 

2 

73 

2.00 

146 

7.3 

1.4 

5 

77 

2.76 

212 

10,6 

2.0 

10 

77 

3,39 

261 

13.0 

2.5 

20 

77 

4.09 

314 

15.7 

3.0 

40 

77 

4.74 

363 

18.2 

3.5 


The direct result of the radiation is the production of vacant lattice 
sites (vacancies) and silicon atoms which come to rest in the interstices 
of the crystal lattice (interstitials). The distribution of vacancies will 
not be uniform, because the vacancies from secondary displacements will lie 
relatively close to the associated primary vacancy. 

The experimental studies must be reviewed to gain a more complete 

model of displacement damage in silicon. Vacancies and interstitials are 

reactive chemical species and are particularly mobile and unstable at room 

temperature. In n-type silicon, it has been shown that vacancies react 

with oxygen impurities to form close coupled vacancy-oxygen pairs 

(V-0) 3 ' 6 " 3 ’ 9 (see Figure 3.3), and with impurity donor atoms, such as 

phosphorus and arsenic, to form close coupled vacancy-donor pairs 

(V-P, V- As) 3 ’ 6, 3,10 (see Figure 3.3). Both defects are electrically 

active and can become negatively charged by accepting an electron from 

the conduction band. The acceptor energy levels of the V- 0 and V-P pairs 

3 11 3 12 

are 0.17 eV and 0.4 eV below the bottom of the conduction band. ’ ’ 

These defects are recombination centers and their formation during elec- 

3 13 

tron irradiation of n-type silicon reduces the minority carrier lifetime. ! 
Since these defects are formed from single vacancies, considerations of 
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A model of the oxygen-vacancy complex or 
the Si-A center. 


A model of the phosphorus-vacancy complex 
or the Si-E center. 




Modei of the divacancy. 


RCA tentative model for the K-center. 


Figure 3.3 Atomistic Models of Radiation Defects in Silicon 

Reference, M. M. Sokol oski, "Structure and Kinetics 
of Defects in Silicon" NASA TN D-4154, 1967 
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mass action indicate that the formation of these defects might have the 
same variation with incident electron energy as that for the formation 
of single vacancies (av). This relationship has been verified experi- 
mentally.^'^ The V-P pair anneals rapidly near 150 C 3,3 ^ and the V- 0 
pair anneals rapidly near 350°C. 3 '^ The introduction rates (change in 
defect concentration per unit fluence) for these defects are in the 
range of 0.1 to 0.3 cm 1 for 1 MeV electrons. Since the calculated 
displacement rate is 5.2 cm~\ it appears that many of the vacancies 
are involved in other reactions at room temperature, such as recombination 
with interstitial atoms. 

The electron irradiation of p-type silicon at room temperature 

3 11 3 12 

results in a defect structure with net donor characteristics. ' ’ ’ 

This defect can donate an electron to (i.e,, accept a hole from) the 
valence band. The energy level of this donor defect is located 0.27 
to 0.30 eV above the top of the valence band. The room temperature 
introduction rate of this defect in silicon by 1 MeV electrons is 
roughly 0.03 cnT 3 . This value is considerably lower than those of 
defects found in n-type silicon. In addition, the introduction rate 
of this defect by 10 MeV electrons is about 16.5 times greater than 
that for 1 MeV electrons. 3 '^ Since the single displacement rate 
increases by only a factor of 2.5 with that electron energy increase, 
this defect appears to involve a more complex structure. It has been 
shown that defects involving the coupling of more than one vacancy will 
result in defects with introduction rates which increase more rapidly 

O 1C 117 

with electron energy than does the displacement rate. * ’ Two 

defect structures (divacancy 3 '^’ and "K" center' 3 ’^ 13 ), which have 

been studied by electron spin resonance techniques, may explain this 
behavior. These defects, shown in Figure 3.3, involve the coupling of 
two vacancies in each defect. Several attempts to determine the domi- 
nant recombination center in electron irradiated p-type silicon have 
yielded conflicting results, 3 '^ 3, 3 ' 3 ^’ 3 *^ 3, 3 ’^ Recent experiments 
have indicated that a defect with an energy level in the range of 0.27 

± 0.02 eV above the top of the valence band controls recombination in 

3 22 

electron irradiated p-type silicon. ' This conclusion is consistent 
with the known energy dependence of p-type silicon in that the diffusion 
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3 23 

length damage coefficient has been shown to vary with electron energy * 

in the same manner as the introduction rate of the E + 0.3 eV level 

3 1? v 

defect. ' 

The production of displacement damage in silicon by energetic protons 

is considerably different because the displacement cross sections are 

several orders of magnitude larger than those for fast electrons and vary 

rapidly with proton energy. The calculated silicon displacement cross 

-20 3 1 

section for 1 MeV protons is 3.5 x 10 cm 2 . * The average number of 

atomic displacements (v) resulting from such a primary displacement 
3 5 

is 4.8. ‘ Using equation 3.1.3, the displacement rate is found to be 
8500 cm - ^ for 1 MeV protons in silicon. The range of a 1 MeV proton in 
silicon is only 17.5 ym; therefore its energy and displacement rate will 
change rapidly after it enters a silicon crystal. The variation of the 
displacement rate with proton energy has been calculated by several 
authors. These results are shown in Figure 3.4. Although there are 
some differences in results, the displacement rate is proportional to 
(lnE)/E for protons of energies between 1 MeV and lOMeV. Above 
10 MeV, the various models differ as to the relative influence of Ruther- 
ford scattering, nuclear scattering, and inelastic processes of spalla- 
tion. Experimentally measured defect introduction rates for proton 
irradiation of silicon are less than one tenth of the calculated dis- 
placement rates. The defect energy levels in proton irradiated silicon 

3 29-3 31 

are those previously discussed for electron irradiated silicon. 

The proton damage, nowever, will be highly inhomogeneous because the 
secondary displacements occur near the site of the primary displacement. 

Neutron displacement damage in silicon is characterized by two 

important differences. The silicon displacement cross section for a 

-24 

1 MeV neutron is 2.4 x 10 cm 2 . This value is well below those for 
1 MeV protons and 1 MeV electrons. For this reason, the number of pri- 
mary displaced silicon atoms will be relatively small. The second 
difference involves the amount of energy transferred to the displaced 
silicon atom by the neutron. Since the 1 MeV neutron-silicon interaction 
is a hard sphere rather than coulombic collision, an average of about 
70 keV is transferred to the recoiling silicon atoms. The subsequent 
secondary collisions between silicon atoms will displace about 1500 
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DISPLACEMENTS (cm 



Figure 3.4 Displacement Rate in Silicon for Protons; 

Curve 1 (ref. 3.24), Curve 2 (ref. 3.25), 
Curve 3 (ref. 3.26), Curve 4 (ref. 3.27), 
Curve 5 (ref. 3.28) 
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silicon atoms. This displacement damage will be clustered near the site 
of the primary displacement. The defects identified in neutron irradi- 
ated silicon include those previously discussed for electron damage. 
Theoretical models of the neutron damage indicate that the high concentra- 
tion of electrically active defects in the cluster causes the center of 

3 32 

the cluster to behave as intrinsic silicon. ‘ This intrinsic silicon 
core is separated from the bulk silicon by a layer of space charge. 
Extensions of this model have been used to explain the majority carrier 


removal 
si 1 i con. 


and minority carrier recombination behavior of neutron irradiated 
3.33-3.35 


The main importance of the displacement defects produced by the 
irradiation of silicon solar cells is in their effect on the minority 
carrier lifetime of the silicon. In particular, the lifetime in the 
bulk p-type silicon of an n-p solar cell is the major radiation sensitive 
parameter. Since minority carrier lifetimes are inversely proportional 
to the recombination rates, the reciprocal lifetime contributions caused 
by various sets of recombination centers can be added to determine the 
inverse of the lifetime as follows: 



(3.1 .4) 


where t = minority carrier lifetime 

t q = minority carrier lifetime before irradiation 

t = minority carrier lifetime due to electron irradiation 
e J 

ip = minority carrier lifetime due to proton irradiation 

One of the most commonly used analytical tools for the determination 
of the particle type and energy dependence of degradation in silicon 
solar cells has been developed from the basic relationship for lifetime 
degradation: 


- = — + K $ (3.1 .5) 

T T o T 
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where x = final minority carrier lifetime 


x = initial minority carrier lifetime 
o 

$ = irradiation fluence 

- damage coefficient (lifetime) 

Minority carrier diffusion length is a more applicable and more easily 
determined parameter for solar cell analysis than minority carrier 
lifetime. Using L 2 = Ox, the above expression becomes: 

— = — + K. $ (3.1 .6) 

L 2 ,2 L 

o 


where L = final minority carrier diffusion length 


L q = initial minority carrier diffusion length 
<s> = particle fluence 

« L = damage coefficient (diffusion length) 

= K /D 

T 


When the fluence is sufficiently high so that L << L q we have: 

K l = 1/L 2 $ (3.1.7) 

If the L vs $ relationship exhibits the proper -1/2 slope, the damage 
coefficient, K^, can be used to uniquely define the particle type and 
energy dependence of silicon solar cell degradation. 


The minority carrier lifetime or diffusion length in an irradiated 
solar cell may be a function of the concentration of excess or nonequili- 
brium minority carriers present in the semiconductor. In solar cells, 
this oehavior is referred to as injection level dependence. This behavior 

is usually associated with solar cells damaged with high energy protons 

3 33 

or neutrons. Gregory ' has shown that the injection level dependence 
of lifetime in neutron irradiated solar cells does not follow classical 
predictions and has proposed a model based on the behavior of clustered 
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damage. The methods of measuring minority carrier lifetime or diffusion 
length often involve the injection of excess minority carrier concentra- 
tions which are many orders of magnitude smaller than those found in 
solar cells operating in space. Such methods are inadequate for the 
generation of data for the prediction of proton and neutron irradiated 
solar cell performance in space. 

3.2 Theory of Silicon Solar Cell Damage 

The basic solar cell equations can be used to describe the changes 
which occur during irradiation. This method would require data regarding 
the changes in the light generated current (i.e., I ), and data regard- 
ing changes in the series resistance, shunt resistance, and the basic 
diode parameters of saturation current and diode quality factor. 

Although such a method would be a logical analysis, most investigations 
have not reported enough data to determine the variations in the above 
parameters. The usual practice in the study of solar cell damage has been 
to reduce the experimental data in terms of changes in the cell short 
circuit current, open circuit voltage, and maximum power. 

It is also possible to characterize solar cell damage in terms of 
the changes in the minority diffusion length. Since the diffusion length 
can be measured experimentally and is a measure of the amount of displace- 
ment damage in the base of the solar cell, this method has been suggested. 
There are several practical and fundamental limitations to this scheme. 

The most serious limitations involve the evaluation of low energy proton 
damage in terms of diffusion length. Very low energy protons do consider- 
able displacement damage within the junction space charge region of a 
solar cell. This nonuniform damage increases the diode saturation current 
(I 0 ) and quality factor (n) by mechanisms which are not related to minor- 
ity carrier diffusion. This damage can cause serious reduction in solar 
cell V Qc without changing the cell diffusion length. In addition, the 
relation between diffusion length and the solar cell output parameters 
is not well defined, diffusion length is more difficult to measure than 
cell output parameters (particularly in the case of proton irradiated 
cells), and accurate measurement of diffusion length in thin or drift 
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field cells is extremely difficult. Because of these problems, methods 

have been evolved to evaluate solar cell radiation effects in terms of 

common engineering output parameters. Experience has shown that the 

variation of common solar cell output parameters during irradiation can 

be descirbed as shown for I in the following case: 

sc 3 



(3.2.1 ) 


The constant represents the radiation fluence at which the I behavior 
changes from constant to a direct function of the logarithum of the fluence. 
The constant C has dimensions of mA/decade fluence and represents the 
decrease in I per decade increase in radiation fluence in the logarithmic 
region. Although the above relationship is only empirical, there is 
some theoretical justification for the expression. Several observers 
have reported that the relation between the solar cell short circuit 
current and the diffusion length is as follows: 3 ’ 36 ’ 3 '^ 


I = A In L + B 
sc 


(3.2.2) 


The constants A and B are dependent upon the spectral content and inten- 
sity of the light source used to measure I . Tada has shown that the 
above expression is theoretically valid over a wide range of diffusion 
lengths for tungsten illuminations and to a lesser range under solar 

O OO 

illumination. A previously discussed relation, equation (3.1.6) can 
be transformed as follows: 



and substituted in equation (3.2.2). The resulting expression 


(3.2.3) 


'sc ■ 6 - ? 


In ( K.* + — i 

V 1 t»V 


has the same form as equation (3.2.1), 


(3.2.4) 


The variation of solar cell V during irradiation also may be 
empirically characterized by an expression similar to equation (3.2.1). 



- C' log 



(3.2.5) 
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In general the open circuit voltage of a silicon solar cell can be repre- 
sented by the following equation which was discussed in Chapter 1: 


I 


V = r ,n i 


sc 


(3.2.6) 


In using this expression, it is assumed that the saturation current (I Q ) 
is dominated by the diffusion component. In such cases the saturation 
current density is given by equation (1.2.5). If this expression is 
combined with equation (3.2.3), the following expression for the saturation 
current as a function of radiation fluence is obtained: 


I 


q D n S 
n n p 


/ 1\ 1/2 

f K 4> + — | 

V L L 0 V 


(3.2.7) 


where S is the cell area. Equations (3.2.4) and (3.2.7) can be substi- 
tuted into equation (3.2.6) to obtain the following expression: 



(3.2.8) 


The radiation fluence term ($) appears twice in the above expression. 

The fluence term in the numerator will have a much lesser effect on V QC 

than that in the denominator because it appears as the natural logarithm 

of the fluence rather than as the square root of the fluence. It appears 

therefore that the V variation with radiation fluence is dominated 

oc 

by the denominator of equation (3.2.8) and can be approximated by 
equation (3.2.5). 


The maximum power ( p max ) of a solar cell can be represented as the 

product of I , V , and a constant as follows: 
sc oc 


max 


= F • I 


sc 


oc 


(3.2.9) 


where F is the form or fill factor. The fill factor (F) is relatively 

insensitive to radiation which penetrates uniformly through a solar cell. 

For this reason, the variation of with irradiation is the same as 

max 
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that for the product of I and V QC> Equations (3.2.1) and (3.2.5) can 
be substituted into (3.2.9) and the resulting expression approaches the 
form of: 

P = P - C" log (] + (3.2.10) 

max maxo a l s> I ' 1 

Expressions of this form are found to closely describe the variation of 
P mav during irradiation. 

ITlaX 

3.3 The Concept of Damage Equivalence 

The wide range of electron and proton energies present in the space 
environment necessitates some method of describing the effects of various 
types of radiation in terms of a radiation environment which can be pro- 
duced under laboratory conditions. Since the changes in most solar cell 
parameters due to irradiation are in some way related to the minority 
carrier diffusion length, it is possible to determine an equivalent 
damage based upon this parameter. In Figure 3.5, the diffusion length 
changes are shown for lOft-cm, n-p silicon solar cells which have been 
subjected to several different types of irradiation. The results are 
described by equation (3.2.3) where the constant is dependent upon 
the radiation type. 

The concept of damage equivalence can alternatively be based on 
common solar cell parameters. The variation of short circuit current 
density for 10 ohm-cm n-p solar cells irradiated in various environments 
is shown in Figure 3.6. The J $c variation in each environment is 
described by equation (3.2.1). In this case two constants, C and <& , 
are required to describe the changes in d sc> Experience has shown that 
the constant C, under solar simulator illumination, does not vary 
greatly for different radiation environments. For electron irradia- 
tions in the 1 MeV and greater range, C is approximately 4.5 mA/cm 2 
decade. For proton and neutron irradiations, C approaches 6 mA/cm 2 
decade. For solar cells with the same starting J , the constant $ 

J v X 

is a measure of the damage effectiveness of different radiation environ- 
ments. The constant $ x for a particular radiation can be determined 
graphically at the intersection of the starting 0 sc and the extrapolation 
of the linear degradation region. 
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Figure 3.5 Variation of Solar Cell Diffusion Length with Fluence 

for Various Radiations 
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Since the value of $ is dependent upon the starting J value, 
it is not a good practical measure for relative damage effectiveness, 

It has been the practice to define an arbitrary constant referred to 
as the critical fluence ($ c ). One method of defining this value is that 
fluence which degrades a solar cell parameter 25% below its unirradiated 
state. Such a parameter is valid only when comparing cells with 
similar unirradiated parameters. To eliminate this problem, critical 
fluence may be defined alternatively as that fluence which will 
degrade a cell parameter to a certain value. 

By use of the critical fluence ($ c ) or the diffusion length damage 
coefficient (K^), it is possible to construct a model in which the various 
components of a combined radiation environment can be described in terms 
of a damage equivalent fluence of a selected monoenergetic particle. 

One MeV electrons are a common and significant component of space radia- 
tion and can be produced conveniently in a test environment. For this 
reason, one MeV electron fluence has been used as a basis of the damage 
equivalent fluences which describe silicon solar cell degradation. 

The use of the damage equivalent fluence scheme involves two sepa- 
rate problems. The first problem is to adequately describe the degrada- 
tion of an unshielded silicon solar cell under one MeV electron irradiation 
under laboratory conditions (i.e., normal incidence). The second problem 
is to reduce the effect of the space radiation environment (i.e,, continue 
ous energy spectrums of electrons and protons, isotropic incidence) on 
a shielded silicon solar cell to a damage equivalent fluence of one MeV 
electrons under laboratory conditions. 

3.4 One MeV Electron Irradiation of Silicon Solar Cells 

In this section the effects of one MeV electron laboratory irradia- 
tion of solar cells will be reviewed and discussed. Data will be presented 
which will form the basis for estimating solar cell performance, after 
the space radiation environment is reduced to a damage equivalent one MeV 
fluence. A very large volume of work has been reported concerning the 
effects of 1 MeV electron irradiation on silicon solar cells. This section 
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will limit the data to evaluations done with solar simulators. The data 
will also be limited to types of solar cells which are currently in 
common use on spacecraft. 

Currently n-p solar cells are in use as a primary power source on 

nearly all earth orbiting satellites. Variations in base resistivity 

and cell thickness cause significant differences in the response to 1 MeV 

3 23 3 39 

electron irradiation. ’ ’ ' Other variables such as the irradiation 

3 40 

temperature in the range of 200 to 370 K, * and p-type base dopant 

(boron vs. aluminum) have been shown to have little or no effect on the 

3 41-3 44 

solar cell response to radiation. 

The variation of n-p solar cell response with base resistivity has 

3 23 3 45 

been studied and reported for the range of 1 to 20 ohm-cm. ’ * 

Current n-p solar cell usage is confined to the ranges of 1 to 3 ohm-cm 
and 7-13 ohm-cm. Cells in the base resistivity range of 1-3 ohm-cm 
have greater initial maximum power output than cells in the 7 to 
13 ohm-cm range. The radiation hardness of n-p cells in the 7 to 13 ohm-cm 
range is greater than that of the 1 to 3 ohm-cm range, when the hardness 
is determined by parameters such as the critical fluence ($> ) or diffusion 
length damage coefficient (K^). As a result, 10 ohm-cm cells have greater 
maximum power output after a certain electron fluence is reached; however, 
the 2 ohm-cm cells produce greater maximum power at lower fluences. This 
crossover fluence depends upon cell thickness but is approximately 1 x 10 14 
one MeV electrons per cm 2 . 

Solar cell thickness has been shown to have a strong effect on the 

3 39 

output parameters of irradiated cells. ’ Cell thickness does not affect 
measures of inherent hardness such as the critical fluence (if properly 
determined) or the diffusion length damage coefficient. The thickness 
does, however, significantly affect the cell output parameters during 
the initial or low fluence stage of an irradiation. A set of graphs 
presenting J , V , P , and V n as a function of 1 MeV electron 
fluence are shown in Figures 3.7 through 3.14. Separate graphs are 
included for n-p cells in the 1-3 ohm-cm and 7-13 ohm-cm ranges. The 
changes in each cell parameter are shown for several solar cell thicknesses 
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Illumination Intensity, 28°C. 
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Figure 3.9 Voltage at Maximum Power Versus 1-MeV Electron 
Fluence for 1-3 ohm-cm N/P Silicon Cells. At 
135 mW/cm^ AMO Illumination Intensity, 28°C. 
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Figure 3.12 Open-Circuit Voltage Versus 1-MeV Electron Fluence for 
7-13 ohm-cm N/P Silicon Cells. At 135 mW/cm 2 AMO 
Illumination Intensity, 28°C. 
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Figure 3.13 Voltage at Maximum Power Versus 1-MeV Electron FI 
for 7-13 ohm-cm N/P Silicon Cells. At 135 mW/cnr 
Illumination Intensity, 28°C. 
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Figure 3.7a Normalized Short Circuit Current Density vs. 1-MeV Electron 
Fluence for 1-3 ohm-cm N/P Silicon Cells. At 135 mW/cm 2 
AMO Illumination Intensity, 28° C. 
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mumination, 28°C. 
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Figure 3.9a Normalized Voltage at Maximum Power vs. 1-MeV Electron Fluence 
for 1-3 ohm-cm N/P Silicon Cells. At 135 mW/cm^ AMO 
Illumination Intensity, 28°C. 
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Figure 3.10a Normalized Maximum Power vs. 1-MeV Electron Fluence for 1 
ohm-cm N/P Silicon Cells. At 135 mW/ cm2 AMO Illumination 
Intensity, 28°C. 
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Figure 3.11a Normalized Short-Circuit Current vs. 1-NleV Flectron 

Fluence for 7-13 ohm N/P Silicon Cells. At 135 mW/cm' 
Illumination Intensity, 28°C. 
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3.12a Normalized Open-Circuit Voltage vs. 1-MeV Electron 

Fluence for 7-13 ohm-cm N/P Silicon Cells. At 135 mW/cm' 
AMO illumination Intensity, 28°C. 
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Figure 3.13a Normalized Voltage at Maximum Power vs. VMeV Electron 

Fluence for 7-13 ohm-cm N/P Silicon Cells. At 135 mW/cm' 
Illumination Intensity, 28°C. 
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on each graph. These data are a composite of data from several sources, 

viq -3 aa.-? 47 

and unpublished TRW data. ’ ' It is emphasized that these 

data are intended to represent the mean behavior of current n-p solar 
cell production in the United States. These production practices involve 
use of quartz-crucible grown silicon, junction depths of approximately 
0.4 pm and absence of drift field collection. Solar cells produced with 
significant changes in the state of the art may require revision or 
addition to this typical radiation environment performance data. 

3.5 Effect of Electron Energy on Solar Cell Degradation 

The concept of damage equivalent 1 MeV electron fluence requires 
some method of evaluating the damage effectiveness of electrons of various 
energies. This effectiveness can be measured by the diffusion length 
damage constant (K^) or solar cell critical fluence ($> c ) for various elec 
tron energies. Experimental data have been reported for the electron 
energy range of 1 to 3 MeV®*^® and from 0.6 to 40 MeV.®’^® The results 
of these studies are in essential agreement and the results of reference 
3.23 are shown in Figure 3.15 (K, ) and Figure 3.16 ($ ) . In this case 
4> is defined as that fluence which degrades I to 19mA/crri under 

C r\ SC 

100mW/crn of tungsten il 1 imination. In both figures, data are shown for 

cells of various resistivities. Some important observations can be made 

from these data. The relative variations of the K. and $ with electron 

L c 

energy are identical. The relative variations of both parameters with 
cell base resistivity are also identical. On the basis of the experi- 
mental data, one can therefore define a relative damage effectiveness 
for each electron energy which will be a measure of the ratio of that 
electron fluence at a given energy to the 1 MeV electron fluence necessary 
to degrade an n-p solar cell to the same output parameter value. For 
instance, if a given 10 MeV electron fluence degrades a solar cell to a 
certain state of damage, then a 1 MeV electron fluence 16.5 times that 
of the 10 MeV electron fluence would be required to degrade the same 
cell to the same output conditions. This relationship will hold regard- 
less of whether 2 or 10 ohm-cm resistivity cells are under consideration. 

Wysocki reported data at 0.8 and 5.8 MeV which indicated that the 

3 49 

relative electron damage constant increased more rapidly with energy. 
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Figure 3,15 Electron Energy Dependence of K, Values for N on P 
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Figure 3.16 Electron Energy Dependence of $ Values for N on P 
Silicon Solar Cells 3 - 23 
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O c n 

Gorodetskii, et al . , reported data in rough agreement with references 

3.23 and 3.48 below 2 MeV, but indicate a much slower rise above that 

3 37 3 43 

energy. More recent studies by Bernard, et al . , ‘ and Lesbre * 
indicate good agreement with the results in references 3.23 and 3.48 
up to 3 MeV and 4.5 MeV, respectively. 

3.6 Effect of Proton Energy on Solar Cell Degradation 

The concept of damage equivalent 1 MeV electron fluence can be 
extended to the effects of proton irradiation. The problem is more 
complex, in the proton case, because the range of protons below 5 MeV 
is less than the thickness of a solar cell. For this reason, low 
energy protons produce nonuniform damage. This situation is further 
complicated by the fact that the damage produced per unit path length 
increases as the proton energy decreases. As a result, when a low 
energy proton is stopped in a solar cell, a large amount of damage is 
concentrated at the end of the proton track. 

When radiation damage is uniform throughout a solar cell, the 

relative effectiveness of various energy particles is the same when 

measured by the diffusion length damage coefficients, or critical flu- 

ences determined by cell parameters such as I, V _ or P , This 

sc oc max 

fact was graphically demonstrated by comparison of Figures 3.15 and 
3.16. In the case of protons with energies greater than 5 MeV, the 
damage to solar cells is relatively uniform. In this high energy range, 
the general concept of equivalency is directly applicable. At lower 
proton energies, the general concept of equivalency is not applicable; 
however, it can be used in a restricted manner as discussed below. 

Early experimental studies of the variation of damage in n-p 

silicon solar cells with higher proton energies indicated conflicting 

3 51 3 52 

results. The results reported by workers at BTL ’ and TRW ’ are 

shown in Figure 3.17 in normalized form. The major difference involves 

the behavior of the damage constant at proton energies greater than 10 

MeV. Recent experimental investigations have confirmed that the varia- 

3 53 3 54 

tion of damage in this proton energy range is very small. ’ ’ 

The results of these recent investigations are also shown in Figure 3.17. 

The degradation of n-p solar cells Irradiated with protons of ener- 
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gies below 3 MeV is more complex because of the nonuniform nature of the 

damage. Several experimental studies of low energy proton effects on 

3 41 3 55-3 61 

unshielded solar cells have been reported in the literature. ' ’ ’ 

Although there are some differences in the reported results, a few general 
observations can be made. Protons in the energy range from 1.5 to 3 MeV 
produce a maximum in relative radiation damage in silicon solar cells. 

The relative damage to silicon solar cell V and P due to low energy 
protons is more severe than that exhibited by the I . Proton damage 
in silicon solar cells can be normalized to the damage produced by protons 
of one energy. The proton energy employed for normalization of relative 
damage should be close to that producing maximum damage in space environ- 
ments, produce relatively uniform damage, and be available for laboratory 
evaluations. The use of 10 MeV proton damage is based on a compromise 
of the above requirements. The results of several studies of proton 

damage have been summarized in terms of relative silicon solar cell 

3 41 1 3 34 3 33 

damage as a function of proton energy. These 

relative damage results, normalized to 10 MeV proton damage, are shown 
in Figure 3.18. The Results in Figure 3.18 have been shown to hold 
for both 10 ohm-cm and 2 ohm-cm solar cells at proton energies greater 
than 10 MeV. 5,33 

It is emphasized that the results in Figure 3.18 are obtained by nor- 
mal incidence laboratory irradiation of solar cells from the front side. 

If similar data were prepared for normal incidence rear irradiations, the 
result would be similar for proton energies above 10 MeV. 3 ' 53 The effects 

due to rear incidence protons with energies below 10 MeV would be much 

3 6? 

lower than shown in Figure 3.18 ' The lower effectiveness occurs because 
rear incident low energy protons have insufficient range in silicon to 
cause atomic displacements in the space charge region of the solar cell. 

The variation of solar cell output parameters with 10 MeV proton 
fluence is described by equations (3.2.1), (3.2.5) and (3.2.10) in much 
the same way as is done for 1 MeV electrons. The values of the constants 
C, C', and C" tend to be somewhat greater than those found for 1 MeV 
electron irradiation. This value determines the decrease in solar cell 
output parameter per decade of radiation fluence. The fact that these 
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Figure 3.18 Relative Damage Coefficient for Proton Irradiation of N-P Silicon Solar Cells 


constants are somewhat different for electron and proton irradiation 
indicates that the concept of equivalency between the different 
types of radiation has limitations and is basically an approximation. 
This equivalence is further discussed in Chapter 6. 

3.7 Additional Effects of Low Energy Protons 

In addition to the low energy proton effects on unshielded cells 
discussed in the previous section, there are two aspects of 
low energy proton damage to be considered. These involve the effects 
of low energy protons on small unshielded gap areas on the front of 
solar cells and on unshielded backs of solar cells, 

When the ATS-1 and Intelsat II-F4 satellites suddenly exhibited 
degradations in power output of the order of 20% in weeks to a month 
after launch, the importance of low energy proton damage was dramati- 
cally demonstrated. Subsequent efforts related this anomalous degra- 
dation to the bombardment of narrow exposed surface areas of the solar 
cells by the intense low-energy proton fluence existing at synchronous 
altitude. The exposed areas resulted from slightly undersized or 
improperly applied cover slides which bared up to a 0.038 cm (15 mils) 
strip of solar cell surface. The high-intensity low-energy proton 
fluence, though incapable of penetrating the solar cell to a depth of 
more than a few microns, was able to produce junction damage which 
would shunt the power producing capability of the whole device. Exposed 
strips as narrow as 0.005 cm (2 mils) were sufficient to drastically 
alter the device's power producing capability. The absence of this 
effect in earlier solar array systems was attributed to shingling 
and overlapping adhesive. 

The results discussed in the previous section clearly indicated 

that low energy proton irradiation has an inordinately greater effect 

upon solar cell \f and as compared to similar irradiations with 
r oc max 

electrons or higher energy protons. The anomalous degradation of the 
ATS-1 and Intelsat II-F4 prompted many investigations into the effects 
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of low energy proton irradiation on partially shielded solar cells. 3 ' 63 " 3 * 67 

Curiously, Brucker and coworkers observed and reported this degradation 

effect in laboratory studies several months before the launch of 
3 63 

ATS-l. The results of these studies indicated that unshielded areas 
amounting to less than 1% of the total cell area can cause significant 
effects on cell power output. As a result of these studies, array 
manufacturers have taken measures to cover all areas of the silicon cell 
front surface with a cover slide and fill any gaps between the metallized 
base and cover slide with adhesive. 

The changes caused by the irradiation of small unshielded areas 
of solar cells with low energy protons can be explained in terms of 
solar cell theory. It was previously mentioned that the range of low 
energy protons in silicon is limited to less than the cell thickness. 
Particles which do not penetrate the cell produce defects only to their 
depth of penetration. This limited penetration results in unusual 
effects in the case of protons because lower energy protons produce 
more displacements per unit path length. The results of this behavior 
are shown graphically in Figure 3.19. In this figure, the calculated 
number of displaced silicon atoms per unit proton path is plotted as 
a function of depth in silicon for a 3 MeV proton (range 92.7 ym). It 
can be seen that the damage rises rapidly to a maximum near the end of 
the proton track. Every proton which is stopped in the silicon pro- 
duces such a damage peak at the end of its track. Protons which enter 
the silicon with energies of 0.5 MeV or less produce damage which is 
concentrated within a few microns of the cell surface. The space charge 
region of a modern cell extends from 0.4 to 1 micron below the cell 
surface. For this reason, low energy proton displacement damage is 
concentrated in the junction region. 

The entire solar cell junction can be considered to be an array of 
small parallel diodes, each having a characteristic described by the parallel 
combination of equations (1.2,3) and (1.2.6). Damage to only a small 
portion of this parallel diode array results in an increased effective 

•5 ri O CQ 

leakage or saturation current for the entire array. * ’ ' In Section 

1,2, the nature of the generation-recombination current was discussed. 
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Figure 3.19 Atomic Displacements as a Function of Depth for a 3 MeV Proton in Silicon 



The saturation current due to generation-recombination in the space charge 
region (equation 1,2.7) increases linearly as the carrier lifetime 
decreases (i.e., displacement damage increases) in the space charge 
region. The increased leakage current of a solar cell reduces the cell 
V Qc because of the relationship of V QC and I (junction leakage current) 
shown in equation (1.3.9). Since cell diode forward current ( I ) is 
increased at all voltages, the cell will also decrease because of 
low energy proton damage to small areas of the junction. 

This effect is illustrated In Figure 3.20, A partially shielded 
solar cell was irradiated with 3 x 10 13 p/cm 2 of 0.250 MeV protons. 

The current-voltage characteristics of this cell are shown before and 
after irradiation. The data indicate that the protons entered the 
silicon through a 0.0076 cm gap between the cover slide and the metal- 
lized bus strip. Although the I sc of the cell was unaffected by the 
irradiation, significant degradations occurred in V and P . Since 
solar cells are usually operated near the maximum power point, such 
changes have grave implications on in-flight performance. 

It has been observed in laboratory studies that the effects of low 

energy protons on small unshielded areas of cells produce a maximum in 

the degradation at a fluence of about 3 x 10 13 p/cm 2 , It has been suggested 

that the reversal of degradation is due to carrier removal effects, 3*61,3.68 

Considerable data exist regarding the effect of proton energy spectrum 

3 65 

and busbar-cover slide gap width on the degradation. ' Most reported 
laboratory studies have been confined to normal incidence proton irra- 
diations. 

In the past, solar cell usage has been confined to body-mounted solar 
cells on spinning satellites. Such applications provide a large measure 
of back shielding to a solar array. The requirements for increased 
spacecraft power and reduced weight have established trends toward the 
usage of oriented solar panels with minimal back shielding. Stofel has 
shown that low energy proton back side irradiation degrades silicon solar 
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Figure 3.20 Low Energy Proton Junction Damage, 0.250 MeV Protons, 

3 x 10 13 p/cm 2 , Partially Shielded N-P Solar Cell 
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3 61 3 65 

cells through carrier removal effects. ' ’ ' The use of thin 

soldered back contacts or other minimal back shielding should greatly 
reduce these effects. 

3.8 Effect of Temperature and Illumination intensity on Radiation 
Damage 

The commonly used solar cell output parameters (I , V QC , and P max ) 
vary with temperature in a near linear manner in the temperature range 
of -50° to +100°C. The expression for V QC (equation 1.3.9) varies 
directly with temperature and contains three other terms which vary 
with temperature (n, J^, and J^). The saturation current J Q increases 
rapidly as a function of temperature, and causes dV QC /dT to have a 
relatively large negative value. The actual value of dV 0C /dT for unirra- 
diated 2 to 10 ohm-cm cells varies from -2.4 to -2.1 mV/°C respectively 
or about -0.4%/°C regardless of resistivity. The temperature variation 
of I is influenced by temperature variation of the absorption coeffi- 
cient (a), and the minority carrier diffusion length (L) as indicated 
by equation (1.3,4). This causes a relatively small positive value of 
dI sc /dT, The normalized value of d I sc /dT found in unirradiated solar 
cells is approximately +0.05%/°C. The variation of current and voltage 
at P max behaves similarly to that of the I gc and V Qc , and because the 
temperature coefficient of V QC is large and negative, this behavior is 
reflected in the temperature coefficient of P mgx which is also large 
and negative. The normalized values of dP max /dT are approximately 0.6%/°C 
for unirradiated solar cells. If dl sc is normalized to the I sc , this 
parameter is invarient with illumination intensity. The normalized or 
fractional dV QC /dT decreases with illumination. 

There are very limited data concerning the variation of the above 

temperature coefficients with irradiation. Luft has reported that the 

dV / dT of silicon solar cells bombarded with 1 MeV electrons does not 

3 69 

change significantly during irradiation. Hayes and Ellis reported 
similar results for cells irradiated with 22 MeV protons or 2.4 MeV 
electrons. * The dI sc /dT, however, is increased significantly by 
electron and proton irradiation, Luft's data indicate a roughly two- 
fold increase in d I sc /dT between unirradiated cells and those irradiated 
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with 10 14 one MeV electrons/cm 2 , and a fourfold increase with TO 16 one 

MeV electrons/cm 2 . The data presented by Haynes and Ellis indicate that 

22 MeV proton Irradiations produce similar increases In dI sc /dT. The 

increase in dI sc /dT is due to the change in dx/dT (i.e., dL/dT) caused 

by irradiation. The absolute dP m= /dT of cells irradiated with electrons 

max 

or protons is decreased with fluence, however, when normalized to the 
unirradiated P max » it exhibits a small negative increase with fluence. 

Luft has also reported that the normalized or fractional temperature 

coefficient of I for both unirradiated and 1 MeV electron irradiated 
sc 

solar cells is independent of illumination intensity. The normalized or 

fractional temperature coefficient of V for both unirradiated and 1 MeV 

oc 

electron irradiated solar cells decreases linearly with increasing 

. , , . . . 3.69 

illumination. 

It should be emphasized that the above observations are based on 
"after the fact" measurements on solar cells which were irradiated at 
approximately 28°C. For this reason, they are not an indication of the 
effect of irradiation temperature on the production of radiation damage. 

A study at TRW indicated that the diffusion length damage constant (K^) 
of n-p solar cells is independent of 1 MeV electron irradiation tempera- 
tures between -&0°C and +130°C, but increases significantly at lower 
3 40 

temperatures. ' Workers at NRL have reported data which confirmed an 
increase in K L at low temperatures. 3,70 ’ 3,71 ’ 3,72 

A study of 1 MeV electron radiation effects on silicon solar cells 

under extremely low temperature and low illumination (Jupiter environ- 

ment) has been reported. * Debs and Hanes reported a study of 3 MeV 

proton damage to solar cells under Jupiter environment conditions. 

Their results indicate that n-p cells have higher starting efficiencies; 

however, after higher proton fluences, p-n 20 ohm-cm cells are more 
3 73 

efficient. 

Although illumination has been shown to effect the evaluation of 
radiation damage in silicon solar cells through injection level effects, 
it has been assumed that the production of displacement type radiation 
damage in silicon solar cells is independent of illumination intensity 
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during irradiation. Reynard has reported that during real time beta ray 

irradiations, silicon solar cells, illuminated and electrically loaded, 

degraded more severely than similar cells irradiated dark without load. 3 '^ 

3 75 

The results of a simile* study did not confirm the above result. ' 

Crabb recently reported that float zone silicon solar cells degraded 

with 1 MeV electrons exhibit a secondary degradation when illuminated with 

3 76 

a 10 sun illumination intensity. 

3.9 Effects of Neutron and Gamma Radiation on Solar Cells 

The radiation associated with nuclear weapons degrades solar arrays 
in the same manner as the radiation of the space environment. Solar 
array designers must allow for these effects when weapon events are 
included in the environment. The radiation from a weapon event is 
delivered at a much higher rate than space radiation. Because of these 
high radiation rates, other aspects of radiation effects become more 
apparent immediately following a nuclear radiation pulse. 

The most important aspect of neutron radiation on silicon solar 
cells is displacement damage which reduces the minority carrier lifetime 
in the same manner as protons and electrons. When silicon devices re- 
ceive neutron irradiation at room temperature, a large fraction of the 
displacement damage anneals within 100 seconds after the irradia- 
tion. The ratio of the damage produced to that which does not anneal is 

. 3 77 3 34 

referred to in the literature as the annealing factor. ' * ' 

Annealing factors larger than 10 have been reported. Such behavior 
is not suprising, because calculated displacement rates for various 
radiations are usually much greater than those found experimentally. 

The transient annealing of neutron damage is not an important considera- 
tion in the design of solar arrays;, however, the nonannealing component 
of neutron damage will contribute to the permanent damage produced by 
space radiation. This aspect of neutron damage has been studied by 
Brucker, 3 ’^ 8 Downing, 3 ^ 9 Morris?’® 0 Stofel, 3,8 ^ and Hicks. 3,82 Most 
of these studies utilized fission neutrons from nuclear reactors. If 
the fission spectrum of such reactors is averaged by weighting each 
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energy component by its theoretical displacement damage factor, the 
mean neutron energy Is very close to 1 MeV. The degradation of n-p silicon 

Q pi 

solar cell parameters with neutron Irradiation is shown In Figure 3.21. 

Neutron fluences may be converted to damage equivalent 1 MeV electron fluences 
by following expression: 

^ 1 MeV e 2000 x 45 1 MeV n (3.9.1) 

When neutron damage is evaluated with a solar simulator and described 

by equation (3.2.1), the constant C is approximately equal to 6 mA/cm 2 

per decade fluence. This value is significantly larger than that found 

for electron irradiation. Similar slope values are found in cells irra- 

3 33 3 81 

diated with high energy protons. Work by Gregory ’ and Stofel has 
shown that diffusion lengths measured in neutron-irradiated solar cells 
increase with the excess minority carrier concentration. This behavior 
is similar to that reported for proton-irradiated solar cells. 

Gamma ray radiation interacts with silicon mainly by the production 

of Compton electrons. These secondary particles have energies high 

enough to cause displacement damage in silicon solar cells. The effect 

3 84 

of gamma radiation on silicon solar cells has been reported by Fang 
and Hicks. 3 ' 82 The results of Cobalt 60 gamma irradiation of n-p silicon 
solar cells are shown in Figure 3.6. The permanent damage aspects of 
prompt gamma radiation from nuclear weapons can usually be neglected. 

The most important aspect of gamma radiation from weapons is the 
transient photocurrent generated in the array during a nuclear event. 

The primary photocurrent can be estimated from the following expression: 


Ipp = 6.4(yA cm" 3 rad"^ sec)* y • A ■ L (3.9.2) 

where y = dose rate (rad/sec) 

L = diffusion length (cm) 

A = cell junction area (cm 2 ) 

The transient rise and fall of the photocurrent has been treated by Wirth 

o oc 

and Rogers. The peak current values developed by solar cells under 
these conditions can be very large and may cause problems in circuits 
interfacing with the solar array. Current limiting by the external load 
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Figure 




and the internal cell series resistance may limit the observed photo- 
currents to values well below the generated current. Under very intense 
pulses of such ionizing radiation at room temperature the cell V Qc satu- 
rates at approximately 0.7 v . 3,86 " 3 - 88 This value appears to be related 
to the barrier potential (V b ) of the junction as determined by capacitance- 
voltage measurements. 

3.10 Lithium Doped Solar Cells 

Interest in this field began with Vavilov's report of a radiation 

3 89 

resistant diode made with lithium-doped, crucible grown silicon. 

Wysocki later reported lithium-doped solar cells which degraded under 

electron irradiation, but rapidly recovered at room temperature, 3 

Float zone silicon, with a characteristic lower oxygen concentration, 

was used to achieve this result. Subsequent work indicated that recovery 

also occurred in lithium-doped, quartz-crucible silicon solar cells. 

Since this initial work, the general subject was studied in two ways. 

Empirical changes in the manufacturing techniques for lithium-doped solar 

3 91 3 Q2 

cells were evaluated with the aim of optimizing the recovery effect. ' ’ * 

Other studies were directed at the development of a physical model of the 
degradation and recovery processes in lithium-doped silicon. 

Some of the more pertinent facts gained during these studies are 
as follows. The lithium concentration in a solar cell is not uniform, 
but increases in a linear or near linear manner with distance from the 
solar cell junction. This characteristic can be used to advantage to 
produce cells with exceptionally high open-circuit voltages. Solar cells 
with low or insufficient lithium concentrations do not recover in a satis- 
factory manner. Float zone silicon solar cells with exceptionally high 
lithium concentrations lose efficiency during storage in the unirradiated 
condition. These same cells, when irradiated and recovered, also exhibit 
a time-dependent loss of efficiency. This loss has been related to the 
room temperature diffusion of lithium into the active area of the cell. 

It has also been observed that higher lithium concentrations cause faster 
recovery rates. Because of the recovery rate dependence of the radiation 
damage in lithium-doped solar cells, it was difficult to evaluate cell 
performance by accelerator irradiations. Real time irradiations of 
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lithium-doped solar cells have been done with beta particle sources. 

The results of these beta irradiations indicated that some 
types of lithium-doped solar cells are slightly superior to n-p cells 
under some temperature conditions. The major potential advantages 
of lithium-doped solar cells over conventional n-D solar cells are in 
regard to proton ’ * u and neutron damage. * * * Figure 3.22 

shows that lithium doped solar cells are clearly superior to conven- 
tional cells. The long recovery period following a neutron exposure 
would probably be a severe limitation in military spacecraft. The 
most advantageous uses of lithium-doped solar cells would be for space- 
craft in proton dominated orbits with high proton fluxes. At present, 
such orbits are not commonly used. A summary of the current state of 

3 94 

the art in lithium-doped solar cells was recently published by Berman, 

3.11 Radiation Effects on Shielding Materials 

The degradation due to radiation effects on solar cell cover slide 
material in space is difficult to assess. The different radiation com- 
ponents of the environment act individually and synergistically on the ele- 
ments of the shielding material and also cause changes in the interaction 
of shielding elements. The complexity is illustrated in Table 3.2, where 
the various effects reported for commonly used cover materials are summar- 
ized and referenced. In addition to the data in Table 3.2, a large volume 
of data ha/e been presented in the literature regarding materials currently 
not in use for shielding solar cells. In this section, the emphasis will 
be on solar cell shielding material currently used in array construction. 

The cover glass shielding currently in use in most spacecraft con- 
struction is usually fabricated from Corning #0211 Microsheet or Corning 
#7940 fused silica. Where thin covers are desired, the usage tends toward 
Microsheet, because it is relatively inexpensive in thin sections. Where 
thicker covers are desired, Corning #7940 fused silica is used to avoid 
the darkening due to radiation. Cover glasses are always used with a 
MgF 2 antireflecting front coating and an ultraviolet rejecting filter on 
the rear surface. Cover glasses are usually attached to solar cells with 
silicone elastomers. 
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Figure 3.22 Recovered Power Output of Irradiated Conventional 
and Lithium Doped Solar Cells 3 ’ 53 ’ 3,93 





TABLE 3.2 

RADIATION EFFECTS ON SHIELDING MATERIALS 
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Most experimental assessments of radiation effects are based on 
accelerated testing in which a complete space environment is not simu- 
lated. This may account for some of the differences between darkening 
of cover glass material observed in laboratory radiation studies and 

space flight data for covered solar cells which indicated that radiation 

3 111 

effects in cover materials were insignificant. ' 

The radiation effects observed in cover materials can be charac- 
terized as ionization damage rather than displacement damage. In general, 
ionization effects are usually dependent upon the absorbed dose and to 
that degree are independent of particle type or energy. Some exceptions 
to this rule occur in the case of highly charged massive particles. In 

such cases, the ionization effects may be concentrated along the particle 

3 112 

track rather than uniformly distributed. * It is reasonable to assume 
that the ionization damage produced in cover materials by space electrons 
and protons is related to the total absorbed dose. This assumption allows 
the various radiation components of the space environment to be reduced 
to a total dose, without a laborious determination of degradation constants 
for each energy and particle. It also allows the use of experimental data 
from a single ionizing environment such as 1 MeV electrons. 

The most significant radiation effects in cover materials involve 

changes in the transmission of light in the visible and near infrared 

region. These data are commonly reported as spectral transmission data. 

The use of cover-glass spectral-transmission data in determining changes 

in solar cell output is rather cumbersome. This procedure was outlined by 
3 98 

Campbell. ‘ An alternate approach to the reporting of the data is the 

use of so-called "wide band" transmission loss. In this method, solar 

cell short-circuit currents are measured under sun simulated conditions, 

with cover slides attached. The cover slides are attached with a thin 

liquid film with an index of refraction (n = 1.4) similar to that of 

silicone adhesive. Cyclohexane and n-amyl alcohol have been used for 

this purpose. The "wide band" transmittance is defined as the solar 

cell I with an irradiated cover slide in place divided by the solar 
sc 

cell I with the unirradiated cover slide in place. Such measure- 
ments are influenced by solar cell spectral response. Results 
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determined with unirradiated solar cells will not be representative of 
those for irradiated solar cells. This error is probably negligible 
compared to the uncertainty of the available experimental data. 

Since the "wide band" transmission loss is a measure of the loss 

in light transmitted, it directly affects the light generated current 

(I L ) and likewise the short circuit current (I sc )- It is desirable to 

use the "wide band" transmission data to estimate the change in solar 

cell P max - Equation (3.2.9) indicates that cell P max is proportional 

to the product of I and V . Because V is proportional to In I , 

sc oc oc sc 

the following relation can be developed to estimate the change in P mgx 
due to cover slide darkening from transmission data: 


P 

max 

P 

maxo 


1n ' W 

m Ci sc ) J 


(3.11.1) 


where P m ,„/P„„„„ - the fractional change in P 

max' maxo 3 max 

T * the "wide band" transmission of irradiated 
cover glass 

I = the short circuit current of cell with unirradiated 
cover glass 

To aid in the estimation of solar array losses due to reduced trans- 
mission from radiation effects in cover slide materials, data relating 
transmittance to absorbed dose is required. In Figure 3.23, "wide band" 
transmittance is shown for various absorbed doses. The absorbed doses 
were produced by 1 MeV electron irradiations in a room temperature, air 
environment which included no ultraviolet illumination. This electron 
radiation is sufficiently penetrating to produce a relatively uniform 
dose through the entire cover slide, coating, and filter. The P max /P max0 
data shown in Figure 3.23 was calculated from the "wide band" transmittance 
value by use of equation (3.11.1). The data in Figure 3.23 include 
0.0152 cm (0.006 in.) 7940 fused silica and 0211 Microsheet cover slides 
with anti reflecting coating and blue filter. It is an established fact 
that Corning #7940 fused silica exhibits little or no darkening due to 
radiation in the visible region. Since the transmission loss for 7940 
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cover glasses must be assumed to be due to changes in the filter, the 
data can also be used for thicker cover slices. For thicker 0211 Micro- 
sheet cover glass, the data in Figure 3.23 cannot be used. 

The dose-depth profiles experienced by cover glass shielding in 
space are highly non-uniform due to the low energy protons stopped in 
the front surface, To accurately estimate the transmission through a 
cover glass with such a dose-depth profile, would require the integra- 
tion of absorption coefficients (as a function of dose) through the 
cover glass and its thin film layer. The lack of absorption coefficient 
data for these materials for various doses in a total space environment 
does not allow such evaluations at this time. 

The diversity of technical opinions on transmission loss in cover 
glass due to space radiation also includes those who do not include this 
factor in array power estimates and those who simply allow for a 2% to 
4% initial loss due to cover glasses and adhesive darkening due to radia- 
tion and ultraviolet effects. Recent studies by Luedke at TRW indicated 
that nearly all darkening produced in 0211 Microsheet by a dose of 10 7 
rad(Si0 2 ) was bleached by a relatively short ultraviolet light exposure. ^*107 
Such results indicate that the use of data such as that in Figure 3.23 is 
probably an overly conservative practice and emphasizes the importance of 
performing cover glass darkening studies in a realistic environment. 

Some investigations have reported results which indicate that cerium 
doping of glass reduces or eliminates darkening due to irradiation. 3 - 104 » 3,113 

Other studies indicated that hydrogen impregnation of glasses reduces 

3 114 

transmission losses due to irradiation effects. 
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CHAPTER 4 


4.0 RELATIVE DAMAGE COEFFICIENTS FOR SPACE RADIATION 

A large volume of experimental data is available for normal incidence 
irradiation of unshielded solar cells. These data are not directly applic- 
able in the prediction of space radiation effects because of the omni- 
directional nature of the space radiation and because of the energy 
degrading effects of cover glass (shielding). In this section, the 
analytical methods of calculating the damage effectiveness of each com- 
ponent of the space radiation will be detailed. The damage effectiveness 
of space radiation is calculated relative to normal incidence 1 MeV elec- 
trons and 10 MeV protons on unshielded solar cells. This concept of the 
damage effectiveness or relative damage constant (D) is an extension of 
the previously discussed concept of equivalent fluence. It will allow 
the reduction of all components of the space radiation to an equivalent 
laboratory (normal incidence, monoenergetic) irradiation. In this way, 
laboratory data can be used to predict the behavior of shielded solar 
arrays in space. In addition, the similar problem of calculating energy 
deposition at various depths in shielding will be discussed. 

4.1 Geometrical Aspects of Radiation Fluences 

An omnidirectional flux is defined as the number of radiation par- 
ticles of a particular type and energy which isotropically traverse a 
test sphere of unit cross-sectional area per unit time. The commonly 
used sources of space radiation literature tabulate the environment in 

O _ I 

terms of omnidirectional fluxes with units of particles cm day . A 

commonly repeated derivation in the literature regarding the conversion 

4 1 

of omnidirectional fluxes to unidirectional fluxes is as follows. 

Assume a unit of plane area in space with an incident omnidirectional 
flux of particles. 

$ = the component of the omnidirectional flux which is 
n normal to a surface 

$ = the omnidirectional flux 
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4tt = solid angle of test sphere (steradians) 
e = angle of radiation incidence (from normal) 
dfi = an increment of solid angle 

= 2n sinede (for rotational symmetry) 
cos e = projected area of unit plane area 


$ 


n 



cos e dn + 



(4.1.1) 



2n sin e cos e de + 



sin e cos e de 



The above derivation implies that the unidirectional fluence is equal in 
intensity or "equivalent" to the omnidirectional flux divided by 2. 
Likewise, if the unit plane area has infinite back shielding (i,e., 
integrate e from o to -n/2 only), one-fourth of the omnidirectional 

fluence is equal to the intensity of the unidirectional normally inci- 
dent fluence. The above expression determines the normal component of 

an omnidirectional flux. The conversion of an omnidirectional flux to 
an equivalent unidirectional flux must properly weight the damage 
effectiveness of all angular components. 

The expression for the effectiveness or relative damage constant, 
weighted for all angular components of an omnidirectional monoenergetic 
flux and assuming infinite back shielding, is as follows: 

/•it/2 

D(E,t) * i I D(E , 6 ) • 2tt sine de (4.1 -2) 

c * / o ° 

where D(E,t) = relative damage coefficient of omnidirectional 
radiation particles with energy E, relative to 
unidirectional 1 MeV electrons or 10 MeV protons 

D(E 0 , e) = damage coefficient of unidirectional radiation particles 
with angle of incidence (e) and energy (E 0 ) relative 
to unidirectional 1 MeV electrons or 10 MeV protons 


t = shielding thickness-, for the case of t=0, E=E 0 
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The quantity 2ir sin e de is an increment of solid angle as in equation 
(4.1.1). Equation (4,1.2) must be further modified to reflect the energy 
degradation in the cover glass shields used on silicon solar cells (t^ 0). 

4.2 Effect of Shielding on Radiation 

A common solar cell configuration involves infinite back shielding 
and an optically transparent finite shield covering the front surface of 
the cell. The assumption of infinite back shielding is not always valid, 
and the differences in both shield thickness and material require separate 
treatments for front and back radiation. If an omnidirectional flux of 
radiation particles with energy E is incident on a solar cell shield 
of thickness t, the particles not stopped in the shielding will exit the 
shielding (i.e., enter the silicon) with an energy of E . The energy E q 
will be a strong function of the angle of incidence because of varying 
path length in the shield. The particle track length in the shield is 
equal to t/cos e. By subtracting the particle track length in the shield 
(t/cos e) from the range of the particle, R(E), in the shield material, 
one can determine the residual range, RCE q ) * of a particle with energy E Q . 
Thus: 

E 0 (E,e,t) - R ' 1 [r(E) - asl-o ] 

where R is a convenient form used to represent an inverse function 
of the range-energy relation R. Proton and electron range-energy data 
suitable for this calculation have been conveniently tabulated by 
Janni^'^ and Berger and Seltzer.^’ 

4.3 Electron Space Radiation Effects 

The evaluation of D(E,e) is necessary to complete the integration of 
equation (4.1.2). The data regarding the experimental evaluation of the 
relative damage coefficient for n-p silicon solar cells, D(E) for various 
electron energies at normal incidence is presented in Figure 4.1 (dashed 
line). Electrons in the MeV energy range penetrate silicon solar cells 
thoroughly enough that the damage produced by an electron can be considered 
uniform along its track. For this reason, the amount of displacement 
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damage produced by a high energy electron is proportional to the total 
track length produced in a solar cell, and hence: 

D(E Q ,e) oc (particle track length) • (projected cell area) (4.3.1) 

The length of an individual electron track in a solar cell is proportional 
to sec 0 or 1/cos e. The number of electrons intercepted by the cell is 
proportional to its projected area normal to the direction of the radia- 
tion or cos e. The net result of these two factors on equation (4,3.1) 
is cancellation of the terms involving the angle of incidence (0), and 
D(E,e) is shown to be independent of 0 or equal to D(E), The fact that 

fast electron damage of unshielded silicon solar cells is independent of 

4 5 

the angle of incidence was experimentally confirmed by Barrett. * 

Equation (4.1.2) for the case of electron space radiation can be 
modified to the following expression: 

D(E,t) = j I D(E ,0) 2 tt sin 0 d0 (4.3.2) 

Equation (4.3.2) can be evaluated with the aid of equation (4.2.1) to 
evaluate E Q and the data in Figure 4.1 to evaluate D(E Q ,0). The integra- 
tion of equation (4,3.2) has been performed by machine and the results 
are also shown in Figure 4.1. The results are also tabulated in Table 
4.1. Because of electron straggling, there might be some question re- 
garding the suitability of equation (4.2.1) to determine E Q , however use 
of alternate Monte Carlo methods yielded results identical to those in 
Figure 4.1. Rosenzweig published similar space electron damage factor 
curves.^*® Barrett also published a similar analysis based on the diffu- 
sion length damage coefficient and empirically fitted analytical expres- 
sions to the data. 4-5 

The evaluation of ionization dose in solar array materials due to 
omnidirectional space electron fluences is analogous to that just com- 
pleted for silicon solar cell degradation. In the case of absorbed dose, 
the energy deposited by the radiation in the shielding is determined in 
terms of rads or joules per kilogram. To evaluate this energy deposition 
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1. 51E-02 

2.75E-03 

0. 

0. : 

.45 

7.206-14 

1.226-01 

7.806-02 

5.38E-02 

2. 526-02 

7.22E-03 

3.00E-04 

0. 

.50 

3.00E-14 

1.47E-C1 

9. 996-02 

7.27E-02 

3.826-02 

1.41E-02 

2.11E-03 

0. 

.55 

3 . 8 1 E - 1 4 

1 . 75E-01 

1.24E-01 

9.386-02 

5.41E-02 

2.36E-02 

5.89E-03 

0. 

.60 

9 .6 IE- 14 

2.05E-01 

1. 50E-01 

1.176-01 

7. 23E-02 

3.60E-02 

1. 19E-02 

0. 


1 .046-13 

2.36E-01 

1.77E-01 

1.416-01 

9. 14E-02 

4.98E-02 

1.98E-02 

0. 

*1 

1 .! 2E-1 3 

2.69E-01 

2.056-01 

1.666-01 

1. 126-01 

6.56E-02 

3.00E-02 

7.896-06 


1 .266-13 

3.40E-01 

2. 696-01 

2.246-01 

1. 60E-01 

1.036-01 

5.61E-02 

2.416-03 

.90 

1.44E-13 

4. 17E-01 

3.38E-01 

2.876-01 

2. 13E-01 

1 .45E-01 

8.79E-02 

l. 016-02 

1 .00 

1 .6CE-13 

5. C0E-G1 

4. 146-01 

3.57E-01 

2. 74E-01 

1.95E-01 

1.27 E-0 1 

2.426-02 

1.10 

1 .76E-13 

5.836-01 

4.91E-01 

4.30E-01 

3.376-01 

2.49E-01 

1.69E-01 

4.3 16-02 

1 .20 

1 .926-13 

6. 70E-0 1 

5. 736-01 

5.076-01 

4.066-01 

3.08E-01 

2.18E-01 

6.71E-02 

1.40 

2.24E-13 

8.59E-01 

1 

7.50E-01 

6.756-01 

5. 596-01 

4.42E-01 

3.32E-01 

1.316-01 

1 .60 

2. 56?-! 2 

1 . 06E+0C 

9.39E-01 

3.56E-01 

7.256-01 

5.91E-01 

4.62E-01 

2.13E-0 1 

1.80 

2.386-13 

1.266+00 

1. 13E+00 

1. 046*00 

9. COE-O i 

7.51E-01 

6.04E-01 

3.116-0 1 

2.00 

3-206-13 

1 • 476 +CC 

1. 34E+00 

1.246+00 

1. 096 + 00 

9.23E-01 

7.60E-01 

4.246-0 1 

2.25 

3.606-13 

1 . 7 ? E + 0 0 

l. 59E+00 

1 .49E+00 

1. 32 E+0 0 

1.14E+00 

9.61 E-0 1 

5.766-01 j 

2.50 

4.00E-13 

2.0CE+00 

1.856+00 

1.746+00 

1. 57E+00 

1.37E+00 

1.17E+00 

7.466-01 

2.75 

4.40E-13 

2.27E+00 

2. 126+00 

2 .016+00 

1. 92E + 00 

1.61 E+00 

1.40E+00 

9.30E-01 

3.00 

4.30F-13 

2.5 IE + 0C 

2. 376*00 

Z.26E+00 

2. C7E+00 

1.86E+00 

1 . 64 E+00 

1.136+00 

3.25 

5.206-13 

2.76E+00 

2.616+00 

2 .496+00 

2. 30E+00 

2.08E+00 

1.86E+00 

1.33E+00 

3.50 

5.606-1 3 

3. QQE+00 

2. 956+00 

2. 73E+00 

2. 53E+00 

2.31E+00 

2.07E+00 

1.53E + 00 

3.75 

6.00c” 13 

3. 25E+0C 

3.09E+0 0 

2 .976 + 00 

2. 77E+00 

2.54E+00 

2.30E+00 

1.73E+00 


6.4CF-13 

3.50F+00 

3. 34E+00 

3 .226+00 

3. 01E+00 

2.77E+00 

2.52E+00 

1.93E+00 


6.30E-13 

3;?2E+00 

3.576+00 

3 .456+00 

3.24E+00 

3. 006+00 

2.74E+00 

Z.13E+b0 


7.20E-13 

3.95E •*•00 

3. 796+00 

3 .676+00 

3.46E+00 

3.22E+00 

2.96E+00 

2.33E+00 

4.75 

7 . 6 C E - 1 i 

4 . i 7E + QC 

4. 026+00 

3 .896+00 

3. 68E+00 

3.44E+00 

3.17E+00 

2. 536+00 

5.00 

3.006-13 

4.40E+0C 

4. 24E+00 

4.12E+00 

3. 90E+00 

3.65E+00 

3.39E+00 

2 .746+0 0 

5.25 

3.41E-13 

4 • 626+00 

4.46 E+00 

4.346+00 

4. 12E+00 

3.37E+00 

3.60E+00 

2.94E+00 ! 

5.50 

8 .8 IE- 1 3 

4.84E+00 

4.68E+00 

4.556+00 

4.34E+00 

4.08E+00 

3.81E+00 

3.14E+00 

5.75 

9.216-13 

5.C6E+00 

4. 90E+00 

4.77E + 00 

4. 556 + 00 

4.30E + 00 

4.02E+00 

3.34E+00 

6.00 

9.61E-13 

5.286+00 

5. 12E+00 

4.99E+00 

4. 776+00 

4.51 E+00 

4.23E+00 

3.54E+00 

6.2 5 

l.OCE-12 

5. 49E + 00 

5. 33 E + 00 

5.20E+00 

4. 57E+00 

4.71E+00 

4 .43 E+00 

3.73E+00 

6.50 

1 .34E-12 

5.71E+C0 

5.55E+00 

5.41E+00 

5. L9E+00 

4.92E+00 

4.63 E+00 

3.92E+00 

6.75 

1 • 0 8 E - 1 2 

5. 93E+00 

5. 776+00 

5 . 63E +00 

5. 406+00 

5.1 3E+00 

4. 84 E+00 

4.11E+0C 

7.00 

1.126-12 

6. 15E+0C 

5.98E+00 

5.856+00 

5.61E+00 

5.34E+00 

5.04E+00 

4.31E+00 

8 .00 

1.2 BE- 12 

6 .886+00 

6.73E+00 

6 • 60E+00 

6.38E+00 

6.12E+00 

5.83E+00 

5.08E+0C 

9.00 

i . 4 4E- 1 2 

7. 606+00 

7.446+00 

7.32E+00 

7.09E+00 

6.83E+00 

6.536+00 

5.786+00 i 

10.00 

t .60E- 12 

9. 30E+0G 

8. 156 + 00 

8.02E+00 

7. 79E+00 

7.52E+00 

7.22 E+00 

6.46E+00 f 

15.00 

2.40E-12 

1. C6E+01 

1. 05E+Q1 

L .046+01 

1.02E+01 

1 .006+0 1 

9. 75 E+00 

9.06E+00 

20.00 

3.2CE-12 

1.23E+01 

1.236+01 

1.22E+01 

1. 20E+01 

1.18E+01 

L.16E+01 

1.10E+01 | 

2 5.00 

4.00E-12 

1.36E+01 

1.35E+01 

1.35E+01 

1.336+01 

1.3 1E+01 

1.296+01 

1.24E+01 

30.00 

4.8CE-12 

1.47E + 01 

1. 476+01 

1.466+01 

1.45E+01 

1.43E+01 

1.41E+01 

1.36E+01 i 

40.00 

6.40F-12 

1.65E+01 

1.64E+01 

1.64E+01 

1.62E+01 

1.61E+01 

1.59E+01 

1.546+01 t 


Table 4.1. Electron Damage Coefficients 
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at various depths in the shielding, an expression similar to equation 
(4.3.2) can be used. Equation (4.3.2) is modified to the extent that 
the electron stopping power 


l- 

\ p ^/collision 

replaces D(E q ), and D(E,t) becomes the absorbed dose per unit fluence. 

The results of this integration are shown in Figure 4.2 and in Table 4.2. 

4 6 

Rosenzweig has published similar curves. 

4.4 Proton Space Radiation Effects 

For proton space radiation, the evaluation of equation (4.1.2) is 
more complex than that previously discussed for electrons. Two problems 
arise in the treatment of space protons with energies less than about 
10 MeV, because of their limited penetration and increased damage produc- 
tion. One problem exists because the relative damage constants based 

on silicon solar cell I , V , and P are different and diverge at 

sc oc max 

low proton energies. The second problem is that low energy proton damage 
has been experimentally characterized only for normal incidence irradia- 
tion, and basic considerations indicate that the damage is a strong func- 
tion of the angle of incidence. The normal incidence proton coefficients 
for energies of 10 MeV and greater can be assumed to be independent of 
the angle of radiation incidence for the same reasons discussed for elec- 
tron irradiation in the previous section. 

The physical distribution of low energy proton damage was discussed in 
section 3.6. The most significant aspect of the low energy proton damage 
is the fact that a majority of the displacements are produced at the end 
of the proton track, as illustrated in Figure 3.19. The high damage con- 
centration near the end of the proton track allows the construction of 
a simple damage model for the prediction of the effect of angle of inci- 
dence on low energy proton damage in silicon solar cells. It is assumed 
that the effect of a low energy proton, of arbitrary angle of incidence 
and energy, is roughly equal to that of a normally incident proton with 
a range equal to the perpendicular penetration of the non-normal incident 
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ELECTRON STOPPING POWER 

RAOISILICONI/UNIT OMNIDIRECTIONAL FLUX 



SHIELO THICKNESS, GM/CM2 


1.84E-02 3.69E-02 7.38E-02 1.23E-01 


1.84E-01 3.6SE-01 


9.616-15 

1.28E-14 

1.606-14 


2.296-08 
2. 86E-08 
2.67E-09 


4.406-14 

4.8CE-14 

5.28E-14 


1.60E-0e 

1.54E-C8 

1.49E-08 



I.92E-14 

2.44E-08 

1. 57E-08 

0. 

0. 

2.40E-14 

2. I5E-08 

1.83E-08 

0. 

0. 

2.72E-14 

2.01E-Q8 

1. 8 IE- 08 

9.48E-09 

0. 

3.20E-14 

1.85E-08 

1. 74E-08 

1.31E-08 

0. 

3*6 

1.74E-08 

1.67E-08 

1.39E-08 

0. 

*.00£-14 

1 • 66E-Q8 

1.61E-08 

1 •<» IE-08 

4. 81E-09 


1.55E-08 

1.5LE-08 

1.46E-0B 


1 .40E-08 
1 • 396-08 
1 .366-08 


8.046-09 
9.51E-09 
1. 04E-08 


5.76E-14 

6.4QE-14 

7.2QE-14 

1.44E-08 

1.40E-08 

1.35E-0E 

1.426-08 

1.376-08 

1.33E-08 

1.34E-08 
1.31E-08 
1 .28E-08 

1.08E-08 
1. 11E-08 
1.12E-08 

4.756-09 

7.07E-09 

8.34E-09 

0 . 

0 . 

1.54E-09 

8.00E-14 

1.32E-08 

1.30E-08 

1 .256-08 

1. 13E-08 

9.02E-09 

4.95E-09 

8.816-14 

1.30E-08 

1.27E-08 

1 .23E-08 

1. 12E-08 

9.40E-09 

6.38E-09 

9.6 IE-14 

1.27E-C8 

1.25E-08 

1.22E-08 

1. 12E-08 

9.64E-09 

7.246-09 


1.04E-13 

1.12E-13 

1.28E-13 


L.92E-13 

2.24E-13 


1.26E-08 

1.256-08 

1.236-08 


L.21E-08 

1.20E-08 


t ) 7 

1,60 

1.80 

2.00 

2.2S 

2.50 

2.75 



2.366-13 
2.88E-13 
3.2QE-13 
3.60E-I 3 
4.00E- 13 
4.40E-13 
4.80E-1 3 


5.2QE-13 

5.60E-13 

6.0CE-13 


1.21E-08 

1.21E-08 

1.22E-08 

1.22E-08 

1.23E-08 

1.24E-08 

1.24E-08 


1.25E-0e 

1.26E-08 

1.26E-08 


. 24E-08 
1.23E-08 
1.21E-08 


8 

8 

8 


1.19E-08 

1.196-08 


I.19E-08 

L.20E-08 

1.20E-08 

I.21E-08 

1.22E-08 

1.23E-08 

1.23E-08 


L.24E-08 

1.25E-08 

I.26E-08 


1.20E-08 

1.19E-08 

1.18E-08 


1.16E-08 

1.17E-08 


I. 12E-08 
l.UE-08 
1.11E-08 


8 

8 


l.UE-08 

I.I2E-08 


9.78E-09 

9.90E-09 

1.01E-08 


1.06E-08 


7.76E— 09 
8.I4E-09 
8.64E-09 


1.17E-08 
I.18E-08 
1 .18E-08 
l . 19E-0B 
1.20E-08 
L« 2 IE -08 
1.22E-08 


1.23E-08 

I.24E-06 

1.24E-08 


I. 13E-08 
1. 14E-08 
1. 15E-08 
1. 16E-08 
1. 17E-08 
1. 186-08 
1.19E-08 


1. 20E-08 
1. 21E-08 
1.22E-08 


1 .0 7E-0 8 
I.09E-08 
1.10E-08 
1.1IE-0B 
1.13E-08 
1. 146-08 
1.16E-08 


1.17-6-0 8 
1.186-08 
1.19E-08 


1.00E-08 

1.02E-08 

1.04E-08 

1.06E-08 

1.08E-08 

1.10E-08 

l.UE-08 


1.13E-08 

I.14E-08 

1.15E-08 


I.04E-12 

I.08E-12 


1. 31E-08 
1.32E-C8 
1.32E-C8 


1.31E-08 
1.3 IE- 08 
I.31E-08 


1.30E-08 
1.30E-08 
1 .31E-08 


1.28E-08 
1. 29E-08 
1. 296-08 


1.266-08 

1.26E-08 

1.27E-08 


1.23E-08 

1.24E-08 

1.246-08 



0 . 

0 . 

2.86E-09 


7.46E-09 


7.97E-0 9 
8.37E-09 
8.71E-09 
9.05E-09 
9.35E-09 
9.6 IE-09 
9. BSE— 09 


1.01E-08 

1.02E-08 

1.04E-08 


6.40E-13 

1.27E-08 

1.266-08 

1.25E-08 

1. 23E-08 

1.20E-08 

1.16E-08 

1.066-08 

6.8CE-13 

1.286-08 

1.276-08 

1 . 26E-08 

1.23E-08 

1 .216-08 

1.17E-08 

1.076-0 8 

7.206-13 

1.28E-08 

1.27E-08 

1 . 266- 08 

1.24E-08 

1.21E-08 

1.18E-08 

1.09E-08 

7.6CE-13 

1.29E-08 

1.28E-08 

1.27E-08 

1. 25E-08 

1.22E-08 

1.19E-08 

1.106-08 

8.00E-13 

1.29E-08 

1.28E-08 

1.27E-08 

1.25E-08 

1.23E-08 

1.20E-08 

l.UE-08 

8.416-13 

1. 30E-08 

1.29E-08 

1.28E-08 

1, 26E-08 

1.24E-08 

1.216-08 

1.12E-08 

8.816-13 

1.30E-08 

1.29E-08 

1.286-08 

1.27E-08 

1.24E-08 

1.21E-08 

1.136-06 

9.216-13 

1.306-08 

1.30E-08 

1.296-08 

1.27E-08 

1.25E-08 

1. 226-08 

1.14E-08 

9.616-13 

1.31E-08 

1.30E-08 

1.29E-08 

1.28E-08 

1.25E-08 

1.23E-08 

1.156-08 


1.16E-0 8 
1.17E-08 
1 . 1 7E-08 



1.126-12 

1.326-08 

1.32E-08 

1.31E-08 

1.296-08 

1.27E-08 

1.256-08 

1.166-08 


1.286-12 

1. 346-08 

1.33E-08 

1.326-08 

1.316-08 

1.296-08 

1.27E-08 

1.216-08 


1.446-12 

1.35E-08 

1.346-08 

1.34E-08 

1.326-08 

1.31E-08 

1.296-08 

1.23E— 0 8 

10.00 

1.606-12 

1. 36E-08 

1.35E-08 

1.35E-08 

1. 33E-08 

1.32E-08 

1.30E-08 

1.256-08 

15.00 

2.40E-12 

1. 39E-08 

1.39E-08 

1.386-08 

1. 38E-08 

1.36E-08 

1.356-08 

1.316-08 

20.00 

3.206-12 

1.426-08 

1.426-08 

1 . 41E-08 

1. 406-08 

1.40E-08 

1.38E-08 

1.35E-08 

25.00 

4.006-12 

1.44E-08 

1.43E-08 

1.43E-08 

1.426-08 

1.42E-08 

1.416-08 

1.38E-08 

30 .00 

4.806-12 

1.45E-08 

1.456-08 

1.45E-08 

1.446-08 

1.436-08 

1.426-08 

1.406-0 8 

40.00 

6.406-12 

1.47E-08 

1.47E-08 

1.47E-08 

1.466-08 

1.466-08 

1.45E-08 

1.436-08 


Table 4.2. Electron Stopping Power, Rad/Unit Omnidirectional Flux 
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proton. To partially correct the inaccuracies of this proposed model, 
a factor is employed which relates the ratio of the total displacements 
produced by the non-normally incident proton to those of a normally inci- 
dent proton with the same perpendicular penetration in the silicon solar 
cell. The low energy proton relative damage coefficient given by the 
above model can be expressed as follows: 

N td (E ) cos e 

0(E o ,e) = D(E n’ 0) ' ^fig — ' 4 - 4 -’> 

where D(E ,e) = relative damage coefficient for a proton entering 
silicon solar cell with energy E q at an angle e 

D(E ,0) = relative damage coefficient for a proton of normal 
incidence (e = 0) with range equal to R(E ) cos e 
or energy E n 

N, .(E ) = the total number of silicon displacements created 
u " by a proton entering the silicon with energy E q 

cos e = the projected area of a unit cell area 

E n = R -1 [R(E ) cos e] 

n »0 

When the range of a proton incident at angle e exceeds the product 
of the thickness of the cell and the secant of e, D(E 0 ,e) is calculated 
as follows: 


D(E 0 ,e) = D(E q ,0) 


(4.4.2) 


Equations (4.4.1) and (4.4,2) allow the evaluation of equation (4.1.2). 

This integration has been done by machine using the D(E Q ,0) values shown 

in Figure 3.18. Separate integrations were done for D(E o ,0) values based 

on I , V and P 
sc’ oc max 


Evaluation of equation 4.1.2 for cell thicknesses of 0.0254 cm 
(0.010 in.) and 0.0457 cm (0.018 in) has shown that, for practical pur- 
poses, the results can be considered independent of cell thickness. The 
results of these integrations for several coverslide thicknesses are 
shown in Figures 4.3 and 4.4. The same data are printed in tabular form 
in Tables 4,3 and 4.4. 
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Figure 4.3 Relative Damage Coefficients for Space Proton Irradiation 
of Shielded N/P Silicon Solar Cells (Based on I ) 
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Figure 4.4 Relative Damage Coefficients for Space Proton Irradiation 

of Shielded N/P Silicon Solar Cells (Based on P m ,„ or 1/ 

' max oc 




'/STEMS 


PROTON OAMAGE COEFFICIENT FOR JSC 

OMNIDIRECTIONAL TO EOUIV 10 HEV PROTON 


SHIELD THICKNESS. GN/CM2 


3.35E-02 6.T1E-02 1.12E-01 


3.35E-01 



3 .046-13 

1.G4E+00 

0. 

4. L6E-13 

1. 02E+00 

3.25E-05 

4 .4 8E- 1 3 

9.95E-01 

4.01E-02 

4.80E-1 3 

9.63E-01 

1.79E-01 

5.12E-13 

9.28E-01 

3.46E-01 

5.44E-13 

8.93E-01 

4.80E-01 



5.76E-13 

8.59F-01 

5.78E-01 

0 . 

6.08E-13 

8.26E-01 

6.45E-01 

0. 

6.40E-13 

7.95E-01 

6.87E-01 

1 .25E-03 

6.72E-13 

7. 71E-01 

7.10E-01 

7.20E-02 

7.04E-13 

7.48E-01 

7.18E-01 

2.07E-01 

7.36E-13 

7.24E-01 

7.18E-01 

3.27E-01 

7.68E-13 
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Table 4.4. Proton Damage Coefficients for V Qc and P f 
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The values of relative damage constants for omnidirectional fluences 
of protons on shielded solar cells allow a space proton environment to be 
reduced to an equivalent fluence of normally incident 10 MeV protons on 
unshielded silicon solar cells. Experimental studies of silicon solar 
cells have indicated that a fluence of normally incident 10 MeV protons 
produces damage which can be approximated by a fluence of 1 MeV electrons 
which is 3000 times that of the 10 MeV proton fluence. 


The evaluation of the absorbed dose in shielding materials due to 
space protons requires an analysis similar to that done for space electrons. 
For this evaluation an expression similar to equation (4.3.2) is used. 

The quantity D(E q ) is replaced by the stopping power (- ^ ^-) for protons 
of energy (E q ) and the quantity D(E,t) becomes the absorbed dose per unit 
incident omnidirectional -flux protons of energy E and at shielding depth t. 
The result of this integration for several shielding thicknesses of fused 
quartz are shown in Figure 4.5 and Table 4.5. Rosenzweig has published 
similar data. 4 *® 

4.5 Alpha Particle Space Radiation Effects 

Solar flares have been shown to have a component of energetic alpha 
particles (helium nuclei). The evaluation of the effects of solar flare 
events on solar arrays requires alpha particle data similar to that for 
electrons and protons. Smith and Blue compared effects of 10.5 MeV pro- 
tons and 42 MeV alpha particles on silicon solar cell degradation. 4 * 7 
The results showed that the 42 MeV alpha particle flux degraded the 
silicon cells 3.8 times as fast as a similar flux of 10.5 MeV protons. 

These results were in good agreement with a theoretical damage ratio of 4. 

Based on the experimental results of Smith and Blue, the proton 
damage constant curve shown in Figure 4.4 can be translated a factor of 
four higher in energy and a factor of four higher in relative damage 
constant to represent a similar family of relative damage constants for 
alpha particles in space. Although the relationship found by Smith and 
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Figure 4.5 Absorbed Dose Per Unit Fluence of Space Proton for 
Various Depths in Planar Fused Silica Shielding 
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Table 4.5. Proton Stopping Power, Rad/Unit Omnidirectional Flux 
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Blue may not extend to lower particle energies, a set of effective damage 
constants for alpha particles, obtained by the above two translations, 
is shown in Figure 4.6. Data are shown based on P and V . 

II IQ A ULi 

Data based on I may be obtained similarly for alpha particle effects. 

4.6 Alternative Approaches 

The previously discussed methods of analyzing the effects of space 
radiation on shielded silicon solar cells can be summarized as follows. 

An effective damage constant was developed which represents the relative 
effect of an omnidirectional (isotropic) fluence of monoenergetic particles 
on a shielded solar cell in terms of an equivalent unidirectional irradia- 
tion by 1 MeV electrons or 10 MeV protons. The concept involves energy 
degradation of particles by the shielding and the effects of the energy 
and angle of incidence of particles entering the silicon solar cell. 

4 9 

An alternate approach to the problem has been used by Carosella ‘ 

4 10 

and Piccianno and Reitman. ' This method involves the calculation of 
the spectrum of particle energies exiting the cover slide (i.e., entering 
the silicon) of a solar cell with infinite back shielding. The calculation 
is extremely laborious and must be repeated for every environment con- 
sidered. The calculation of this "modified" energy spectrum destroys all 
angular information in the data and yields a spectruri which is neither 
monodirectional , omnidirectional or isotropic. The problems relating 
to the angular content of the "modified" spectrum emerging from the 
shielding are of no consequence in the calculation of absorbed dose, 
electron damage, or high energy proton damage. It has been shown that 
these quantities are independent of the angle of incidence. It is there- 
fore justified to weight the "modified" spectrum with relative damage 
coefficients to evaluate electron damage in terms of a damage equivalent 
monoenergetic normal incidence fluence. In the case of low energy pro- 
tons, the use of the referenced methods incorrectly assumes that proton 
damage is independent of the angle of incidence. This shortcoming is 
particularly serious in the case of many common space environments in 
which the lower energy proton damage dominates the solar cell degradation. 
The procedures used in this chapter are based on techniques previously 
described by Tada.^'^’ 
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CHAPTER 5 


5.0 THE SPACE RADIATION ENVIRONMENT 

The radiation environment near the earth, detrimental to solar 
cell performance, consists of electrons and protons trapped in the geo- 
magnetic field, corpuscular radiation associated with large solar flare 
activity, and to a lesser extent, galactic cosmic-ray radiation, Near 
Jupiter, an environment similar to the earth's trapped particle radiation 
is expected, but the intensity could be far greater than that near earth, 
due primarily to the anticipated large magnetic field. In the following 
sections, each environemt is qualitatively described to assist the reader 
in determining the proper environment for use on making solar cell degra- 
dation estimates. Quantitative, or detailed, descriptions of each envi- 
ronment are beyond the scope of this manuscript. 

5.1 Geomagnetical ly Trapped Radiation 

The geomagnetic dipole field is responsible for the existence of a 
radiation belt near the earth, holding the trapped charged particles for 
long periods of time. It is a plasma confined in an inhomogeneous mag- 
netic field, and the understanding of transport, loss and capture mechan- 
isms has improved considerably over recent years. The dynamics of this 
radiation environment are greatly influenced by solar activity. 

Geomagnetical ly trapped radiation may be either of natural origin 
or of artificial origin, such as high-altitude nuclear explosions. A 
particle has to possess a charge to be trapped in a geomagnetic field, 
and the constituents are electrons and protons. Regardless of the origin, 
the particle with just the right momentum and pitch angle will be trapped 
in the field. The particles will then spiral about a field line with 
varying pitch angle and curvature in the inhomogeneous field. They con- 
tinue the motion until they reach the mirror (or reflection) point where 
the pitch angle becomes zero, and then bounce back into the other hemi- 
sphere. They continue to bounce back and forth between the mirror points 
(latitudinal motion), and at the same time drift in the longitudinal 
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direction as the result of forces due to the gradient of field strength 
and the curvature of field lines. During a quiescent state (periods of 
normal solar activity), therefore, the trapped particles can be charac- 
terized by three periodic motions: (a) cyclotron oscillation about the 

field line with Larmor frequency, (b) latitudinal motion between mirror 
points, and (c) longitudinal drift. The direction of motion for electrons 
is obviously opposite to that of protons because of an opposite charge. 
Near the mirror points, the particles collide with upper atmospheric 
gases, gradually losing their energy and changing trajectory until they 
are lost in the lower atmosphere. 

At some distance from the earth, the field is distorted by the "solar 

wind" as shown in Figure 5.1. It is a plasma frcm the sun, consisting 

mostly of protons with average energy of a few keV and of density on the 
3 

order of 10/cm . The solar wind interacts with the geomagnetic field 
resulting in the formation of a shock wave. As the solar plasma passes 
the shock wave, the random speeds of the particles increase, and the 
magnetic field becomes 'turbulent. There is a region of hot plasma near 
the earth-sun line on the day side. The solar wind deforms the geomag- 
netic field and squeezes it into a cavity called the magnetosphere. The 
magnetosheath ranges from 10 to 14 earth radii. 

The geomagnetic field lines just behind the magnetosheath are quali- 
tatively similar to those associated with the simple dipole model and 
trap corpuscular radiation as described above. During quiescence, a 
relatively steady flow of solar wind blows the field away from the sun, 
contributing to an asymmetric shape of the radiation belt, compressed on 
the sun's side and forming the tail of the magnetosphere and the thin 
neutral layer on the dark side of the earth. 

5 1 

Mcllwain ’ has proposed a coordinate system consisting of the mag- 
netic field B and the integral invariant I which can adequately relate 
measurements made at different geographic locations. He introduced a 
parameter L = f (B , I ) , analogous to a physical distance in a dipole field 
(the equatorial radius of a magnetic shell), thus reducing the number of 
variables needed to describe the physical situation of trapped charged 
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particles and presenting field data in a manner which facilitates its 
physical interpretation. For a radial distance of R and a dipole moment 
of M, the transformation using the dipole relation is expressed as 
follows: 


R 3 L 

where R is L cos 2 x, M is the geomagnetic dipole moment, and \ is the 
magnetic latitude, Mcllwain expanded the parameter L into a polynomial 
function of a variable which is a function of I, B, and M and elegantly 
represented the physical phenomena of trapped particles. 

Since its introduction, numerous particle field data were presented 
in this (B,L) coordinate system. Vette and coworkers have concentrated 
efforts on the compilation of particle field data reported by numerous 
investigators and have constructed models of the radiation environment. 

These data are regarded as the best consolidated source of information 
available on trapped radiation environments, and are used as the single 
source of data on this subject throughout this manuscript. The reader 

c pc o 

may consult the referenced publications for detailed and quantitative 

discussions of the trapped electron and proton environment models. 

5.1.1 Trapped Protons 

The most recent description of the trapped proton environment is 
presented in references 5.2 through 5.4. The largest proton concentra- 
tion of intermediate energies is near the earth within an L-value of 
four (geocentric) earth radii, peaked at about two earth radii. The 
high energy protons concentrate even closer to the earth, peaked at 
1.5 earth radii whereas the distribution of the lower energy protons 
extends nearly to synchronous altitude (L = 6.6 R g ). Generally speaking, 
the energy spectrum becomes softer as the L-value increases. At synchro- 
nous altitude, the spectrum is so soft that practically no protons with 
energy greater than two MeV exist. 
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5.1.2 Trapped Electrons 


Trapped electrons with energies of a few hundred keV extend to the 
outer boundary of the magnetosphere, which fluctuates at 8 to 10 earth 
radii. There are two intense regions: an inner one covers the L-values 

in the range of 1.2 < L < 2.8 and peaks about 1.4 earth radii, whereas 
the outer zone ranges 3 < L < 11 and peaks at around 4 to 5 earth radii 
with the flux about 10 7 eler.trons/cm 2 -sec for both zones. 

The outer zone is a very dynamic region of space, and the particles 
are considered to be pseudo-trapped because the lifetimes are shorter 
than the drift time around the earth. However, powerful sources supply 
electrons to this region of space, and thus substantial fluxes are al- 
ways present. In this zone, the flux has large short-term temporal varia- 
tions related to the local time as well as a long-term change in the 
average flux associated with a solar cycle. 

In the inner zone, the effect of storms on the average flux is 
significant at high L-values and higher energies. A long-term increase 
in the inner zone flux is correlated with the increased solar activity. 
Another source of temporal variation is due to a decay of residual elec- 
trons from the Starfish nuclear explosion. These temporal variations 
are accommodated in the recent compilation of data and publications on AE3, 
AE4 and AE5 by Vette. 5 ’ 5 ' 5 * 8 

5.2 Orbital Integration 
5.2.1 Circular Orbits 

Vette and coworkers have time integrated both the trapped proton 
and electron environments for convenient energy ranges and tabulated 
the average daily fluence for various altitudes and inclinations. There 
are two forms of spectra in his data: one is of the form of integral 

flux and another difference flux, the latter of which should not be 
confused with the differential flux. 

If <f> CE ) is a differential flux at energy E in MeV, normally expressed 
in terms of particles/cm 2 -sec-MeV, and $(>E) is an integral flux with an 
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energy greater than E. expressed in particle?/cm 2 -sec, tne relationship 
of these two quantities is 



= E 4>(E.) A E. 

■i J J 

U 

On the other hand, the difference fluence a$ is simply 

A$. = $ (>E . ) - $(>E . + aE) (5.2.2) 

J J J 

5.2.2 Trajectories Other Than Circular Orbits 

For the spacecraft trajectories other than circular orbits tabulated 
in references 5.2 through 5.7, the radiation environment encountered by 
the spacecraft may be determined by some other method. 

One approximation for hand calculations is to divide the trajectory 
into small segments with suitable time intervals At(r,i), so that during 
this flight time interval the environment can be regarded the same as that 
of a circular orbit at altitude r and inclination i. The environment 
weighted by flight time then becomes 

*(>E,I) = E *'(>E,r,i) • At(r,i) (5.2.3) 

R 

The difficulty of this approximation is that (a) the $' is averaged over 
a circular orbit so that the o' is not equal to an instantaneous flux at 
(r,i), and (b) the i is constantly changing if i f 0. Thus an error can 
be very large; in fact, it will be appropriate only for an order of mag- 
nitude calculation if i f 0. 

A more accurate and simpler way of determining the environment is 
to determine a trajectory on an isoflux contour map of energy E p plotted 
on geographic coordinates. By knowing instantaneous flux at (r,i) and 
the time interval At(r,i), the integral flux can be time integrated by 
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(5,2.4) 


*(> E n ) 



At ( r . i ) 


This calculation will determine the flux of an integral spectrum at E . 

If isoflux contour maps for different energies are available, a series 
of such calculations leads to several points on an integral spectrum. 

The most sophisticated way to determine the environment is to make 
use of the physically significant coordinate system (B,L) so that uncer- 
tainties and inaccuracies attributable to the geographic coordinate 
system are eliminated. A set of state vectors or classical orbital ele- 
ments can be used to solve Kepler's equation and generate a trajectory 
with suitable time intervals. These geographic coordinates are then 
transformed into geomagnetic shell coordinates (B,L) on which isoflux 
contour maps are plotted. This approach is computationally involved 
and hence is practical only with the aid of a computer. 

The instantaneous flux 4> ( >E ,t.) is thus determined and is time 

* J 

integrated on each flux map of energy E p in the following manner: 

<f(>E p ) = Z ( > E „ » "t j ) • Atj (5.2.5) 

J 

Upon time integration of instantaneous flux throughout a given trajectory, 
performed on one isoflux contour map of specified energy E n and particle 
type, one point is finally determined in an integral flux-energy spectrum 
$(>E). If similar calculations are performed on a number of maps of 
different energies, exactly the same number of points can be determined 
in the final spectrum. 

An energy spectrum at an arbitrary point in space, in general, is 
a function of both B and L coordinates and can be expressed in either 
exponential or power form. If such a distribution function in either 
form is applicable to an entire energy region for all points in space, 
only one isoflux contour map is required to determine the time integrated 
flux-energy spectrum. 
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If FN(>E ,>E,,Lt) is an energy distribution function for a parti- 
n k j 

tioned energy greater than E k 

+(>E k ,t j ) = 4>(>E n ,t j ) • FN(>E n ,>E k ,Lj) (5.2.6) 

The L'. in FN is the L value in the input table and is the nearest to the 
J 

computed value L^. The E n is an energy specified in an isoflux contour 
map. The flux <f>(>E k ,tj) is integrated in each energy shell E k as the 
time integration proceeds throughout a trajectory. 

The distribution function for an exponential spectrum with a param- 
eter E o (>E n ,B,L) is defined as 

FN(>E n ,>E k ,t j ) = exp [(E n - E^/E^E^ ,L . )] (5.2.7) 

" ex > C( E n ' Wh&n'W 


The FN is thus normalized at E„, and hence the E„ should agree with the 

n n 

energy specifying an isoflux contour map. For a power form, the distribu- 
tion function is defined as 

-P(>E 

™< >E n* >E k‘ t j> = < E k /E n> J J 


where P(>En, Bj L j ) = Exponent of power form which depends on the energy 

specifying the isoflux contour map, as well as 
special location of trajectory in B-L coordinate 
(Bj.Lj) 

and again the FN is normalized at E . As an example of the above method, 
the geomagnetical ly-trapped proton environment is machine calculated for 
a highly elliptical orbit with an inclination of 63.5 degrees and is shown 
in Figure 5.2. 
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INTEGRAL PROTON FLUX (>E), p/cm 







5.3 Cosmic-Ray (Galactic Cosmic-Ray) Radiation 

Galactic cosmic rays are a highly penetrating radiation originating 
beyond the solar system. They possess energies greater than 1 BeV 
(some may exceed 10 8 BeV) and are capable of extraordinary interactions 
with matter in the upper atmosphere such as spallation, fission, frag- 
mentation, and the subsequent secondary processes. The local cosmic-ray 
radiation in the atmosphere contains protons, neutrons, pi-mesons, mu- 
mesons, electrons, photons, and strange particles. 

Near the upper limits of the atmosphere, the primary radiation, 
consisting of 79 percent protons and 20 percent alpha particles, pre- 
dominates over the products of nuclear reactions and the decay products, 
thus the components change with altitude. 

The ability of charged particles to penetrate a magnetic field is 
limited by the Lorentz force and is measured by a quantity called the mag- 
netic rigidity, defined by the ratio of the momentum to the charge. The 
radius of curvature of charged particles in the field is then related 
to the magnetic rigidity, and hence the ability of particle penetration. 
The magnetic cutoff momentum, and hence the cutoff energy, for a given 
vertically incident particle at a given altitude is closely related to 
the latitude of the geomagnetic field. Only protons with energies 
greater than about 15 BeV can penetrate tne earth's magnetic field at 
the equator. 

One remarkable characteristic of cosmic rays is their isotropy. 

The average diurnal effect is very small; however, there is a definite 
relationship between the fluctuation and solar activity in general; 

27-day effects, an 11 -year fluctuation cycle, and the Forbush decrease 
associated with the magnetic storms are examples. Although the energy 
is very high, the flux is negligibly small compared with other environ- 
ments considered, and this environment is ignored in solar cell array 
degradation cases at present. 

5.4 Solar Flare (Solar Cosmic-Ray) Radiation 

Solar flares occur in the neighborhood of sunspots, very seldom emit 
white light, and cause a sudden increase in intensity of the hydrogen 
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alpha line (6,563 A). After its inception, the flare rapidly expands 
over an area of a few million to a billion square miles of the solar 
disk, reaching a peak intensity and gradually decaying and completely 
disappearing within several minutes to several hours, depending on the 
size of the flare. 

Within half an hour or more following the appearance of large solar 
flares, energetic particles, consisting mostly of protons, are detected 
at the earth, particularly in the polar regions inside the auroral zones. 

The radiation dies away with a time constant of one to three days. The 
constituent particles are electrons, protons, alpha particles, and very 
small numbers of medium nuclei (C, N, and 0). The ratios of protons to 
alpha particles, and of protons to medium nuclei vary considerably between 
solar events, whereas the ratio of alpha particles to medium nuclei remains 
relatively constant. 

Although the fluctuation in flux intensity is much more severe and 
random than those of galactic cosmic rays, the following phenomena were 
observed in the past: (a) there may be an 11 -month cycle in the peak 

number of events, (b) there is a semiannual variation which has maxima 
in March and September, probably near the equinoxes, (c) the maximum 
number of events occurs on the average near the September equinox and 
the minimum during December or January, (d) the number of flares varies 
with the 11 -year solar cycle, and (e) there is a definite tendency for 
flare events producing a large proton fluence to occur during the increase 
or decrease of sunspot activity rather than during the maximum. Predicted 
and observed sunspot numbers for the previous solar cycle are shown in 
Figure 5.3. 5,10 

Solar flare particle fluxes arriving at the earth are highly time 
dependent in intensity, spectrum and isotropy. The rise time naturally 
varies with individual event and is strongly energy-dependent, reaching 
the maximum intensity first at higher energies and thus showing a harder 
spectrum at the beginning. After the peak of radiation, the integral 
flux decays with time at a rate approximately proportional to t n , where 
t is time and n is a number, roughly equal to 3. The particle flux arriv- 
ing in the upper atmosphere is for the most part isotropic; however, 
significant anisotropies frequently exist for shorter durations, arriving 
from a highly preferred and fairly narrow direction in space from 30° to 
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60° west of the earth-sun line for a period of a few minutes. 

5 11 

A model described by McCracken ‘ is largely based on experimental 
observation. The following is an excerpt from his article: 

* 

'Figure 5.4 shows the model for the magnetic regime 
created by the plasma disturbance originating in a large 
solar flare (not necessarily a flare that results in a 
large proton event). The plasma ejected by the flare 
carries the lines of force of the sunspot with it, the 
lines of force being stretched outward from the sunspot 
in a quasi -radial fashion. The sun's rotation causes 
the lines to curve westward. The configuration of the 
lines of force near the leading edge of the plasma dis- 
turbance is not yet known; however, it is known to exclude 
cosmic rays from outside, and to inhibit the escape of 
cosmic rays injected at points inside the magnetic regime. 

In the Forbush decrease, the arrival of the leading edge of 
the plasma disturbance at the earth initiates the magnetic 
storm, and once the earth is inside the magnetic regime, 
some galactic cosmic rays are screened away from the earth. 

This phenomenon now provides a direct magnetic connection 
from the earth to the sunspot group. Consequently, if 
another flare were now to produce cosmic rays, they would 
travel rapidly along the magnetic lines of force to the 
earth. They would, therefore, arrive at the earth soon 
after the occurrence of the flare (about 20 minutes), and the 
maximum intensity would rapidly be reached. The divergent 
nature of the magnetic lines of force implies that the 
cosmic rays would tend to become collimated, eventually 
travelling roughly parallel to the lines of force. Also, 
the particles would be partially trapped within the mag- 
netic configuration; and so after a period of anisotropy, 
a period of isotropy may be observed. 

★ 

At a point outside the magnetic regime, Figure 5.4a .there 
is no direct connection to the sunspot, and hence cosmic 
rays cannot arrive rapidly at the earth. They can only 
arrive by diffusion across the lines of force-a process 
that tends to delay and isotropize them. Therefore, an 
appreciable time delay exists between particle production 
and arrival at the earth (30-120 minutes), and the intensity 
rises slowly to a maximum some hours after the flare. The 
maximum omnidirectional intensity of radiation is less than 
that which would be observed if the earth has a direct mag- 
netic connection to the sunspot. The radiation may be 
mildly anisotropic, with the maximum intensity oriented 
along the lines of force leading to the sun, but not to 
the sunspot group in which the flare occurred." 


★ 

Figure numbers altered from original text. 
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Figure 5.4 Changes of the Interplanetary Magnetic 
Field Regime Model with Time^'^ 
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Until recently, observations of solar flare particles were made 
only for those with relatively high energies (10, 30 and 100 MeV), 
and much higher than the energy range of normal interest in connection 
with solar cell degradation. The time integrated spectrum normally 
exhibits an exponential form with respect to rigidity and is customarily 
expressed in terms of the characteristic rigidity R q as follows: 

O-R/RJ 

*(>R) = $(>R 0 ) • e 0 (5.4.1) 

where R = rigidity / J°^ e \ . /volt-amp-secX or (volt) 

\coulomb/ \ coulomb J 

= pc/zq = yj E 2 - (m Q C 2 ) 2 /zq 

E = total energy 

T = kinetic energy 

p = momentum ( Me V/c),( joule sec/ m), or (newton. sec) 
m Q c 2 = rest mass energy, 938 MeV per proton 
zq = atomic charge 

<2>(>R) = integral flux having rigidity 
greater than R 

The R q varies not only with each event but within the spectrum of 
an event. Integral solar proton flux is tabulated in Table 5.1 for 
selected flare events from the year of 1956 through 1962 together with 
the characteristic rigidity. The characteristic rigidity R Q ( 1 0-30 ) is 
computed from the flux at 10 and 30 MeV and R Q ( 30-1 00 ) at 30 and 100 
MeV respectively. The annual integral flux is shown in Table 5.2. The 
R q computed for the annual flux is smaller during the years near sunspot 
maximum (50 ^ 70 MV) but the total annual fluence is higher during these 
years. 

Since solar flare particle fluxes are rich in low rigidities, a 
strong cutoff phenomenon is expected. During the quiescent state, the 
cutoff rigidity at low latitude is a strong function of direction as 


/T(T + 2m o c 2 )/zq 
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TABLE 5.1 INTEGRAL PROTON FLUX AT 10, 30, AND 100 MeV AND 

CORRESPONDING CHARACTERISTIC RIGIDITY R ^ 5 ‘ 12 ^ 

o 


DATE 

(> 10 MeV) 

$(> 

30 MeV) 



(protons/can 2 ) 

2/23/56 

1.8 x 10 9 

1.0 

x 10 9 

3/11/56 

- 


- 

8/ 3/56 

- 

2.5 

x 10 7 

11/13/56 

- 


- 

1/20/57 

- 

2 

CO 

O 

X 

4/ 3/57 

- 


- 

6/22/57 

- 


- 

7/ 3/57 

- 

2 

x 10 7 

8/ 9/57 

- 

1.5 

x 10 6 

8/29/57 

- 

1.2 

x 10® 

9/21/57 

- 

1.5 

x 10® 

10/20/57 

- 

5 

x 10 7 

11/ 4/57 

- 

9 

x 10 6 

2/ 9/58 

- 

1 

x 10 7 

3/23/58 

2 x 10 9 

2.5 

x 10® 

4/10/58 

- 

5 

x 10® 

7/ 7/58 

1.8 X 10 9 

2.5 

CO 

O 

X 

8/16/58 

4 x 10® 

4 

x 10 7 

8/22/58 

8 x 10® 

7 

x 10 7 

8/26/58 

1.5 x 10 9 

1.1 

x 10® 

9/22/58 

9 x 10 7 

6 

x 10® 

5/10/59 

5.5 x 10 9 

9.6 

x 10® 

6/13/59 

- 

6.5 

x 10 7 

7/10/59 

4.5 x 10 9 

1.0 

x 10 9 

7/14/59 

7.5 x 10 9 

1.3 

x 10 9 

7/16/59 

3.3 x 10 9 

9.1 

x 10® 

8/18/59 

- 

1.8 

x 10 6 


>(> 100 MeV) 

R Q (30-100) 

R o (10-30) 

3.5 x 10® 

(MV) 

195 

173 

- 

- 

- 

6 x 10® 

144 

- 

- 

- 

- 

7 x 10 6 

61 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

3 x 10 6 

56 

- 

- 

- 

- 

1 x 10 7 

127 

- 

- 

- 

- 

_ 

_ 

- 

1 x 10 7 

64 

49 

- 

- 

- 

9 x 10® 

62 

52 

1.6 x 10® 

64 

44 

1.8 x 10® 

56 

42 

2.0 x 10® 

51 

39 

1 x 10 s 

50 

38 

8.5 x 10 7 

84 

58 

- 

- 

- 

1.4 x 10® 

104 

68 

1.0 x 10® 

80 

58 

1.3 x 10® 

105 

79 
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TABLE 5.1 (Continued) 


DATE 

(> 10 MeV) 

$(> 

30 MeV) 

«{> 100 MeV) 

R Q (30-100) 

R o (10-30) 





(protons/cm 2 ) 



(MV) 


1/11/60 




4 

X 

10 5 

- 


- 

- 

4/ 1/60 

1.5 

X 

10 7 

5.0 

X 

10 6 

8.5 x 

10 s 

116 

93 

4/ 5/60 


- 


1.1 

X 

10 6 

- 


- 

- 

4/28/60 

1.3 

X 

10 7 

5.0 

X 

10 6 

7 x 

10 s 

104 

102 

4/29/60 


- 


7 

X 

10 6 

- 


- 

- 

5/ 4/60 

1.2 

X 

10 7 

6 

X 

10 6 

1.2 x 

10 6 

127 

147 

5/ 6/60 


- 


4 

X 

10 6 

- 


- 

- 

5/13/60 

1.5 

X 

10 7 

4 

X 

10 6 

4.5 x 

10 5 

94 

77 

6/ 1/60 


- 


4 

X 

10 5 

- 


- 

- 

8/12/60 


- 


6 

X 

10 s 

- 


- 

- 

9/ 3/60 

9 

X 

10 7 

3.5 

X 

10 7 

7 x 

10 6 

127 

108 

9/26/60 

2 

X 

10 7 

2.0 

X 

10 6 

1.2 x 

10 5 

73 

44 

11/12/60 

4 

X 

10 9 

1.3 

X 

10 9 

2.5 x 

10 8 

124 

91 

11/15/60 

2.5 

X 

10 9 

7.2 

X 

10 8 

1.2 x 

10 8 

114 

82 

11/20/60 

1.4 

X 

10 8 

4.5 

X 

10 7 

8 x 

10 6 

118 

90 

7/11/61 

1.7 

X 

10 7 

3 

X 

10 6 

2.4 x 

10 s 

81 

59 

7/12/61 

5 

X 

10 8 

4 

X 

10 7 

1 X 

10 6 

56 

40 

7/18/61 

1 

X 

10 9 

3 

X 

10 8 

4 x 

10 7 

102 

85 

7/20/61 

1.5 

X 

10 7 

5 

X 

10 6 

9 x 

10 5 

120 

93 

9/28/61 

5 

X 

10 7 

6 

X 

10 6 

1.1 X 

10 6 

121 

48 

11/10/61 


- 



- 


- 


- 

- 

2/ 4/62 


- 



- 


- 


- 

- 

10/23/62 

6 

X 

10 5 

1.2 

X 

10 s 

1 X 

10 1 * 

83 

63 
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TABLE 5.2. OBSERVED ANNUAL INTEGRAL SOLAR PROTON FLUX 


Year 

Number of Events 

Inteq 

2 

ral Flux (protons/cm ) 

R o 

$(>10 MeV) 

$(>30 MeV) 

$(>100 MeV) 






(MV) 

1956 

4 

2.0 x 10 9 

1.0 x 10 9 

3.5 x 10 8 

250 

1957 

9 

— 

4.0 x 10 8 

2.0 x 10 7 

60 

1958 

8 

7.0 x 10 9 

7.8 x 10 8 

2.4 x 10 7 

50 

1959 

6 

2,2 x 10 10 

4.2 x 10 9 

4.6 x 10 8 

70 

1960 

15 

6.8 x 10 9 

2.2 x 10 9 

3.8 x 10 8 

130 

1961 

6 

1.6 x 10 9 

3.5 x 10 8 

4.2 x 10 7 

90 

1962 

2 





1963 

1 






52 






Solar-cycle total 

3.9 x 10 10 

9.0 x 10 9 

1.3 x 10 9 

100 (av.) 
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well as of latitude (approximately proportional to cos 4 X for large 
geomagnetic latitudes), and hence of L. Galactic cosmic rays follow 
this normal Sttirmer cutoff as do the flare particles just before the 
plasma cloud hits the geomagnetic field. After the impact of the plasma 
front, the field is disturbed (magnetic storm) in such a manner that 
the field due to a time-dependent ring current appears to superimpose 
on the normal geomagnetic dipole field, causing the disturbed line of 
force to stretch further out of the earth at a given latitude. As a 
result, the particle rigidity necessary to penetrate at a given latitude 
is greatly reduced, and the cutoff energy becomes time dependent. A 
recent satellite observation indicated that the cutoff energy at syn- 
chronous altitude seems to be much less than that many expected, and that 
flare protons with energy as low as a few hundred keV were observed 
during the storm. If this is the case, the cutoff energy due to the geo- 
magnetic field becomes insignificant at this altitude, because the cut- 
off due to a solar cell cover shield is normally far greater than the 
magnetic cutoff during a storm. If both altitude and latitude are low, 
the field perturbation due to the storm may be insignificantly small 
compared with that of the quiescent state, and the Stormer cutoff approxi- 
mation may prevail. The geomagnetic shielding phenomena are shown in Figure 
5 13 

5.5 ’ for protons in a class three flare on July 18, 1961. 

For the purpose of predicting the size and spectrum of solar flare 
proton events, many statistical analyses have been made on proton events 
observed near or on the earth. Unfortunately, the correlation between 
the prediction and observations has been rather poor. A Poisson distri- 
bution may be appropriate for sunspot numbers and solar flares on the 
sun, but not for solar flare proton events. The flares which are large 
enough to emit a large number of energetic particles and further satisfy 
the requirements of protons to reach the earth obviously belong to a 
special class of solar flare events. Phenomena observed during 
solar cycle 19 are enumerated below for review, placing particular em- 
phasis on those which appear to be dependent on solar activity. 
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Figure 5.5 Solar Flare Proton Environment at 200 n. mi. Circular 
Due to Flare Event on July 18, 1961, Class Three Flare 



a. The flares capable of producing large proton events tend 
to occur when the rate of change in annual sunspot 
number becomes greater. 

b. The characteristic rigidity of solar flare protons is 
randomly distributed throughout an 11 -year cycle, 
but both the annual expectation value and variance 
are not. During a period of increasing or decreasing 
sunspot activity, the R 0 becomes larger on the average 
than that during the maximum, and the variance becomes 
smaller during the solar maximum. That is to say, the 
solar flare proton events are relatively steady and 
confined in a smaller rigidity range during the solar 
maximum, whereas the size and spectrum become erratic 
when the rate of change in sunspot activity becomes 
severe. 

c. The size of each event, as measured by an integral 
proton flux of energy greater than 30 MeV, is 
almost randomly distributed over an 11 -year cycle, 
but a line connecting the successive annual fluence 
plotted against sunspot number is not a single- 
valued function. 

A model was developed to implement the above observations. With 

the aid of the model, annual free space solar flare proton fluxes of 

energy greater than 10 and 30 MeV are determined for each year of cycles 

20 and 21, and are tabulated in Tables 5.3 and 5.4, respectively. These 

predictions are compared with those described in reference 5.13 as well 

5 15 

as with some experimental observations compiled by King. ' These pre- 
dicted flux values agree with the experimental values very well for most 
of the years of cycle 20 except for the years of 1966 and 1970. The dis 
crepancy between the predicted and the observed $>( >10 MeV) can be as 
large as a factor of 7, as in the case of 1966 and 1970, when the annual 
flux variation is considered. However, the overall discrepancy between 
observed and predicted 4> is within a factor of 2 when the flux is inte- 
grated over the 11 -year solar cycle. The August 1972 event is a severe 
fluctuation in magnitude and time of occurrence with respect to the 
11 -year cycle; yet the predicted annual value is within a factor of 4 of 
the observed value. Despite the crudeness of this technique and a disap 
pointing prediction for some years of cycle 20, this technique will 
provide a much better estimate than those cited in reference 5.13, 2n 
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TABLE 5.3 BOTH PREDICTED AND OBSERVED ANNUALLY INTEGRATED SOLAR FLARE PROTON FLUXES FOR CYCLE 20 
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$(>10) Integral proton flux of energy greater than 10 MeV, protons/cm 2 -yr 

R Q ( 1 0-30) Characteristic magnetic rigidity calculated at energies of 10 and 30 MeV, MV 

CF (> 10 ) Correlation factor (CF) defined as: 

CF {>10) = $(>10) (0bserved)/$(>10) (Predicted) 



TABLE 5.4. PREDICTED ANNUALLY INTEGRATED SOLAR FLARE PROTON FLUXES FOR CYCLE 21 
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TOTAL 2.6E11 5.5E10 4.9E9 8.5E8 





















Figure 5.6 


Predicted Smoothed Sunspot Number for Solar Cycle 
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improvement of better than a factor of 10. The environment specified 
in reference 5.13 is shown in Figure 5.7, together with data taken from 
references 5.11 and 5.12 for comparison. 

5.5 Jupiter Trapped Radiation Environment 

Earth observations of decimetric RF noise emanating from the vicinity 
of Jupiter suggest the following as a Jovian environment: 

a. Radio emission generated near Jupiter is synchrotron 
radiation. 

b. Charged particles are responsible for the RF noise and 
are trapped by a magnetic field as in the case of near 
earth environment. 

c. The magnetic field is probably due to a dipole in large 
part and is very strong as compared with that of earth 
(probably % 10‘ 5 of telsa at the Jovian equator). 

d. Consequently, the intensity of the Jovian radiation 
belt is far greater than that of the earth's radiation 
belt which can be scaled upward accordingly to satisfy 
the observed synchrotron radiation. 

Prospective models for the Jovian radiation belt were postulated and 
discussed at the Conference on Jupiter Radiation Belt Workshop held at 
JPL, Pasadena, California, July 1971. A few well-known experts on the 
earth's radiation belt raised questions about some postulates and assump- 
tions made in the proposed model. These questionable aspects will become 
clarified in the not too distant future as space probes fly over Jupiter. 
The first in situ data will probably be available from the Pioneer 10 
mission in December 1973. The post-workshop radiation models are shown 
in Figure 5.8. 
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INTEGRAL SOLAR FLARE PROTON FLUX (>E), PROTON/cm 2 - YEAR 



Figure 5.7 Solar Flare Proton Environment of the Last Solar Cycle, 19 
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30 


Figure 5. 


2 3 4 56789 10 

MAGNETIC SHELL PARAMETER, JUPITER RADII 


8 Fluxes of charged particles in Jupiter's trapped radiation 
belts, as functions of distance from the magnetic dipole 
in the magnetic equatorial plane. Local values of the 
characteristic energy E 0 are shown in MeV.5-16 
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CHAPTER 6 


6.0 SOLAR ARRAY DEGRADATION CALCULATIONS 

In the previous sections, the three basic input elements necessary 
to perform degradation calculations were developed. The first of these 
elements is degradation data for solar cells under normal incidence 1 MeV 
electron irradiation. The second input element is the effective relative 
damage coefficients for omnidirectional space electrons and protons of 
various energies for solar cells with various thickness of cover glasses. 
The third input element is space radiation environment data which has 
been integrated for the orbit of interest. The section will cover the 
use of the above inputs to complete a solar array degradation estimate. 



The effective relative damage coefficients allow the conversion of 
various energy space electrons and protons into equivalent fluences. 

The equivalent fluences are based on normal -incidence monoenergetic 
irradiations for which the degradations of the solar cells of interest 
are characterized. The process of weighting an integral energy spectrum 
of electrons for a given orbit can be described as follows: 


*1 MeV e = ^2 [ * (>E) ' * (>E + AE)] ‘ D(E ’ t} (6J J } 
E=0 

$ 1 MeV e = the damage equivalent 1 MeV electron 
fluence (e/cm 2 -year) 

$(>E) - <J>(>E + aE) = the isotropic particle fluence having 

energies in a small energy increment 
greater than energy E (e/cm 2 -year) 

D(E,t) = the relative damage coefficient for 
isotropic fluences of space particles 
of energy E on solar cells shielded 
by cover glass of thickness t 
(dimensionless) 
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The quantities 4>(>E) - $(>E + AE) for a range of energies are also known 
as the difference spectrum. This spectrum can be generated from an 
integral energy spectrum for any energy increments desired. For the case of 
space protons, equation (6.1.1) can also be used with the exception 
that D(E,t) values for protons are based on 10 MeV proton fluences 
rather than 1 MeV electrons. The calculated equivalent fluence will 
therefore be a damage equivalent 10 MeV proton fluence. The equivalent 
10 MeV proton fluence can be converted to equivalent 1 MeV electron 
fluence as follows: 


$ 1 MeV e 


$ 10 MeV p 


x 3000 


( 6 . 1 . 2 ) 


The above relationship is an approximation which must be made for the 
purpose of combining electron and proton damage. In Section 3.3, the 
differences between electron and proton degradation were discussed. 

Since the slope of the degradation curve (the constant C in equation 
3.2.1) is different for 1 MeV electron and 10 MeV proton irradiations, 
the constant in equation (6.1.2) will differ depending on the level of 
degraded cell output at which this constant is determined. At present, 
the best information available indicates a value equal to .3000 when cell 
output parameters are degraded by 25%. In cases when the cell degrada- 
tion is entirely dominated by proton damage, such as circular earth orbits 
below 1000 km (550 n mi), the cell degradation could be estimated more 
accurately by calculating the equivalent 10 MeV proton fluence, and 
referring to 10 MeV proton cell damage curves which are not currently 
available. 


An additional problem arises in calculating equivalent fluences for 

proton environments. The results shown in Figures 4.3 and 4.4 have 

shown that different values of D(E,t) for proton irradiation are found 

when this damage constant is based on cell I or P and \l . This 

sc max oc 

differs with the results of electron irradiation where one value of 
D(E,t) describes the behavior of all cell output parameters. Because 
of the two sets of D(E,t) values for proton irradiation, two different 
equivalent 10 MeV proton fluences must be considered. One of these will 
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The other will 


describe the variation of solar cell P mav and V . 

Ilia X UL 

describe the variation of solar cell I S(; . 

The values of D(E,t) have been calculated by assuming infinite back 
shielding. Although this condition is often approached in the body- 
mounted solar arrays of spinning spacecraft, it is not generally true. 

The designer must also evaluate the contribution of equivalent fluence 
resulting from radiation incident on the back side of the solar cells. 

The result is a front and a back component of equivalent fluence. A 
question arises as to the values of D(E,t) to be used for back irradia- 
tions. In the case of trapped space electron irradiation, it is reason- 
able to use the same values of D(E,t) for both front and back irradiations. 
The only problem in this case is to reduce the backshielding offered by 
the panels, satellite, etc., to an equivalent in gm/cm 2 of planar shielding. 

The case for space protons is considerably more complex because of 
the nonpenetrating nature of low energy protons. There is a need for 
D(E,t) values calculated from rear irradiation experimental results. If 
such data were available, it is expected that the values of D(E,t) below 
10 MeV would be considerably smaller than those found for front irradia- 
tion with protons. This difference would be expected because rear inci- 
dent low-energy protons are shielded from the space charge region of the 
solar cell by the bulk region of the cell. Since the desired proton 
back irradiation D(E,t) is not available, the only alternative is to use 
the front irradiation data. To allow for the self-shielding effect for 

cells irradiated with protons from the rear, the thickness of the cell, 

o 

minus 0.011 gm/cm (0.002 in) should be included in the total back- 
shielding. 

The various contributions and variations of equivalent fluence which 
can be encountered in a natural space environment are summarized in Table 
6.1. Columns in the right side of the table indicate the contributions 
from the various radiation components to the two different types of equiv- 
alent fluence. Although the most general case can involve all the con- 
tributions shown in Table 6.1, in a typical earth orbit only a few of 
these contributions may be significant. 
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TABLE 6.1 

SUMMARY OF EQUIVALENT FLUENCE CONTRIBUTIONS 

Contributions 

I P V 

sc max* oc 


1. 

Trapped electrons, front, ( I sc • P max , V Qc ) 

X 

X 

2. 

Trapped electrons, back, (I sc> P max , V QC ) 

X 

X 

3. 

Trapped protons, front , U sc ) 

X 


4. 

Trapped protons, back, (I ) 

X 


5. 

Trapped protons, front, (P max » V QC ) 


X 

6. 

Trapped protons, back (P max , V Q(; ) 


X 

7. 

Flare protons, front, (I ) 

j L 

X 


8. 

Flare protons, back, (I S( J 

X 


9. 

Flare protons, front, (P max , V QC ) 


X 

10. 

Flare protons, back, (P max » V QC ) 


X 


Thus, 1, 2, 3, 4, 7, and 8 contribute to the I total equivalent fluence 
and 1, 2, 5, 6, 9, and 10 to the P max > V QC total equivalent fluence. 

The calculation of equivalent fluence and subsequent estimation of 
degraded solar cell output from the data in Figures 3.7 through 3.14 
yields data which are valid for temperatures of 28°C and solar illumina- 
tion power densities of 135 mW/cm 2 . When degraded solar cell outputs are 
desired for temperatures other than 28°C, corrections can be made by use 
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of the temperature coefficients discussed in Section 3.8. The evaluation 
of changes in solar cell response due to reduced light transmission in the 
cover slide materials will be covered in the next section. 

6.2 General Procedure, Absorbed Dose 

To use the cover glass darkening data previously presented, a 
procedure is necessary to evaluate the absorbed dose produced by the 
various radiation components of the space environment. This can be done 
by use of the data developed in Section 4.0. The procedure is similar 
to that used for equivalent fluence, with the exception that the absorbed 
dose is a point function and therefore varies with depth in the cover 
material. To calculate the absorbed dose at a particular depth in the 
cover materials, the following expression is used: 

OO 

Dose(d) =^0(>E ) - $(>E + AE)] • I(E,d) (6.2.1) 

E=0 

where Dose(d) = the absorbed dose in the cover material at a depth d 

I(E,d) = the absorbed dose per unit fluence for isotropic 
space radiation particles of energy E at depth d 
in the shielding material. Figs 4.2 and 4.5, Table 
4.2 and 4.5 

The absorbed dose must be calculated at several depths in the cover mate- 
rial and the electron and proton portions of the environment must be 
summed to determine the dose-depth profile. The necessity of including 
contributions from back radiations must also be considered. In practice, 
the dose deposited will decrease greatly with increasing depth into the 
cover materials. The greater dose near the surface is due largely to low- 
energy trapped protons, and contributes little to the average dose 
deposited in the cover materials. Because of the current uncertainties 
in evaluating cover material transmission loss in space, there is little 
object in an extremely accurate evaluation of the surface dose. When the 
average dose deposited in the cover material is known, the degradation in 
transmission can be estimated from the data in Section 3.11. These loss 
factors may then be applied to the estimated solar cell output parameter 
values. 
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6.3 Rough Degradation Calculations 

A rough determination of the equivalent fluence can be made by 
manual calculation. The calculation follows the procedure described by 
equations (6.1.1) and (6.1.2). The examples to be shown are designed 
for simplicity rather than accuracy. The energy increments (E l5 E 2 ) 
used in these rough calculations are those commonly tabulated in circu- 
lar orbit integrations " ' The D(E,t) values used are taken from 
Tables 4.1, 4.3 and 4.4 for the mean energy value of the energy incre- 
ment. Calculations are shown for cover glass thickness of 0.0335 gm/cm 2 
(0.006 in fused silica), 0.0671 gm/cm 2 (0.012 in fused silica), and 
0.1675 gm/cm 2 (0.030 in fused silica). The details of such an equivalent 
fluence calculation are shown in Table 6.2 for trapped proton radiation 
in a circular orbit at 835 km (450 n mi) altitude and 90° inclination. 

The D(E,t) values used in Table 6.2 are those based on P max and V . 

Several observations can be made regarding the calculations in 
Table 6,2. The largest contribution to the equivalent fluence for 
0.0335 gm/cm 2 (0.006 inch fused silica) shielding occurs in the flux 
increment between 4 and 6 MeV. The equivalent fluence contributions 
from protons with energies greater than 30 MeV appear to be negligible. 

The use of the D(E,t) value for 5 MeV (1.25) leads to serious equi- 
valent fluence errors in the energy increment of 4 to 6 MeV because D ( E , c ) 
changes very rapidly with energy in this region. The equivalent 1 MeV 
electron fluence calculated for 0.0335 gm/cm 2 shielding by this rough 
method is 6.81E13*e/cm 2 -yr. A similar detailed machine calculation 
(to be discussed) employing much smaller energy increments yielded an 
equivalent 1 MeV electron fluence of 5.57E13 e/cm 2 -yr. This difference 
is entirely due to the use of smaller energy increments in a machine 
calculation. The accuracy of the manual calculation can be improved 
by this procedure, but additional values of <f(Ei) must, be obtained by 
interpolation. It should be noted that this equivalent fluence is the 
front radiation contribution only. The back contribution must be calculated 


*Throughtout this section, the floating point notation will be used to 
represent exponential quantities. 6 .81 El 3 = 6.81 x 10 13 
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Table 6.2 Manual Calculation of Equivalent Fluence (Trapped Protons) (P 
Circular Orbit 835 km (450 n mi), 90 Degree Inclination 
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separately and added to the front contribution. It should be noted that 
the particle flux is not reduced by a factor of V 2 to allow for the 
assumed infinite rear shielding. This factor is allowed for in the 
determination of D(E,t) . In Table 6.3, the above calculation is repeated 
using D(E,t) values based on I . The procedure and problems are 
identical to those previously discussed. 

In Tables 6.4 and 6.5, the trapped proton contributions of equivalent 
1 MeV electron fluence are calculated by manual methods for a circular 
earth orbit of 4650 km (2500 n mi) altitude and 90° inclination. Such an 
orbit penetrates the region with the most intense trapped proton flux. 

As in the previous case, the major equivalent fluence calculations occur 
in the lower proton energy increments, and protons of energies greater 
than 30 MeV can be ignored without significant error. In this and the 
previous example, proton energies below 4 MeV have been omitted because 
they are "cut off" by the lightest shielding considered [i.e., D(<4,6) = 0 ]. 
The calculated equivalent fluences for this orbit are approximately one 
thousand times greater than found at the lower altitude previously con- 
sidered. In the circular orbit of 4650 km (2500 n mi) altitude, 90° 
inclination, there is also a relatively high flux of trapped electrons. 

The rough evaluation of this contribution to the equivalent fluence is 
shown in Table 6.6. The values of $(<Ej) shown are taken from integrated 

, f. C 

orbit tables from maps AE4 and AE5. ' The values of D(E,t) are taken from 
Table 4.1. The calculation procedure for trapped electrons is exactly the 
same as that for trapped protons, with the exception that one equivalent 
fluence value will describe the variation of the solar cell parameters 
I an d V,. As in the case of the trapped proton evaluations, the 

major equivalent fluence contributions occur in a few lower energy incre- 
ments. For cover glass shielding of 0.0335 gm/cm 2 (0.006 in fused silica), 
an equivalent fluence of 3.14E13 equivalent 1 MeV el ectrons/cm 2 -yr is deter- 
mined by these rough methods. A detailed machine calculation of this 
value indicates 1.99E13 equivalent 1 MeV electrons/cm 2 -yr. Although this 
fluence is large enough to produce significant solar cell degradation, if 
considered separately.it is only one-thousandth of the previously calculated 
trapped proton equivalent fluence contribution for this orbit. On this 
basis, it is reasonable to ignore the trapped electron contribution to 
equivalent fluence in this orbit [4650 km (2500 n mi), 90°]. 
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Table 6.3. Manual Calculation of Equivalent Fluence (Trapped Protons) (I 
Circular Orbit 835 km (450 n mi) 90 Degree Inclination 
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Table 6.4 Manual Calculation of Equivalent Fluence (Trapped Protons) (I 
Circular Orbit 4650 km (2500 n mi), 90 Degree Inclination 
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Table 6.5 Manual Calculation of Equivalent Fluence (Trapped Protons) (P , 
Circular Orbit 4650 km (2500 n mi), 90 Degree Inclination ' d ' 
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Equivalent 1 MeV electron fluence, P_, v ,V . 

uiqX OC 

(Machine Calculation) e/ cm -year 4.35E16 1.27E16 2.09E15 



Table 6.6. Manual Calculation of Equivalent Fluence (Trapped Electrons) 
Circular Orbit 4650 km (2500 n mi), 90 Degree Inclination 
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At altitudes greater than 4650 km, the trapped proton contribution 
to the equivalent fluence decreases rapidly and the trapped electron 
contribution becomes more significant. At synchronous or geostationary 
altitude, the trapped proton contribution can be neglected. An example 
of a rough calculation of equivalent fluence for three different shield- 
ing materials for synchronous (35,900 km, 19400 n mi) altitude and 0° 
inclination is shown in Table 6.7. The rough calculated equivalent 1 MeV 
electron fluence for cells with 0.0335 gm/cm 2 of shielding is 4.17E13 
e/cm 2 -yr. Detailed machine calculations of this quantity indicate 2.54E13 
equivalent 1 MeV electrons/cm 2 -yr. The reason for the higher value found 
by the rough manual calculation is again related to the size of energy 
increments in the lower energy range and the rapidly changing values of 
D(E,t) in these ranges. 

The calculation of absorbed dose in shielding materials is very 
similar in outline to the equivalent fluence calculation. This procedure 
is described mathematically in equation (6.2.1). In Table 6.8 a rough 
calculation of cover material dose is shown for three different depths 
in a fused silica cover glass. The rough calculated absorbed doses are 
quite close to the values calculated by detailed machine methods. The 
same procedure can be used for calculation of the absorbed dose deposited 
by trapped protons. The I(E,t) values in Figures 4.2 and 4.5 and Tables 
4.2 and 4.5 may be used for this purpose. Although the dose contributed 
by trapped protons is often very high in the surface layers of shielding, 
this is usually not a significant contribution to the average absorbed 
dose in the shielding. 

6.4 Computer Calculated Equivalent Fluence 

The aforementioned rough calculations can be improved in accuracy 
and speed with the aid of computer processing. Although the quantity 
computed is exactly the same -as before, the selection of difference flux 
and the corresponding damage coefficient can be programmed to achieve 
higher accuracy and more consistent results. The increased accuracy of 
calculated fluence is achieved mainly by use of finer energy increments 
for a given environment. 
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Table 6.7. Manual Calculation of Equivalent Fluence (Trapped Electrons) 
Circular Orbit 35,900 km (19,400 n mi), 0 Degree Inclination 
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Table 6.8. Manual Calculation of Absorbed Dose (Trapped Electrons) 

Circular Orbit 35,900 km (19,400 n mi), 0 Degree Inclination 
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♦Reference 6.4 
















For circular orbits around the earth, Vette, et al.^*"*”®*® have 
time integrated both electron and proton environments for convenient 
energy ranges, various altitudes, and for inclinations of 0°, 30°, and 
90°. The average daily fluences are tabulated in references 6.1 through 
6.9. Electron environment is taken from reference 6.5 in which both inner 
zone (AE-5) and outer zone (AE-4) electron models are integrated and 
tabulated. 

For the trapped proton environment, three maps are utilized: AP5 

(reference 6.1) for the energies from 0.4 to 4 MeV, AP6 (reference 6.2) 
for the energies from 4 to 30 MeV, and AP7 (reference 6.3) for the ener- 
gies greater than 50 MeV. For a given altitude and inclination, an 
integral spectrum of the AP7 was extrapolated back to 30 MeV and the 
intensity was normalized to the AP6 in order to eliminate a discontinuity 
there. Similarly, the intensity at 4 MeV from the AP5 was normalized 
to that of the AP6. The API data (reference 6.8) for the energies from 
30 to 50 MeV was not utilized because it is obsolete data, with the spec- 
trum not readily available and the energy interval covered being too 
small on a logarithmic scale to add any significant information. 

Although lower energy proton fluxes are more damaging than higher 
energy proton fluxes, most of the low energy component described in the 
AP5 (reference 6.3) is eliminated by the cover glass. Thus, the most 
important portion of the energy spectrum from the standpoint of solar 
cell damage is in the neighborhood of a few MeV, and hence, a normaliza- 
tion to the AP6 is of physical and practical significance. 

The assessment of solar-flare proton effects is complicated by 
several problems. They are: 

a. The unpredictable nature of future solar flare proton fluxes 
and energy spectrums. 

b. The undefinable nature of geomagnetic cutoff energy during 
a flare event, and hence, the evaluation of the near earth 
flare environment. 

c. The uncertainty in the isotropy of flare fluxes. 
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A few assumptions and limitations must be imposed on the calculation 
of equivalent fluence to overcome the above problems. As for the first 
problem, the equivalent fluence is calculated from two alternate methods 
of estimating the flare environment. The first estimated future flare 
environment is based on reference 6.9 and is conmonly used in connection 
with many military satellite systems. An alternate environment estimate 
is determined from a recently developed model. Thus, the reader can 
have a choice of selecting an appropriate environment to fit the circum- 
stances on hand. The environment in reference 6.9 leads to a worst case 
degradation estimate, whereas that described by the proposed model is a 
probable estimate. 

The magnetic cutoff energy varies with both altitude and latitude 
even during quiescent periods, and hence becomes time dependent for space- 
craft moving with respect to the earth. Further complications are caused 
by the plasma disturbance and magnetic field regime sweeping through the 
earth, the magnitude of which depends in part on the size and location 
of flares on the solar disk. Therefore, it is impossible to generalize 
all these conditions; however, there are two (2) distinct cases in which 
certain assumptions are valid as previously discussed: 

a. at high altitude and latitude, the geomagnetic field 
makes almost negligible contribution to the cutoff 
phenomena during the storm, and 

b. at very low altitude and latitude, the Stormer's 
cutoff approximation may prevail. 

The damage coefficient for omnidirectional flux is exclusively 
utilized with the followinq understanding: 

a. If the solar flare proton flux is omnidirectional through- 
out the event, the equivalent fluence computed with use 

of the omnidirectional damage coefficient described in 
section 4.0 will not result in any error due to the 
directionality of proton flux. 

b. If the flux is unidirectional throughout the event, 
though such an event is very rare, the computed equivalent 
fluence will be in error by a factor of two. 
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Therefore, the uncertainty in flux directionality can be removed 
by the use of omnidirectional damage coefficient with the provision 
that the estimate can be very reasonable for most of the events with 
a very small probability of a factor of two underestimate. 

The annual equivalent fluences thus calculated by computer are 
tabulated in Tables 6.9 through 6.30. The contents of these tables 
are summarized in Table 6.31 to assist the readers in finding a table 
for appropriate environments and parameters. 

The damage produced by back radiation is, for the first order 
assumption, regarded as the same in nature and magnitude as that pro- 
duced by the front radiation. In this context, an equivalent fluence 
attributable to the back radiation can be added to the front contribu- 
tion by estimating an effective shield thickness of back shielding. 

This bold assumption is not valid when higher order effects are con- 
sidered. If a composite back shielding material is similar to the 
front cover glass, a correction factor for a stopping power due to 
atomic number Z is small (proportional to Z) and only a density correc- 
tion is required for the estimate. This is done by shifting a curve 
of equivalent fluence vs. cover glass thickness by a density factor. 

If the Z is vastly different, the equivalent fluence should be recom- 
puted based on effective damage coefficient for the new shielding 
material . 

6.5 Solar Array Degradation 

The process of calculating equivalent 1 MeV electron fluence reduces 
the space radiation environment to a laboratory electron environment for 
which solar cell degradation has been evaluated. When the damage equiva- 
lent fluence is known, the estimation of solar array degradation can be 
completed. Up to this stage in the process of estimation of array degra- 
dation, the variables of base resistivity and cell thickness have not 
entered into the calculation. 

The necessary solar cell data and the equivalent fluence value allow 
the estimation of solar cell output parameters through use of the data 
in Figures 3.7 through 3.14. 
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Table 6.10. Annual Equivalent 1 MeV Electron Fluence Due to Trapped Electrons, Circular Orbits, 
Inclination 30 Degree, Infinite Backshielding Assumed 




















^4 

o o 


—5 cm 

! 

i 

eg <n cm eg rg cm 

rg cm cm 

rsi <M (VI 

eg eg eg 


eg eg cm 

cm eg cm 

CM CM -f 



0 

1 


— ■ 

1 •"* 

H 


^4 H 

^4 -4 

—4 


—4 —4 — 

—4 4 

-4 

»"4 

-4 H 

c 


c 

— 4 - 4—4 

-4 

^4 

-4 

-4 -4 - 



♦ 

♦ 

♦ 

♦ 

4 

♦ ♦ 

4 

* 

4 

♦ 

♦ * * 

4 4 

♦ 

4 

4 4 

c 

n 

C 

444 

+ 

4 

♦ 

444 



UJ 

UJ u 

UJ 

UJ 

UJ 

UJ UJ 

UJ 

UJ 

UJ 

UJ 

UJ UJ UJ 

UJ UJ UJ 

UJ UJ UJ 

UJ 

uu uj 

UJ UJ UJ 

UJ 

UJ 

UJ 

UJ UJ u 



in 

in -4 

m 

cm ro 

O' m co: e- 

O CO 

O so o 

h- CO 

rg 

m 

in M3 

eg in 

O' 

0 m cm 



O' 

r- r-.ir 



en 

o co 

QOCO 

f- O ■£; !\J 

CO 

ao 

O' <o eg 

(VI > 03 

m © oo 

«o 

^ eg 

©me- 

in 


in 

h- —4 if 



co 

eg in 

H 

h- 

CM 

M3 h 

O'je- 

i 

<M 

—4 —4 H 

-4 -4 

—4 

CM CO CO 

in 

M> 

r- 

7 , 

6. 

5. 


CO 

CM 

h —4 e 



H 

O -4 

j 

^(UN 

CO tO rO 

CO 

eg rg 

eg cm cm 

ea cm eg 

cm cm eg 

ro 

CO 

cn 

CO CO CO 

<0 

CM 

CM 

rg rg n 



o 

-4 


-4 

-4 


•4 H 

^4 

*-4 

—4 

—4 

H H H 

-4 r4 

-4 

*-* 

—4 —4 

—4 

-4 

—4 

^4—4—4 

-4 


H 

— 4 — 4 *“ 



J 

♦ 

4* 

D 

0 

C 

4 4 

4! 4 

4 

4 

♦ * ♦ 

4 4 

4 

4 

4 4 


♦ 

4- 

♦ 4- ♦ 

♦ 

♦ 

♦ 

4 ♦:< 



UJ 

■ 1 1 1 1 1 i 


c! 

UJ UJ 

UJ 

UJ 

UJ UJlUJ UJ UJ 

UJ UJ 

UJ 

UJ 

UJ UJ 

UJ 

UJ 

UJ 

UJ UJ UJ 

UJ 

UJ 

UJ 

UJ UJ u 



(30 

m -J* 

E 

rim; a 

-4 rg r*—! rg 

ro ro; m o om m 

CM 

r- 

CO 


M3 

-J* 

■» 00 ^ 

-4 

O' 

O' 

>0 O M 



>0 

o o 

E* 

m 

n* 

f0 O 00 ; 

in O' 

o o m vO o 

*|0 

CM 


4 WN 

0 

r* 

O' 

— CO «< 



-4 


-H 

m 

in 

1 

2. 

1. 

1. 

o 

m 

- • - 
co eg 

c\J ro 

ro 


r 

H 

w* 

H —4 —4 

-4 


in 

4 rg .*■ 



»-4 

O — 


CM 

CM 

co co co co 

ro rg 

eg rg «\j- eg eg 

CM 

1 

1 

rg cm cOico 

co 

CO 

(Tl pr, rO 

CO 

CO CM 

(VI (VI »V 



o 

—4 


•—4 

—4 

r4 

—4 r4 

—4 

-4 

—4 

—4 

H —4 r4 

~4 —4 

f4 

H 

—4 ^4 


—4 


*"4 —4 ^4 

—4 

—4 

-4 

—4 —4 — 1 



i 

♦ 

4- 

♦ 

4 

4 

4 4 

4 

4 

4 

4 

♦ ♦ + 


4 

4 

4414 

♦ 

4 

444 

♦ 

♦ 

<4 

4 4 4 



UJ 

UJ UJ 

UJ 

UJ 

UJ 

UJ Uj 

UJ 

UJ 

UJ 

UJ 

uj uj uj 

UJ UJ 

UJ 

UJ 

UJ UJ 

UJ 

Ul 

UJ 

UJ UJ UJ 

UJ 

UJ 

UJ 

UJ UJ u 



CM 

O' 

H 

O' 

O' 

r- 

eg co -4 

4- 

in 

m 

O' ^4 

O in n 

ico o ini cm *4 

M3 

O' CO >0 

<0 

4- 0 

-4 *0 If 



*■4 

4 1 4-j — 



s-oino 

co 

O' 

in <m o 

ro 

in 

co .o Oj'* r- 

O' 

O' O' o* 

>* 

^4 

0 

(n m 0 


CM 

X 

a 

^4 

in 

—4 

4- 

■4 

1^ 

*-4 CM 

N, 

—4 

CO 

»o in >t 

-t 4- 

in 


GO -4 

—4 

-4 

—4 

•4 —4 -4 

-4 

-4 

© 

M3 'fr e 

UJ 

























UJ 

'V 

CM 

o 


H 

CM 

<M 

ro co 

m; 

CO 

ro 

CO CO CM 

PM CM 

eg 

co 

ro co 

CO 

CO 

W 

CO (O <0 

CO 

CO 

CO 

eg cm n 

Q C 

X 

o 

H 

—4 

H 

—4 

i—4 

—4 —4 

r*4 

—4 

—4 

—4 

H —4 —4 

-4 -4 

—4 

—4 

—4 H 


—4 

— H 

—4 —4 «H 

—4 

-4 

—4 

H H- 

O 

o 

1 

♦ 

♦i 4 

4 

4 

4 4 

4 

4 

4 

4 

4 4 4 

4 4 

4 

4 

4 4 

♦ 

4 

♦ 

♦ ♦ ♦ 

♦ 

♦ 

♦ 


Ul 


UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ UJ 

UJ 

UJ 

UJ 

UJ 

UJ UJ UJ 

UJ UJ 

UJ 

UJ 

UJ UJ 

UJ 

UJ 

UJ 

UJ UJ UJ 

UJ 

UJ 

UJ 

UJ UJ u 

o 



in 

H 

-* 

m 

O' 

M> >0 

*-4 in 

in 

O' 

-4 O' 

M3 *■« 

in 

CM 

rg 4- 

CM 

CM 

4 

^ co r- 

h- 

n 

0 

00 4 0 




4- 

O' 

•4- 

co 

CO 

cm in 

vr r- 

o 

* 

CM O CM 

-4 -4 

GO 

o 

rg 4 

O' 

co 

s0 

r~ <o 

0 

^0 

CO 

4-^-4 

o 

</) 

• 

• 

• 

• 

• 

• 

• • 

• 

• 

• 

« 

• • • 

* • 

• 

• 

• • 

• 

• 

• 

• • 4 

• 

• 

• 

• • < 

o 

UJ 

M) 



M3 

CM 

O' 

eg co 

CO . CM 

eg 


l 

l 

9 

GO CO 

00 

—4 

—4 *4 

-4 

CM 

CM 

CM CM CM 

CM 

—4 

—4 

O' -0 4 1 

II 

o 

























o 

-4 







cn ro 


















z 

X 

CM 


H 


CM 

CO 

COiCO 

co 

(O 

(0(0(0 

CO CO 

CO 

CO 

CO CO 

CO 

CO 

CO 

<0 <0 co 

co 


CO 

cn cn rv 

M 

h- 

o 


— - 


—4 

r4 

H H 

— J 


—4 

-H 

—4 —4 —4 

-4 *4 

—4 

»— ♦ 

—4 —4 

—4 

—4 

-4 

-4-4-4 

—4 


-4 

-4 -4 - 



» 

♦ 

♦ 

4- 

4 

4 

4 4 

4 

+ 

4 

4 


4- ♦ 

4 

4 

4 4 

♦ 

4 

4- 

444 

♦ 

■f 

•4 

444 

« 

o 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ UJ 

UJ 

UJ 

UJ 

UJ 

UJ UJ UJ 

UJ UJ 

UJ 

UJ 

UJ UJ 

UJ 

UJ 

UJ 

UJ UJ UJ 

UJ 

UJ 

UJ 

UJ UJ u 

1- 

-J 

in 

O 

O' 

M3 

co 

rg, 

in co 

—4 

r- 

o- 

CO 

eg O' o 

rg h- 

v 0 

* 

CM M3 

e- 

0 

r- 

O in on 

0- 

—4 

nO 

O' rg or 

— « • 

UJ 

cn 

o 

m 

in 

o 

nj 

O' e- 

o 

CO 


>t H O'liO -t 

»!■ 

m 

1- O' 

m 

—4 

in 

C-» M3 co 

© 

(0 

© 

cn 0 <v 

cc in 


• 

• 

• 

• 

• 

• 

• • 

• 

• 

• 

• 

• • • 

• • 

• 

• 

• • 

• 

• 

• 

• • • 

• 

• 

• 

« • < 

UC UJ 

X <n 

H 

CM 

co 

fO 

H 

rg J- 

4 

4 

CO 


2 

2 

1 

1 

1 

—4 

H 

—4 -4 

CM 

CO 

co 

<n cn m 

CM 

CM 

-4 

— hN 

O <t 

■Si 

























aC Q 
< Z 


























-J < 


eg 

-4 

—4 

rg 

CM 

ro 

cn <n 

CO ro 

CO 

CO 

ro rn nv rn f> 

CO 

CO 

m (n 

fO 

CO 

CO 

m «i m 

CO 

(O 

CO 

cn cn rg 

3 


o 

-4 

-4 

•-4 

*4 


—* 

—4 

—4 

—4 

—4 

-4 *4j -4 -4 

-4 

—4 

—4 —4 

—4 

r4 

—4 

—4 —4 —4 

-4 

-4 

-4 

-g — — 1 

o 4- 


1 


4- 

4* 

4 

•4- 

4 4- 

4 

4 

4 

4 

+ + + 

♦ + 

4 

4 

4 4 

+ 

•f 

4 

♦ ♦ 4 

♦ 

4 

♦ 

♦ 4 ♦ 

Qt UJ 


UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ UJ 

UJ 

UJ 

UJ 

UJ 

UJ UJ UJ UJ UJ 

UJ 

UJ 

UJ UJ 

UJ 

UJ 

UJ 

UJ UJ UJ 

UJ 

UJ 

UJ 

UJ UJ UJ 

M < 


CO 


CO 


co 

in . —4 *~4 


CM 

>*■ 


co in Oi O' go 

•o 

-T 

r- 

—4 

co 


CM 0 

—4 

> 0 

CM 

g) rg h 

U 


*0 

eg 

—4 

O 

>0 

■tjinN 

00 CO 

•o 

© 

^ .t — 

o 

O 

-4 4 - 

—4 

h- 

CO 

— 

in 

© 

co 

^ CO 'O 

sc - 

REF 


H 

H 

ro 

^4 

cn 

-4 

co in 

in 

in 

** 

co ro co rg eg 

CM 

CM 

CM CM 

ro 

CO 

>r 

^ ^ -f 

cn 

eg 

CM 

-4 -4 O' 

























. j 

■ j 

“3 ► 

in 



rH 

rH rg eg 

ro 

co co 

CO; ro 

ro co 

co co co 

co co 

CO 

CO 

co co 

CO ro 

ro 

rn CO ro 

m 

co 

CO 

<0 co r*i 

CC Z 




H 

f4 

•-4 

—4 

H H 

»— 4 

—4 

—4 



r*4 —4 

—4 

—4 

—4 —4 


—4 

—4 

- 4 - 4-4 


r4 

—4 

-4 r4 — J 

o a 



4 

4 

4 

4 

4 

4 4 

4 

4 

4 

+ 

+ + + 

+ 4 

4 

4 

4 4 

-4 

■f 

♦ 

444 

4 

♦ 

♦ 

444 

LL ac 



UJ 

UJ 

UJ 

UJ 

UJ 

UJ UJ 

UJ 

LU 

UJ 

UJ 

UJ UJ UJ 

UJ UJ 

UJ 

UJ 

UJ UJ 

UJ 

UJ 

UJ 

UJ UJ UJ 

UJ 

UJ 

UJ 

UJ UJ UJ 




r- 

r- 


CO 

rO 


O' 

o in 

eg ro e- 

rg —4 

in 

rg o ro 

co >o 

—4 

rg m M3 

CO 

0 

—4 

O' — Cl 




M3 

in -4- 

0" 

O' 

h- CO -O 

in 

CM 

© 

o eg in 

CM 

GO 

■r m >r 

eg © 

© 

O O' m 

t • • 

© 

O 

CO 

in 0 in 



• 

• 

• 

• 

• 

• 

• • 

• 

• 

• 

• 

• • • 

• • 

• 

• 

• • 

• 

• 

• 

• 

• 

• 

• • * 



o 


•-» >$• 

—4 

-t r- 

33 

co 

00 

C- 

r- e- 

in >r 

(O 

ro 

CO co 

4 * in 

in 

M3 in u\ 

'T 


CO 

rg cm ^ 






























rg 

CM 

CM 

CO (O 

<0(0*0 

co 

CO 

rO 

CO CO CO 

co ro 

CO 

m >r 




>#• «t >t 


■j- 


4-4-4 




o o 

o o o 

O O OiO 

o 

O 

o o o 

© o 

o 

o o o 

0 0 

0 

000 

O 

0 

0 

O O C 




♦ 

♦ 

4 

4 

4 

4 4 

♦I 

4 

4 

4 

+ + + 

4 4 

4 

4 

♦ ♦ 

♦ 

♦ 

4- 

444 

*4 

4 

*4 

4 4 4 




UJ 

UJ 

v0 

LU 

UJ 

UJ 

UJ UJ 

UJ 

UJ 

UJ 

UJ 

UJ UJ UJ 

UJ UJ 

UJ 

UJ 

UJ UJ 

UJ 

UJ 

UJ 

UJ UJ UJ 

UJ 

UJ 

UJ 

UJ UJ U 



31 

CO 

rO 

p4 

,co 

in h a);g- 

o 


co O' *0 

CO 1-4 

CO 

M3 CM *-4 

O co 

e- 

in M' eg 

—4 

O' 

© 

M3 in cr 



hZ 

r» 

m 

m 


co ro r— 

CM 

e- 

—4 

v0 O in 

■4- -4- 

CO 

(vi o — 

CO 4- 

M) 

00 0 rg 

* 

in 

e- 

(3 H(C 

LLJ o 


— 

• 

• 

• 

• 

• 

• • 

• 

• 

• 

• 

• • • 

• • 

• 

• 

• • 

• 

• 

• 

• • • 

• 

• 

• 

* • 


UJ 1 


eg m 

co 

H 

—4 

-4 rg 

rg 

ro 

CM 


sT m in 

sy r- 

UJ 

VT 

—4 —4 

H 

—4 

—4 

-4 CM CM 

CM 

CM 

eg 

rg m n 

uj x 

a i 

























o 

3 

























> UJ 

*- : 

























UJ o 

— • ; 


eg eg 

rg 

CM 

rg 

co co <o:<o 

CO 

m 

ro co ro 

co ro 

CO 

CO CO CO 

3 

3 

3 

^ -t ^ 



■4- 

4-4-4 

SI 

4— 

— 

o o 

o o o 

o o o o o o 

o o o 

© o 

o 

o o o 

0 0 

S 

000 

OOO 

0 0 c 

o 

-J 

• 

♦ 

♦ 

4 

4 

♦ 

4 4 

4- 

4 

4 

* 

4 4 4 

4 4 

4 

4 

4 4 

♦ 

+ 

♦ 

444 

♦ 

♦ 

4 

♦ 4* 4 

-« ►- 

< 

X 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ UJ UJ 

UJ 

UJ 

UJ 

UJ UJ UJ! 

UJ UJ 

LU 

UJ 

UJ UJ 

UJ 

LU 

UJ 

UJ UJ UJ 

UJ 

LU 

LU 

Ul Ul U 



• 

o o, 

o o o 

O m O in 

o in 

o in o 

o o © 

O O o 

0 0 

0 

000 

OOOl 

© O c 

• UJ 


z 

in O 

in o © 

O rg in; 

r- 

o 

rg 

in r* o 

in o m 

o in o 

©. u 

0 

© -4 CM 

CO 


in 

O t- a. 

> 3 

| 

wr 

• 

• 

• 

• 

• 

• • 

•1 

• 

• 

• 

• • • 

• • 

» 

• 

• • 

• 

• 

• 

• • • 

• 

• 

• 

• • ( 

H* O 

a 

UJ 

i 


>- 1 

m 

>r 

o 

CO 

■4 —4 

-41 -4 

i 

_J 

CM CM 

eg cm ro 

CO >fr 

>r 

m 

in M3 


© 

O' 

m 

—4 

—4 

—4 

—4 —4 — > 

j 




6-21 


Table 6.11. Annual Equivalent 1 MeV Electron Fluence Due to Trapped Electrons, Circular Orbits, 
Inclination 60 Degree, Infinite Backshielding Assumed 
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Table 6.12. Annual Equivalent 1 MeV Electron Fluence Due to Trapped Electrons, Circular Orbits, 
Inclination 90 Degree, Infinite Backshielding Assumed 
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Table 6.13. Annual Equivalent 1 MeV Electron Fluence for J Due to Trapped Protons, Circular 
Orbits, Inclination 0 Degree, Infinite BackshieTding Assumed 
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Table 6.14. Annual Equivalent 1 MeV Electron Fluence for J Due to Trapped Protons, Circular 
Orbits, Inclination 30 Degree, Infinite Backshrelding Assumed 









EQU IV. 1 MEV ELECTRON FLUEKCE FOR JSC - CIRCULAR ORBIT, INC = 60 OEGREE 
OUE TO GEOMAG TRAPPED PROTQNS, REF. AP5. AP6 , ANO AP7. 
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Table 6.16. Annual Equivalent 1 MeV Electron Fluence for J Due to Trapped Protons, 
Circular Orbits, Inclination 90 Degree, InfiniteBackshielding Assumed 
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Table 6.17. Annual Equivalent 1 MeV Electron Fluence for V and P Due to Trapped Protons, 
Circular Orbits, Inclination 0 Degree, Infinite Backshielding Assumed 
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Table 6.18. Annual Equivalent 1 MeV Electron Fluence for V and P Due to Trappi 
Circular Orbits, Inclination 30 Degree, Infinite Backsmilding Assumed 
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Table 6.19. Annual Equivalent 1 MeV Electron Fluence for V and P Due to Trapped Protons, 
Circular Orbits, Inclination 60 Degree, Infinite 0 Backsfrreiding Assumed 



EQCIV. 1 NEV ELECTRON FLUENCE FOR VOC AND PMAX - CIRCULAR ORBIT, INC=90 DEGREE 
DUE TO GEOMAG TRAPPED PROTONS, REF. AP5, AP6, AND AP7. 


i 

t 



35E-01 

»h cmcm m m 

H H H ^ 
♦ ♦♦♦♦ 
UJ UJ UJ UJ UJ 

diMOl/iH 

! 

4 4 + : + 4 4 

UJ UJ UJ UJ UJ UJ 

co >0 o,f* cm m 
h O *H N O* 

CO CM O CD CO 
rH rH f-l 'O O 
4 4 4 4 1 

ill UJ UJ UJ UJ 
co co m m 




_■ ^ m .h in 

r4fOin,mfOH 

CM CM 

^ a 0 

0 0 


H 

4N M ■f 

^ >*• m m m -4- 

1 

CA Hiff' O 

cn cn 


O 

*H W rH H H 

w 4 H H H H 

rtrtWoo 

0 0 


; I 

4 4 4 4 4 

444444 

4444 4 

1 1 


!UJ 

UJ UJ UJ UJ UJ 

UJ UJ UJ UJ UJ UJ 

UJ UJ UJ 

UJ UJ 

UJ UJ 


1 00 

O M<1 f i 

o in vO O' o r- 

00 c- 

O' nj 

r- in 


•o 

ro n, SO O' O 

O l/i N O -O ^ 

CM N- 00 

-O O 

n >* 


j rH 

N N JIM H 

OJ 00 (NJ r- f— 

H H fO 

m h- O 

O' M5 


rH 

.4 im m m + 

4* in in m in in 

\ 

cn cm O cm cn 

rH ^H 


o 

r— 1 rH rH H »-H 

H |U rH r— 4 rH H 

rH H rH ! rH O O 

0 0 


1 

4 4 4 4 4 

4 4 4 4 4 4 

+ *■ * * * 1 

I 1 


LU 

UJ UJ UJ UJ UJ 

UJ UJ UJ UJ UJ UJ 

UJ UJ UJ 

UJ UJ UJ 

UJ UJ 


CM 

«o O N N m 

■f sl'tOMH 

-oo 

CO ^ 

*0 n 


rH 

rT) f<\ ^ ^ >0 

4 m min o « 

CM rH CM [CM 0" cm 

m cm 

Csl 

31 

O 

H 

f<i »4 -t H 

4 fAH 

1 

l 

cn m h.h n ^ 

i 

cn cm 


jOJ 

H (Ni n m ^ 

^ n n >o 'O ^n 

m ^ cm 

0 m 0 

rH rH 

X 

o 

pJ H H H H 

H H H r4 H H 

H rH H H O O 

0 0 

o 

: I 

4 4:4 4 4 

4 4 4 4 4 4 

•f ♦ + 

4 + 4 

4 4 


UJ 

iii m in u i ni 

UJ UJ UJUJ UJ UJ 

UJ UJ UJ 

UJ UJ UJ 

UJ UJ 

•> 

rH 

>o <4* m *4* o 

CM in *O N O <M 

m cn n> 

O' m m 

O' ■* 

UI 

|N* 

m o\in o o 

O >* O CM r-4 

- - (M 

O' rH 

rH CM 

oO 

• 


• • • 

• • • 

• • 

UJ 

j-O 

o m oj co m 

<7 

I 

I 

6 

C 

A 

- N -O 

>0 ^ CVJ 

in cn 

Z 






-3C 

o 














X 

CM 


m in o sO o 

in n m 

rH 00 m 

4- 4- 

►— 

o 

*-H H rH «— » H 

^4 r4 H H H r4 

H H H 

-4 0 0 

0 0 


1 

♦ 4 4 4 4- 

4 4 4 4 4 4 

■f * +i+ ♦ + 

♦ + 

Q 

uu 

UJ UJ UJ UJ UJ 

UJ UJ UJ UJ UJ UJ 

UJ UJ UJ UJ UI UJ 

UJ UJ 

-J 

in 

-< in ao .T r- 

oon^ ow 

m co 0 s 

rH rH N- 

CO CM 

UJ 

m 

sO m in in 

co m co 4 * cm n 

CM *H CM 

m rH %o 

(<i •? 


• 

• • • • • 

• ••••• 

• • • 

• • • 

• • 

X 

ro 

rt H in H >0 

^ ^ N H 

tn ^1 >»• 

n -h in 

CM rH 

U) 








cm 

nj in <r in 

in -O 'O N N o 

•O ~0 <T 

ro O 0 

•o >0 


O 

H H *H H H 

rH rH "4 rH f4 H 

H H H 

rH rH O 

0 0 


i 

4 4 + + 4 

444444 

* ♦ + 

4 4 4 

4 4 


UJ 

UJ UJ UJ UJ UJ 

LI 1 III Ui 111 111 LI 1 

UJ UJ UI 

LU UJ LU 

UJ UJ 


30 

>t nj h- m n. 

in o H N lA H 

'O 0 <t 

0 m nO 

rH CM 


■o 

CO H (M C' ^ 

OMHHHrl 

- 0-10 

m n* so 

CM CM 



m ffl H m H 

^ M S.H r4 N 

m rH CO ; CM OO CJn 

O' in 


; 

U> «/i ^ O 

rH H H ^ H 

4 * 4 4 4 

UJ UJ UJ UJ UJ 

cm O' N- •-* o 

o«t 'tHm 

*0 N- CO ao CO O' 

H H H rH rH rH 

4 4 4 4 4 4 

UJ UJ UJ UJ UJ UJ 
O 4 'O;^ fA © 

*-h co o : >o o' O' 

OOO-HrHsf 

CM CM cm cm CM CM 

444444 

UJ UJ UJ UJ UJ UI 

fA >f CO N in CM 

CM CM rH O r- N- 

* 
CM CM 
4 4 
UJ UJ 
cr 

>0 in 


o 

N h ^ h rn 

CO cn rH *-H rH CM 

p- nj 

rH O' rH 

> 0 CM 



• i 

(M n ,m m ra 
O OO O o 
4 4:4 4 4 

i 

t 

m cn mm cn cn 
o o o o o o 

4 4 4 4 4 4 

(*1 >t ^ ^ -f ^ 
O 0 0 O 0 0 
+♦+.+++ 

4’ 4- 
0 0 

4 4 



UJ UJ LU UJ UJ 

LU UJ UJ ; LU UJ UJ 

JJ UJ UJ 

LU LU LU 

JJ UJ 


X 

co sO m «—c* co 

in oo ojm .0 «-< 

v0 00! 

in CM O' 

>o cn 




03 N N'O in ^ 

CM ^4 >4i 

co cm in 

O' cn 


%»» 

• • • • • 

• ••!•••■ 

• • • 

• • • 

• * 

UJ 

O 

z> 

h- 


ru uv co h 

j 

I 

>-> nj m vr m r- 

i 

t 

l 

O' H H; 

i 

rH CM CM 

CM cn 

r-* 


cm cm cm cm rM 

rn cn m rn cn cn 

cn rn m >j- 

4- 4- 

N* 

— % 

o o o o o 

o o oo o o 

000000 

z> 0 

— J 

• 

♦ 4(4 4 4 

+++,++* 

4 4 4 

4 4 4 

4 4 

<x 

X 

LU UJ UJ UJ UJ 

JJ UJ LU LU UJ UJ 

UJ UJ UJ 

UJ LU UJ 

U UJ 


• 

o o o o o 

o o oo o o 

000 

OOO 

0 0 


i: 

inomuo 

o in o in o o 

OOO 

ON 4 

>0 00 



h >o m 

| 

1 

h h M cq ff) sj- 

• • * 

in 'O 00 

• • • 

rH rH »H 

• • 
rH rH 


6-30 


Table 6.20. Annual Equivalent 1 MeV Electron Fluence for V and P Due to Trapped Protons, 
Circular Orbits, Inclination 90 Degree, Infinite Backsinelding Assumed 
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Infinite Backshielding Assumed 
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Infinite Backshielding Assumed 
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Table 6.23. Annual Absorbed Dose Due to Trapped Electrons, Circular Orbits, Inclination 
60 Degree, Infinite Backshielding Assumed 



dSQREEC OCSE( RAD SI). CIRCULAR ORBIT, INC = 90 DEGREE 
DUE TG GECHAG TRAPPED ELECTRONS, REF. AE4 AND AE5. 
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Table 6.24. Annual Absorbed Dose Due to Trapped Electrons, Circular Orbits, Inclination 
90 Degree, Infinite Backshielding Assumed 







ABSORBED OOSEIPAD SI), CIRCULAR ORBIT, INC = 0 DEGREE 
DUE TO GEOMAG TRAPPED PROTONS, REF. AP 5 , AP6, AND AP7 
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Table 6.25. Annual Absorbed Dose Due to Trapped Protons, Circular Orbits, Inclination 0 Degree, 
Infinite Backshielding Assumed 
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Table 6.27. Annual Absorbed Dose Due to Trapped Protons, Circular Orbits, Inclination 60 Degree. 
Infinite Backshielding Assumed 
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TABLE 6.30 EQUIVALENT 1 MEV ELECTRON FLUENCE FOR SOLAR FLARE 
PROTONS, BASED ON TADA MODEL 
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Table 6.31. Summary of Tabulated Data in Tables 6.9 Through 6.30 


Absorbed 
Dose at 
Various 
Depths 

Table 6.21 
Table 6.22 
Table 6.23 
Table 6.24 

Table 6.25 
Table 6.26 
Table 6.27 
Table 6.28 
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The previously calculated equivalent fluence data for synchronous 
orbit will be used in the following example to illustrate the degradation 
of a solar array: 


Solar Cell 


Cover Glass 


Backshielding 


10 ohm-cm resistivity 
0.0305 cm (0.012 in.) thick 

0.015 cm (0.006 in.) thick 
0.0335 gm/cm 2 

fused silica, antireflecting 
coating, blue filter 

Infinite 


Equivalent 1 MeV Electron Fluence 


Trapped Electrons 
Trapped Protons 

Total 

Solar Cel 1 Output 


2.55E13 
0 

2.55E13 e/cm 2 -yr 
Absolute Relative 


At 0 Months 


I 

V 
P 

V 


sc 

oc 

max 

mp 


At 4 Months 


I 

V 
P 

V 


sc 

oc 

max 

mp 


At 1 Year 


I 

V 
P 

V 


sc 

oc 

max 

mp 


Equivalent Fluence, 0 

36.7 mA/cm 2 
0.552 V 

14.8 mW/cm 2 
0.433 V 

Equivalent Fluence, 8.50E12 e/cm 2 

36.2 mA/cm 2 
0.548 V 

14.7 mW/cm 2 
0.432 V 

Equivalent Fluence, 2.55E13 e/cm 2 

35.7 mA/cm 2 
0.543 V 

14.2 mW/cm 2 
0.430 V 


1.00 
1 .00 
1 .00 
1 .00 


0.986 

0.992 

0.994 

0.998 


0.975 

0.984 

0.960 

0.994 
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Absol ute 


Relative 


1 Years 

Equivalent Fluence, 7.65E13 e/cm 2 


'sc 

34.9 mA/cm 2 

0.950 

v oc 

0.535 V 

0.970 

p 

max 

13.4 mW/cm 2 

0.905 

V mp 

0.426 V 

0.985 


The effects of cover glass transmission loss due to radiation 
darkening have been omitted from this estimate. The average absorbed 
dose due to trapped electronics in this orbit is approximately 10? 
rad per year. The data in Figure 3.23 indicate that such a dose 
would cause a transmission loss of about 0.5% reflected in solar cell 
I or P , . The effects of solar flares on solar cell degradation 

S C maX 

have been omitted in this example. Adhesive darkening due to ultra- 
violet illumination has been omitted from this calculation. An 
additional factor which also must be included is the reflection loss 
due to glassing of solar cells with silicon monoxide anti-reflection 
coatings. With conventional state of the art cover glasses and 
solar cells, glassing may cause a 4 to 6% initial decrease in short 
circuit current. As discussed in Chapter 1, cells with improved 
anti -reflection coatings exhibit significantly reduced glassing or 
reflection loss. 
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CHAPTER 7 


7.0 FLIGHT DATA 

Considering the number of satellites in orbit, there is a very 
limited amount of array performance data available. A very large per- 
centage of satellites are placed in low altitude earth orbits (less 
than 400 km) which experience very little degradation from trapped 
radiation. Of the satellites in orbits above 400 km, very few are 
instrumented to sense and telemeter electrical quantities which directly 
reflect basic solar-cell performance parameters. Flight data from prior 
to 1963 often involve p-n solar cells or n-p solar cells with much lower 
efficiency and poorer radiation hardness than that of currently available 
commercial cells. There are two sources of more recent solar cell per- 
formance flight data. These include a few solar cell flight experiments 
specifically designed to evaluate solar cell degradation. Such experi- 
ments have been flown on the ATS-1, 7 ' 1 LES-6 7 ' 2 and ATS-5 7 ' 3 satellites. 
An alternative approach to the evaluation of solar cells in a space 
radiation environment is a general analysis of an entire satellite solar 

array. An example of this type of analysis was reported for IDSCS 
7 4 

satellites. ’ Most of the available flight data are from satellites in 
synchronous orbit, however, a limited amount of data are available for 
lower earth orbits. 

The data in Table 7.1, relating to solar array performance in a 
synchronous orbit, were collected by L. A. Gibson of the Aerospace 
Corporation. 7 ’5 All the solar cells used in these satellites have 
10 ohm-cm base resistivity and cover glass shielding varying from 
0.015 cm (0.006 in) microsheet to 0.076 cm (0.030 in) fused silica. No 
information was reported regarding cell thickness or backshielding. The 
reported degradation in power in most cases is between 2 to 6 percent 
after one year. The power loss estimated in Section 6.6 on the basis of 
trapped electrons alone was 4 percent per year for cells with 0.015 cm 
microsheet shielding. Considering that solar flare proton contributions 
to equivalent fluence and cover material darkening were omitted in this 
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calculation, there is reasonably good agreement between satellite perfor- 
mance and the predictions. The omission of solar flare equivalent flu- 

ence contributions appears justified, as flare activity was relatively 

7 5 

low during the time period of the subject flight data. * The poor 
performance of IDSCS satellite solar arrays resulted from excessive ultra- 
violet degradation transmission loss in the cover glass adhesive due to 
the use of an improper primer. 


Although degradations due to solar flares are often estimated and 
projected over long satellite missions, the flare events are discrete 
and their effects occur as rather abrupt degradations. An excellent 
example of this behavior is shown in Figure 7.1 for two satellites in 
synchronous orbits during the flare events of August 1972. The analysis 
was provided by H. Riess of TRW.^* 8 The solid line in Figure 7.1 is 
based on solar cell degradation predictions based on trapped electrons 
at synchronous altitude. The data indicate that the flares produced 
an abrupt 2% loss in maximum array current ( i . e . , short circuit current) 
in both satellites. It also can be observed that 5 months after the 
flare, the flight 3 array current had recovered to within nearly 1 per- 
cent of the value predicted without solar flares. This indicates that 
considerable annealing of flare radiation damage occurs after termination 
of the event. The behavior of a group of test cells on the ATS-5 satellite, 
during the August 1972 flare events, is shown in Figure 7.2. These cells 
degraded 0.2 percent in V oc , 2.3 percent in I sc , and 3.5 percent in P max 
during the above period. 7-7 

A limited amount of flight data are also available from satellite 
solar arrays operated in circular orbits at lower altitudes. Data from 
several such satellites are tabulated in Table In these 

orbits the equivalent fluence contribution is mainly from trapped pro- 
tons. Two of the orbits, ERS 6 and Explorer 38, pass through regions 
in which the trapped proton flux is extremely large. The equivalent 
1 MeV electron fluences are shown for each orbit. These were determined 
from Tables 6.16 and 6.19 with the assumption that: only trapped proton 

radiation contributed to the equivalent fluence, that infinite back- 
shielding existed, and cover material darkening losses were negligible. 
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SOLAR ARRAY OUTPUT 
SOLAR ARRAY SUN ANGLE (0) = 
1 AU 
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TIME IN ORBIT, MONTHS 

Figure 7.1 Performance of Two Satellite Solar Arrays in Synchronous 
Orbit During the August 1972 Solar Flares 




(MWJ H3MOd wnwixvw 
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Figure 7.2 Degradation of Solar Cell Maximum Power Versus Time in 
Synchronous Orbit, ATS-5 Experimental Cells 7.7 



TABLE 7.2 SOLAR CELL ARRAY DEGRADATION, VARIOUS CIRCULAR ORBITS 
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The equivalent fluence values were used to estimate solar cell parameter 
changes from data in Figures 3.7 through 3.14. The predicted changes 
are shown in Table 7.2 along with observed parameter changes from flight 
data. The predicted degradations are in reasonable agreement with ob- 
served values. 

The results of experiments on ERS 6 included several observations 
which have important consequences in array degradation predictions. 7 ^ 0 
The cells of this satellite were observed to degrade in short circuit 
current at a rate of 5.5 ± 0.2 mA/cm 2 -decade. This value compares well 
with those reported in Section 3.3 for laboratory proton irradiations 
in the 10 MeV energy range. Since the above rates are higher than those 
normally found for 1 MeV electron irradiations, it would be more accurate 
to use experimental proton degradation data for proton dominated orbits 
if such data were available. It was also observed that cells with 
adhesively-attached cover glass shields degraded at the same rate as 
those with mechanically-attached (no adhesive) shielding. It was con- 
cluded that adhesive darkening effects were less than the experimental 
error or negligible. The data also indicated that transmission loss in 
cover glass is not an important factor in array degradation. 

In general , many factors can make array performance degrade to 
outputs lower than predicted. One of the most common causes is the 
effect of low energy protons on small unshielded areas of the solar cell. 
Another problem which has caused serious difficulties is the failure of 
cell interconnections, particularly in the case of solderless solar 
cells. 

Gibson concluded that an initial 2 percent transmission loss, due 

to ultraviolet darkening of cover glass adhesive when added to the solar 

cell degradations, would account for the power degradation observed on 

7 5 

most satellites reported in Table 7.1 In the LES-6 experiments, 

Sarles and Stanley concluded that an initial 4 to 10 percent loss was 

7 2 

caused by cover slide material darkening. 
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Areal Density Fused Silica Microsheet Aluminum 

2.2 g/cm3 2.5 g/cn?3 2.7 g/cro 3 



Table A Shielding Thickness Conversions 










































TABLE B CONSTANTS, PROPERTIES, AND VALUES 


SILICON 


Atomic Weight 
Density 

Crystal Structure 

Lattice Constant 

Atomic Radius 

Atomic Density 

Energy Gap@ 300K 

Energy Gap @ 0 K 

Electron Mobility (intrinsic) 

@ 300K, y p 

Hole Mobility (intrinsic) 

G> 300K, u p 

Electron Diffusion Constant 
(intrinsic) @ 300K, D p 

Hole Diffusion Constant 
(intrinsic) @ 300K, D p 

n . @ 300K 

Dielectric Constant 

Specific Heat, Cp @ 300K 

Thermal Conductivity (a 300K 

Coefficient of Thermal Expansion, 

AL 

LAT 

Debye Temperature 

Activation Energy, Self Diffusion 

Energy of Ionization 

Energy of Sublimation 

Elastic Moduli 


28.09 

2.33 (g/cm 3 ) 

Diamond, 8 atoms/unit cell 

5.43xl0 -10 m, 5.43 (A) 
n o 

1.18x10 m, 1.18 (A) 

5.00x1 0 22 (cm -3 ) 

1 .78x1 0 -1 9 (J), 1.11 (eV) 

1.91xl0 -19 (J), 1.21 (eV) 

1350 (cm 2 /V s) 

480 (cm 2 /V s) 

35 (cm 2 /s) 

12 (cm 2 /s) 

I. 5xl0 10 (cm -3 ) 

II. 7 

0.7 (J/g K) 

1 .5 (W/cm K) 

2.5xl0 -6 (K -1 ) 

658 (K) 

7.7xl0 -19 (J), 4.8 (eV) 
5.76xl0 -19 (J), 3.6 (eV) 
7.80xl0 -19 (J), 4.9 (eV) 

1. 674x1 0 11 (N/m 2 ) 

0. 652x1 0 11 (N/m 2 ) 
0.796X10 11 (N/m 2 ) 


Index of Refraction 
Mohs’ Hardness 


3. 5-4.0 
7 
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TABLE B-continued 

SILICON-continued 

Solar Absorptance 

0.8 

Hemispherical Emittance 

0.3 

QUARTZ GLASS (FUSED SILICA) 
Molecular Weight 

60.8 

Density 

2.2 (g/cm 3 ) 

Energy Gap 

12.8xl0" 19 (J) , ^ 8 (eV) 

Dielectric Constant 

3. 5-3. 9 

Index of Refraction 

1.46-1.51 

Specific Heat, C 

l (J/g K) 

r 

Thermal Conductivity 

0.014 (W/cm K) 

Coefficient of Thermal 
Expansion, , AL 

0.55x1 0" 6 (K" 1 ) 

LAT 

Mohs* Hardness 

4.9 

Young’s Modulus 

7 . 1 6x1 0 1 0 (N/m 2 ) 

Rigidity Modulus 

3.10xl0 10 (N/m 2 ) 

Poisson’s Ratio 

0.16 

Solar Absorptance 

0.01 

Hemispherical Emittance 

0.78 

Solar Absorptance (on array) 

0.75-0.85 

Hemispherical Emittance (on array) 0.78-0.80 

SILICONE ELASTOMERS (TYPICAL) 
Density 

1.1 (g/cm 3 ) 

Index of Refraction 

1 .41 

Coefficient of Thermal 
Expansion, AL 

300x1 0' 6 (K -1 ) 

LAT 

Thermal Conductivity 0 300K 

.0017 (W/cm K) 

Specific Heat 0 300K 

1.0 (J/g K) 

Bond Thickness Between Cover 
Glass and Solar Cell 

75-150 (pm) 

A -3 



TABLE B-continued 
' CONSTANTS 

Boltzmann's Constant, k 
Planck' s Constant, h 
Speed of Light, c 
Electron Charge, e 
Permittivity of Free Space, e Q 
Permiability of Free Space, 
Electron Rest Mass, m g 
Proton Rest Mass, 

Avagadro's Number 


1.38xl0" 23 (J/K), 8.62x10 
6.63x1 0" 34 (J s) 

2.998xl0 8 (m/s) 

1. 602x1 0" 1 9 (C) 

8. 86x1 0" 1 2 (F/m) 

12.6x1 0" 7 (H/m) 

9:ilxl0" 31 (kg) 

1 .67x1 0~ 27 (kg) 

6. 022x1 0 23 (g-mole" 1 ) 
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(eV/K) 


SILICON SOLAR CELL DATA 

Active Area of 2 cm by 
Solar Cell 

Series Resistance, R g 
Shunt Resistance, R $h 


2 cm ? 

3.8 (cm; 

0.2-0. 5 (ohm) 

> 1000 (ohm) 
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